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PRODUCTION AUTOMATION IS THE MOST IMPORTANT LINK 
IN TECHNICAL PROGRESS 


(With reference to the Plenary Session of the CPSU Central Committee) 


A Plenary Session of the Central Committee of the Communist Party of the Soviet Union is being held on 
June 24, 1959. The Session will discuss a range of questions closely connected with automation. 


A radical increase in labor productivity required for establishing the material and technical basis of a 
communist society and the improvement of working conditions demand the transition to a qualitatively new stage 
of production organization, to group mechanization and automation of all the branches of national economy, 
Automation provides for a substantial increase in the productivity of many technological processes, the improve- 
ment in the quality of production and for better working conditions, 


Thus, for instance, automation increased the productivity of soda ovens by 22% and that of the fractionating 
columns in the manufacture of rubber by 50%, The automation of the shop for regenerating acetic acid at the 
Vladimir chemical] plant increased the productivity of labor by 45% and lowered the cost of production by 15%, 


moreover all the expenditure incurred repaid itself in 2.5 months, Automation of a coal-cleaning plant increased 
production per worker by a factor of 2,5, 


Automation considerably curtails relative capital investment. Increasing the output of blast furnaces by 


means of their automation amounts to one tenth or one twelfth of the cost of building new furnaces to produce 
the same additional output, 


The introduction of new technological processes especially insuch branches of industry as the atomic, 
chemical or metallurgical is only possible on the basis of highly mechanized and automatic production groups. 


It is only due to automatic control that the artificial earth satellites and cosmic rockets can be guided along the 
required trajectory with a high precision, 


Our country has achieved considerable success in the automation of several branches of our national 
economy: power production, oil processing and metallurgy. 


At the present time the level of development of our science and technology is such that full possibilities 
exist of introducing automation on a wide scale in all the branches of our national economy and to start the 


transformation of technological processes on the basis of their complete automation, For this purpose, we have 
the following prerequisites, 


Progress in the theory of automatic control and regulating equipment provides the possibility of introducing 


automatic control of production under the most advantageous technological conditions on the basis of self- 
adjusting systems. 


Many scientific research and experimental development organizations have developed and tested models 
embodying new technological ideas in the field of automation, 


Our industry is changing over to continuous, conveyor belt methods of production, 


In other words, our country has accumulated sufficient scientific and technical experience in the sphere of 
automation. 
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On the whole, however, automation in our country lags behind both our requirements and some of the 
capitalist countries, especially the USA, As yet, we introduce automation in separate production processes 
instead of making the whole production cycle automatic, this lowers the efficiency of automation. Such labor- 
consuming operations as loading and unloading, finishing processes and transportation lack not only automation 
but often lack simple mechanization as well. As the result of this in some industries many more workers are 
employed in auxiliary operations than in basic production (for instance, in superphosphate production, bakeries, 
in processing seeds and foodstuffs, etc.). The low level of mechanization of the finishing processes in rolling 
mills and the lack of adaptability in the existing equipment prevents their automation, although 70% of the 
workers in rolling mills are employed in this work. 


The main reasons for the lag in automation of our industry are the following. 


Automation is often carried out without due consideration of the technological processes involved and is 
often applied to processes which are not adaptable for automation. 


When designing automatic installations often appropriate calculations of the expected efficiency are not 
made, 


Another reason is the insufficient number and size of instrument- making plants which produce automation 
equipment, the lack of standardization of components produced by them, the iack of many important types of 
equipment and, obsolescence of a large part of the instruments produced, which complicates the design and in- 
troduction of new automatic systems, The lack in automation equipment restricts the initiative of industrial 
establishments which with an abundance of technical facilities would have considerably assisted automation. 


Bad organization of automation work leads to an increase in the time required for the design, installation 
and assimilation of the automation equipment, to a lack of co-ordination in the development of new tech- 
nological processes and the planning of new plant automation, As the result of this the forces and facilities 
available are not fully utilized. The development of new installations is often duplicated by several organiza- 
tions whereas some of the problems are not dealt with by anyone. 


There is a lack of qualified personnel in the sphere of scientific research, design of instruments, planning, 
installation and maintenance of automatic plants, The training of specialists in automation is at present pro- 
ceeding at such a slow pace that the discrepancy between the number of specialists required and available instead 
of decreasing is actually increasing from year to year, The number and size of scientific research organizations 
obviously does not meet the requirements in this sphere. Insufficient experimental facilities in the Institutes, 
OKBs, and factories leads to delays in research work and especially in its application. As the result of this 
many technical ideas of great economic importance which originated in the USSR were first applied abroad. 


In order to overcome the above defects it is necessary greatly to expand the industry producing equipment 
for automation and mechanization. The growth of this industry must proceed at a pace which does not hold back 
the development of other industries. The increase in the production of automation and mechanization equipment 
can and must be carried out by means of a more efficient utilization of existing production capacity through 
better organization, cooperation, standardization and correct distribution of the various types of equipment, 
which can raise production by a factor of 2-3 and also by the construction of a new plant. 


The instrument-making industry must be directed toward the establishment of a single state system of 
automatic equipment which would include various technical facilities required for the solution of modern com- 
plex problems of industrial automation, This system must be organized on the unit construction principle. The 
relatively small number of types of equipment in this case will provide the possibility of organizing their mass 
production, decrease their cost and raise their reliability, By means of suitable staffing it will be possible to 
satisfy all the requirements in automatic checking and control equipment. 


It is necessary to reorganize the system of designing, manning and installation of automatic equipment, 
concentrating all the efforts of the organizations involved on the solution of the most important nation-wide 
problems of automation, and establish priority in this work according to its economic importance. Moreover, 
work should be organized in such a manner that technologists and designers of the machines take part in it side 
by side with automation experts. The latter condition is a prerequisite for efficient group automation installations. 
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It is also necessary to speed up the organization of model automatic experimental installations, shops, and 
plants in the basic industries as envisaged by the decisions of the 20th Congress of the CPSU, The most qualified 
scientists, engineers, automation specialists, technologists and designers should take part in the work on these in- 
stallations. They should have priority in instruments, machines, cables, etc. 


In order to ensure technical progress in automation, it is necessary to speed up considerably the organization 
of scientific research and designing agencies which are developing new principles for the functioning of automatic 
control systems, new automatic equipment and new automation systems. In the coming seven years, the capacity 
of the scientific research organizations in the sphere of automation must be greatly increased, Their experi- 
mental and production facilities must be strengthened considerably and experimental workshops and laboratories 
set up in every instrument- making plant, 


Organizational measures must be taken in order to prepare a state plan for introducing automation and 
mechanization in all branches of our national economy, 


It is necessary to expand the training of automation personnel, to train automation specialists in all the 
technological institutes, establish faculties of automation and remote control in polytechnical institutes, to train 
scientific personnel for automation in the universities, extend correspondence courses on automation organize 
refresher courses for engineers with the object of raising their qualifications in the sphere of automation, greatly 
to expand the network of technical schools which train automation specialists, 


Side by side with group automation it is necessary to develop small scale automation and mechanization 
of separate installations and processes, Such automation based on the initiative of the plant management can 
be carried out without delay and on a wide scale providing important economic benefits, 


In evaluating the progress of a plant, the amount of new technique it has introduced and in particular the 
amount of automation and mechanization should be taken into account, It is necessary to establish in all the 
industries a system of awards which would provide considerable material benefits to the entire personnel, from 
the director to an ordinary worker, of a plant which introduces automation, 


The fulfillment of above measures will speed up the reorganization of all the main industries on the basis 


of automation according to the decisions of the CPSU 20th Congress on the establishment of material and technical 


foundations for a communist society, 


Automation will be most effective when the whole of the industry is reorganized, It can be compared with 
electrification which has created a new power basis in production and thus changed the technological processes, 
machines and methods of production. Automation will become the main controlling factor in production, it will 
introduce radical changes in industry leading from unit and group automation (to “automatic machine systems” 
according to Marx) to shop and plant automation which will in the long run ensure the highest possible labor 
productivity, The creation of automatically controlled production is the characteristic development of technology 
in all the branches of our industry for the very near future, 


Widespread automation of production in the process of building communism will radically change the 
nature of labor, will raise the cultural and technical level of the workers and peasants and will contribute sub- 
stantially to the elimination of differences between intellectual and physical work, between the town and the 
country. Automation will also reduce the working day, and provide more leisure for the working people in 
improving their scientific and technical knowledge and developing their talents. 


The Soviet people welcome the words of the CPSU 20th Congress resolutions “Our coveted aim is near at 
hand. We shall have to pass through a decisive stage of peaceful economic competition with capitalism and 
win the contest in the shortest possible time. We have all that is necessary to win this prize. When we solve 
these problems and open up new horizons it will be easier to advance still further. For the sake of the great aim 
of building communism we should work as hard as possible," 
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SOME SIMPLIFIED STABILITY CRITERIA FOR NONLINEAR 
CONTROLLED SYSTEMS 

A. K. Bedel*baev 

(Alma- Ata) 


Certain simplified stability criteria which do not require preliminary 
canonization of their initial equations are given for nonlinear controlled 


systems. 


1. Plotting Lyapunov's V-function 
One of the possible ways of plotting Lyapunov’s V~function is given in [1] in a quadratic form: 


Q= > (1.1) 


k=1 


its derivative Q with respect to time, found by means of linear differential equations with constant coefficients, 


te = >) Perm (#=1,2,...,2), (1.2) 
k=1 
has the form: 


This method is extended in the present article to tne case when expression (1.1) satisfies a partial 
differential equation of a more general form: 


+ + Ponta) = —F (1.9) 


s=1 


where F (ny, Ng... » Np) is a predetermined quadratic form with constant coefficients as, = ays t 


F= » (1.4) 


s=1k=1 


which are only subject to Sylvester's inequality 
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As it will be shown later the use of function F provides the possibility of obtaining a more suitable region 
of stability for the choice of parameters for a required system's regulator, 


Let Age Ag» + « » An be the roots of a characteristic equation 


n 
(— 41)" D(X) =| Per — = 0 (1.6) 
of system (1,2) which satisfy the conditions: 
he An (1.7) 
Rer,<0 (1.8) 


Let us denote by nys(t), 1s (©) = 5, the fundamental system of solutions for Eq. (1,2) and assume that: 


(1.9) 
im] 


The definitely positive Lyapunoy's function can be written in the form [2): 


Q=\ de. (1.10) 


The calculation of this function requires as it was pointed out in [1] the knowledge of roots ),. In this 
connection it would be desirable to obtain an expression for o,, in terms of the coefficients of the initial system 
(1.2) and the given formula (1,4) instead of the roots of (1.6), For this purpose the method indicated in [1] can 


be used, 
In fact from (1.9) and (1.10) we obtain: 


a 6 


But since according to [3] 


Dig (4) Dig (—2,) 
\ ma (t) (t)de = >} 
é 


Q assumes the form (1,1) in which 04; is determined by the formulas 


= > > Gap 


a=] s=1 


Assuming 


Any Ang -.- 
n 
0 
i 
) 
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n n—t 
@apDia(d) Dip (—d) = 1,2, 2) (1.11) 
r=0 


and reasoning in the same way as in [4] these formulas can be given the final form: 


( 1)" n—1 
= 2 (i, 7 =1,2,...,n). (1.12) 


Here H is Hurwitz's last determinant 


O 
0 0 0 fares fa 


made up of the coefficients of Eq. (1.6), and Hy,r+, is a minor (including sign) of the first line and (r+ 1)-th 
column element of this determinant; it will be seen from (1.11) that quantities depend in the general case 
on the coefficients of the initial system (1.2) and form (1.4), Thus, with n= 2 they becomes 


qi = — Ary, qs = — Ay + 412 (Pr1P22 + Pr2P21) — (1.14) 


= — qs — P12 + 


The expression for af? is considerably simplified if the coefficients for (1.4) are chosen according to 
condition: 


0, s+k, 
= 
ask k, (1,15) 
where a, is an arbitrary real positive constant. Then with n= 3 the expression for aft) becomes: 
= 0, giz = (My, + + Por Bis) — Psa, 
gis = 0, gis = (Ms, + PisBos) — + 43 (Mis + psi Bi2), 
= + 42M +- 03M 13M 55, (1.16) 
qs: = = — + (2M — Bis) — 


qs = 0, = — + (M go + PasB1s) + (Mos + 
qs = + 2M + 43M 3M 3, 


= 4s, = — is — + (2M — Bin), 
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where Mj, is a minor (iricluding sign) of the i-th line and j-th column element of the characteristic determinant 
D(\) with ) =0 and n=3 and with quantities Biy equal to: 


-= Pir + Poor Big = Pir + Pas, Bas = Poo + 


In conclusion let us note that there is another approach to the selection of coefficients of (1.4), namely, 
let 


a= = fon = 0, 
dy, = fo, = foy =O, (1,17) 
an, = 0, ..., Gun = fay 


where fo, fy,...» fp, are the coefficients of the characteristic Eq, (1.6). Then Lyapunoy's V-function plotted 
for the system (1.2) in the form (1.1) will provide the stability region determined by Hurwitz's necessary and 
sufficient conditions with respect to coefficients fo, fy...» fn of Eq. (1.6). Moreover, the structure of formula 
(1,12) will become considerably simpler.* Thus formula (1.14) will now assume the form: 


= 0, qs = — — 


(0) w (1,18) 
qu = — fo, = foPurP22 — — 
qs) = — fe, = — foPuPre + + 
2. Simplified Stability Criteria of Controlled Systems 
1. Let us examine a controlled system represented by equations [5] 
n—1 
>} + (s= 1,2,...,2— 1), 
(2.1) 


= /(3), >) — ré- 


s=1 


With respect to the nonlinear characteristic f(o) of the seryomotor it is assumed that it satisfies all the 
conditions outlined in monographs (5, 6], and in particular that it can be represented in the form 


(2) = ho + 9(9), (2.2 


where h is a positive constant distinct from zero, and function (qo) is such that for all the values of of 0 
it satisfies the condition: 


op(s) > 0. (2.3) 
System (2,1) on the basis of (2,2) can now be written: 
te = >) + hep (2) (s=1,2,...,), 
k= 
(2.4) 


a= 2 — TP (2), 
=1 


* Here, however, the law of subscript symmetry will no longer hold, i,e., Oy. F Chg. 


667 


| 

n n 

1 


Nn = €, Pen = Nz, Pur = hjx, 


0, sen, 
h, = * n, qk = 2 (2.5) 
In view of (2.3) we shall plot Lyapunov's function for (2.4) in the form: 
V=@+ (2,6) 
0 


where Q is the quadratic form derived by the method given in the previous section, The derivative V of 
function (2,6) with respect to time obtained according to Eqs, (2,4) is represented by the expression: 


—V =F 4+ rg*(s) + 29 (2) >) gun, 


a= — + (xn + Onk + Gx), (2.7) 


and F is the quadratic form whose coefficients satisfy Sylverster's inequality (1.5). Hence, by adopting a 
reasoning similar to the one in [7] we obtain the inequality: 


< A, (2.8) 


‘which in conjunction with (1.8) constitutes the stability criterion of controlled systems (2,4), Here A stands for 


the discriminant determinant of form (1.4) and A,, for a minor (including the sign) of the corresponding element 
of the determinant, 


Thus, the following theorem holds: 


In order to guarantee an asymptotic and absolute [5] stability of undisturbed movement of controlled 
systems (2,4) *at large” it is sufficient to choose their parameters according to inequalities (1.8) and (2.8). 


Inequality (2.8) depends on the coefficients of the form (1.4) whose choice is not limited by any con- 
sideration except Sylvester's inequalities, This circumstance narrows down considerably the practical applicability 
of the criterion, In this connection various modifications of the criterion obtained by means of the choice of 
coefficients of form (1.4) are of practical interest, Below, we shall give some of these modifications. 


2. Let us take the coefficients of Hurwitz's matrix (1,13) as coefficients for expression (1.4), Then the 
inequality will assume the form: 


n n 
Dd < rH. (2.9) 


k=1 


Let us illustrate this case on a system of indirect control [8]: 
Tap + p =0, 


= ag, (2,10) 
b= 7-2 +00), 


| where 

-§ 

| 

| 
| where 

| 
nn 
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9) 


10) 


Assuming ~ = 14, B = t= 7Tag, where r is the relative time with respect to which the derivatives of 
Ny 1, and o are taken, By eliminating variable n Eq. (2,16) is given the forms 


= Put + Piste 
Ne = Put + + hep (9), 


(2,11) 
+ 
= + 92% — 
where 
r, 
Pu = —9, Pu = —1, Po = —T 
h=r=T, fji= =, js=—1, 
+. 
Conditions (1.8) in this case are reduced to the inequalities: 
64+7>0, (2.12) 
and condition (2.9) assumes the form: 
— 8182 + (2,13) 


Since 


+ +9) and ++ a); 


the values of O43, Og;, and Og3 should be determined, From (1.18) we have 


= —i, =.— 7? (04-4), 


=—T(0+-<), = T +07 (04 4) 


Thus, quantities gy and gg which form part of the condition (2,13) are determined from the formulas 


= — +) (207 +4 0-4) +4], 
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and hence from ( 


a 
H=T(0+47)(0+ 4). 
3, Let the choice of the coefficients of (1.4) be subject to conditions (1.15) where a, is determined ‘rom 
relationship 
+ Gx = 0 (A == 1,2,...,”). 


After transformation by means of (1,11) and (1,12) these relations can be written in the form of the 
following system of solving equations with respect to the unknown a,: 


n 
>} = Yo = 1,2,..., 2), (2,14) 
k=l 
where 
% = — Hq, (2.15) 


and c,), are certain constants which depend in general on the parameters of the controlled systems. If we now 
assume that as the result of the solution of Eqs, (2,14) expressions for a;, a2, ..., 4, have been obtained in 


terms of the regulator parameters j,, Int 
| Oy = fay (2,16) 
the stability criterion of the controlled systems is reduced to fulfilling the inequality 
(9 (2.17) 


which deterinines, in conjunction with (1.8),the region in which regulator parameters should be chosen, 
Let us illustrate this on the previous example (2.10) for which the system (2,14) will take the form: 


where 
Cy = —TT OTT, t= T (a + 08) (0+ TY, 
Where we find: 


—7"]. 


Hence, the stability criterion of an indirectly controlled system amounts to the fulfillment of the in- 
equalities: 
(a — 87) > (0 + T) (a + 08) + 86%), 


a (8 + 27) > 87". 


(2,18) 


| 

| where 

| 

| 

| 

| 

| 670 


18) 


Let us note that in the case of a positive self-adjustment of the machine (6 > 0) the fulfillment of just 
these conditions is sufficient for stability of a system (2,10). 


Above results are greatly simplified if the last term is omitted from Lyapunov’s function (2.6) and it is 
taken instead in the forms 


V = Dd 


k=] 


Then, its derivative V with respect to time obtained by means of Eqs. (2.4) will be expressed bys 


—V=F+ op 
if coefficient ag of form F be chosen on the basis of conditions: 
2onn + jx = 0 (k = 4, 
These equalities can again be represented in the form of Eqs, (2,14) where y, is then determined from the 


formulas: 


— Hj, (s=4,2,...,n), (2,19) 


In the above example the stability criterion of system (2,10) will bes 


(2,20) 


4, Let us deal with one more simplified criterion of stability, It is of a more general character than the 


criterion (2.8), it is, however, of interest in many particular cases in solving problems of automatic control 
theory. 


Eliminating variable € from system (2.1) we obtain: 


n—1 
= = 1, 2,.... — 4), 
kel 


n—1 (2,21) 
& = >) qum —(p + rh) o —re(9), 


k=l 


where 
n, 
= + (mast, r#0), 


n—1 n—1 


s= 


Let Lyapunov's function for this system have the form: 
n—1 


V= > + 


k=] 


By repeating the former reasoning it is easily seen that the derivative of this function with respect to time 
obtained by means of Eqs, (2.21) will be represented by an explicitly negative function 
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— >} — (p rh) — 


if these inequalities are fulfilled: 
The conditions given below must be added to (2.22) 
Ree, <0 (s=1, 2,...,m—1), (2,23) 
where €, is the root of the equation 
(— 1)*"*A = = | = 0. (2,24) 


As a first example, let us examine a combination of automatic control system equations (6), p. 107): 


= + — Mee, (2,25) 
= Gim + — (p + rh) —re(o), 


1G?) 
+1 ) q2 a + 


— pG*) 
Vi(T?+iG) 


The condition of the real parts of Eq. (2.24) roots being negative is reduced in this case to one inequality 


—m,=r=i. 


(2,26) 
Further, we find 
a; = 91921» as = Gi + 92922- (2,27) 


Then, system’s (2,25) stability criterion is reduced to the fulfillment of the inequality 


T* + 1G? )+ |>° 


— pG?) 


(2,28) 


Thus, the controlled system (2,25) is asymptotically and absolutely stable “at large” if its parameters are 
chosen according to the inequality (2,26) and (2,28), Let us note that these inequalities were previously deduced 
by A. M, Letov [6] on the basis of Malkin's [7] modified criterion. 
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| where 
| 
| 
| 
| 
| 
m | 


6) 


B) 


As another example, let us examine the problem of a governor with an isodrome* attachment (6), p. 69): - 


+ 


Hence, in A, M, Letov's notation we obtain a system of equations 


= 
Ne = + — M29, 


Ns = 
= + + — (Pp + rh) — re(s), 


where 
qe. = (a — gG*), 


ms = p Vis 


+ 
— — por) 1 (a — (B — por) 


1 


9 (9) = (7? + 1G) V (0) — he. 
Conditions (2,23) in this case amount to the fulfillment of the inequalities: 
k+ al k+al i 
For the coefficients of the solving equations 


+ + = Yi, + Cope + = Yo, 
+ -|- = Ys 


we have 


C1 = —fifs, — fide, 1 = 


Cx = fs, ols, Cs = 12 = 
= 0, Ca = 0, = — 


* Publisher*s note: Isodrome =hydraulic lagging device, 


(2,29) 


(2,30) 


(2,31) 
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e 

= 


al 4 
From (2.31), we obtain H = fa(fi/a— fs): 
= + — a, = 91 — 


Hence, the stability criterion of system (2.29) amounts to the fulfillment of the inequality 


E> 
Vi (T? + U +18 (U + 1E)(E — pG*) 
4 pG*) 
(2,32) 
Vi (Tt + 1G) — _ + E— (U + 16)(E — pG) 
T? + 1G? (+ 
> al (k + al) (E — pG?) 
Vit? + 1G) — 
Another inequality should be added to (2,32) 
(2,38) 


Vi (T+ IG) 
which follows from the second condition of (2.22). 


Thus, the controlled system (2,29) is asymptotically and absolutely stable “at large® if its parameters 
satisfy the inequalities (2.30), (2.32) and (2.33), Let us note that these inequalities in the absence of the isodrome 
action of the governor (N = oo) transform into inequalities obtained for system (2.25), 


As the last example let us examine the controlled system represented by the equations 


$ + 1,23h — 0.6168 +. 8.8¢ = 0, 
6 + 0.1338 —¢ =0, 
E=f(c), o =0.136 a8 + bb + 0.73 ah — 1.93 ak, 


where a and b are the regulator's design parameters chosen on the basis of the system's stability. 
Assuming B = ny, ¥ =Ne v =n let us write Eq, (2.34) in the form of (2,21): 


= 

= + + — (Pp + rh) — 
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. | It is assumed here that 


qu = — 0.133, dis = 1, Gn = —0.004, 


m=—0.0210, q=—3.33b, 
p= — 3.33, r=1,93 a. 


The condition for the real parts of the roots to be negative in Eq, (2,24) in this case consists in fulfilling 
the inequalitys 


2.86 a + 20.786 > 0. (2,35) 
The coefficients of the solving Eqs. (2,31) have the form: 
Cy = 0, = — Cog = 0, 


Here, it is assumed that 


4.69, fo=0.597 44.562, 0.5932. 


Since, in this case 


i 
= — — 99911981 — (72 + 93), 
f . 1 
a, = — a; = + do + 939s), 


it is easily seen that the stability criterion of system (2.34) amounts according to (2,22) to the fulfilling of the 
inequality 


e+-rh>O, 91+ 492 + > 


In the terms of the initial conditions, we have: 


14.53. a? — 3,33ab + 20.79>0, h> (2,36) 
Since for a > 0 and b > 0 the inequality (2,35) and the first condition (2,36) always hold for the assumed 


numerical values, the stability criterion of system (2,34) amounts to the choice of the servomotor time T, from 
the condition 


0<T,< 


3. Choice of Parameters for Stable Controlled Systems 


The simplified criteria for the stability of controlled systems, derived in the previous section, provides the 
possibility of determining the permissible range of the regulator parameter variations, 


If the number of parameters is greater than two, the plotting of a stability region sometimes becomes 
difficult. In such cases it is advisable to proceed as follows. 


where 
3) 
ome 
615 


The values of a, = a%, found from the resolving Eqs, (2.14), are substituted in the expressions for the 
quadratic forms F and Q, Let these expressions be, respectively 


n—1 
F* = >} (3,1) 
n—1 
C= > (3.2) 


Let us now bring into consideration integral: 


If this integral is used to evaluate the quality of controlling the linear part of the system (2,21) represented 
by equation 


1, = > Vox (s=1, 2, ..., a—14), (3.4) 
k=1 


this integral can be found without a preliminary knowledge of the roots of Eq. (2,24), In fact, since we have 
Eq. (2,23) and the equality Q* = —F* we obtain from (3,3) 


> te, (3.5) 


where nj is the initial value of the solution n;(t) of the system (3.4) fort =0, Integral (3,3) is thus expressed by 
means of the regulator parameters jy, ++ 4? 


(3.6) 


The quality of system (3.3) control from the point of view adopted by us will be best when regulator param- 
eters jy, jn... +» Jp-3 are chosen so that integral (3.3) is minimized by n-3 of these parameters [9]. This leads 
to the solution of the equations: 


al 
a, =% 


from which it is possible to find all the j, with the exception of two values j, and j,, for which it is necessary 
to plot a stability region of system (2.21): jg=ig Gy, jy # v, 


If these values of j, are substituted in conditions (2,22) and (2,23), it becomes possible to find after certain 
transformations the required region of stability: 


ju) > 9. (3,7) 


Let us illustrate the above on the system of equations (2,25) for which we haves 


. 1 
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3.7) 


Let us assume for simplicity that nj =1, Then from (3,5), we obtain: 


From the condition of the minimum =0 and = =0, we find, respectively 


i 
E= (pl@—U FIO). 
The stability criteria of (2.26) and (2,28) will now acquire the form: 


2(qy + pV zy) > p*y 4+- 32, V zy > py, 


A4pqy? +- 7x V zy +- p’y V zy > + p*y’, 
py+2ah>a2 (r=l; y=T? + 


The first three inequalities will after transformation acquire the form: 


72? — 3 (p* — 29) 2y — 2pty* > 0 


kot < y < kya, 


where 


k __ 3 (2q — p*) + V9 (2q — p*)* + k __ 3(2q — p*) — V9 — p*) + 56 p* 


Now it will be easily seen that the stability region of system (2,25) determined by these conditions lies 
in the first quadrant of plane x, y, and is bounded by straight line y = kyx and the OX axis. 
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STABILITY OF NONLINEAR CONTROLLED SYSTEMS 


E. N. Rozenvasser 


(Leningrad) 


Certain circumstances which provide the possibility of extending the 
sphere of application of the method proposed by A. 1, Lur’e for investigating 
the stability of controlled systems are outlined. 


1. In this article all A. 1. Lur’e's notations and premises are retained unless otherwise specified, In [1] 
two types of differential equations systems representing the movement of automatic control systems are examined. 


The combination of direct control system equations 


te = beatae + hef (2) (k= 41, 2, ..., 
fete, of (2) >0 
s=1 


is reduced by a normal linear transformation to the form: 
Zp = det, + f (0) (p= 4, 2,..., a), (2 


c= D 1%» (3) 
o=1 


where d, (p =1, 2,..., m) are the roots of the equation: 
D(A) = | bea — Brak | = 0 (4) 


(6, is Kronecker's delta). 
It is assumed below that all quantities Xp are distinct from zero and that their real parts are negative, 
The quantities y, comprising (3) are determined by the relationship 


where polynomial Hy, (,) is adeterminant | | im which the k-th column is replaced by the 
quantity hy. 
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Taking an explicitly negative function 


@=1 B=) 


let us calculate the derivative a by means of Eq. (1)s 


By adding to (7) the quantity 


and equating to zero the coefficient of f(o) we shall express the stability conditions found by A. I. Lure in 


terms of the solubility of the equations: 
me a 
— 2 « (ep =1, 2, ..., a). 
Similarly the combination of the indirect control system equations 
= + mk (k=4, 2,...., 
a=1 


= het, + f (2) =1,2,...,m), 
= — rf (2), 
o=1 


N, 05) 
h 


(6) 


(1) 


(8) 


(10) 


(12) 


and the polynomials N,(\) are set up by means of quantities ny similarly to the Hy(\) polynomials. By means of 


the explicitly negative function 


n n 2,8 a 


it is possible to obtain the stability conditions expressed in the solubility of equations: 


2a, Vr — 2a, +B, =0 (p= 4, 2,..., 


a=1 


(13) 


(14) 


n n 
|_| 
(9) 
a n 
| 
| | 
| (11) 
| where 
: 
| 
| | 
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It is moreover important to note that Eqs. (14) are obtained with the assumption that 


r>0. (15) 


In the case of a direct control system a similar method of obtaining stability conditions can be adopted. 
In fact by differentiating (3) and in conjunctions with (2) we finds 


n 
= — rf (2), (16) 


where 


Combining relationship (16) with Eq. (2) we obtain a system which in a formal manner coincides with the 
canonical form of indirect control system equations, 


Thus, it is possible to obtain stability criteria for directly controlled systems by two methods: by means 
of Eqs, (8) and by means of Eqs, (14) where quantitiesr and Bp are given by the relationship (17). 


It is interesting to compare the two methods and determine the spheres of their application, Below,we 
give certain considerations in this connection, 


2. Let us examine a slightly more general system of equations than the one investigated by A, I, Lur*e: 


n 
Dd) + hyp (3, t) (k = i, 2, n), (18) 


>; (19) 


where bj,» hy, and js are the same as in system (1) and g(o, 1) is a nonlinear function of o and time t such 
that (0, ) =0, and that there is an instant t* which provides that 


t)o>0 for t>?’. (20) 
Moreover, let us assume that the system under consideration has a unique condition of balance 
¢=0, Zp =0 (p =1,2,....8). (21) 


ALL, Lur‘e's transformation is applicable to system (18) which can be reduced tos 


Zo +9 (2, t) =1,2,...,™), 


(22) 
a= 
where and yp are the same as in (2) and (3), 


Further, repeating exactly the deduction of the stability conditions for system (2) and (3) by means of 
Lyapunoy's function (6) we find that the solubility conditions for Eqs. (8) will ensure the stability of the balance 
conditions of system (18). In particular, stability criteria are obtained in this way when 


(2, t) = $(t)o, (23) 
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p(t)>0 for (24) 


if the system under consideration is linear, i.e., if function ¢ (o, t) is determined according to (23), it 
is possible to obtain a more general condition for stability assuming that 


e(t)>—M for t>?’, 


(25) 
M> 0. 
For this purpose, it is sufficient to examine similarly the systems 
= Dy + hyp(t, o) (k=1,2,...,2), 
a=] 
(28 
where 

= bra — Mhja- (27) 


3. It should be noted, however, that the stability criteria obtained for stationary systems of the form of 
(1) by means of Lyapunov's function (6) are often difficult to construct. This is due to one of the necessary 
conditions of solubility of system (8) being: 


n 
r= — > 0. (28) 


Condition (28) often leads to an unnecessarily complicated circuit structure of the required automatic 
control system. 
The application of Lyapunov's function (13) provides the possibility of obtaining stability criteria free of 


this defect. Moreover, the solution of the stability problem of an n-th order system is reduced to solving n-1 
quadratic equations. In fact for system (1) we have on the basis of (18): 


(29) 
e=1 


but one. of the solving equations after transformation proposed by A. 1. Lur‘e becomes [1]: 


2 V — 0, (30) 
n 

where ty. is unknown and = >) is a real constant 

From (29) we find 
682 


TE 


if 

dk 
| 

| 

| 

} 


9) 


0) 


31) 


%* 


Since r = 1, Eq, (30) is always soluble, By substituting this value of tn-4 into the remaining equations, 
we obtain a system of n - 1 equations. 


It should be noted that in the literature on this question, it was erroneously considered up to the present 
that the necessary condition for the solubility of an equation of the form of (14) was in the general case: 


Sees rso, (32) 
o=1 


on this basis an incorrect conclusion was drawn that if the initial system of (10) and (11), before they were re- 
duced to the canonical form, was (1), with roots of the characteristic equation of an open loop system distinct 
from zero, the combination of the quadratic Eqs, (14) is insoluble, It can easily be shown that inequality (32) 

is a necessary condition only in the case of an indirectly controlled system (9); it is obtained, however, from the 
conditions of Raut-Hurwitz for a linearized system [1] instead of the solvability conditions of type (14) equations. 
For a general case of a directly controlled system condition (29) is fulfilled identically and has no direct relation 
to the Raut-Hurwitz conditions for a linearized system, this is also confirmed by the example cited below, 


4, As the simplest example of applying above reasoning let us examine by means of two methods a direct 
control system of the second order represented by the equations: , 


+ Ay f (3), 
= Agi + + (9), (33) 
= + foo. 


The Raut-Hurwitz conditions for a characteristic equation of an open ioop system (33) are: 
My + Gag <0, — > O. (34) 
From (5) we calculate polynomial M(\): 
M (hk) = —ri—P, (35) 
r= — — 
P = — hyQye) + — hy (36) 
The solvability conditions of Eqs. (8) corresponding to system (33) are: 


M M (,) 


By using the summation formulas [1] and in conjunction with (34) it is possible to write conditions (37) in 
the form: 


P>O0, r>-~ (38) 


In the plane of parameters r, P condition (38) represents a portion of the first octant between the axis 


r> 0 and the straight line r= — —t—. 
441 + 423 
By the virtue of the above, conditions (34) and (38) will ensure stability in the case when in system (33) 
f(a) is replaced by function (20). By using the reasoning of section 2 it is possible to show, for instance, that for 
a general asymptotic stability of a trivial solution of equation 


= 
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+ la + bp las =O >0, >0) 


with an arbitrary function ¢(v) it is sufficient for the coefficients of % and x to be larger with t > t* than any 
arbitrarily small positive number, 


Let us now examine system (33) by the second method. In this case in view of the reasoning of section 3 
the problem is reduced to the solution of one quadratic equation. After transformations we obtain the stability 
conditions in the form: 


r>0, rg— (39) 


where r and P are the same as above. 
With r=0, we obtain the single condition: 


(40) 


It also follows from this that a trivial solution of equation %+ ax + f (x) =0 is asymptotically stable with 
any function f (x) represented in the form: 


o(z)z>0, 9(0)=9, 
providing a > 0. 
This result cannot be obtained by the first method, 
In the plane r , P,condition (39) represents a part of the first octant below the straight line 


P (a3, + 


431429 — 21949) 


Thus, the two methods give in general different regions of stability which for 


(441429 — > (4, + 
do not cross, 
If, however, 


01109 — < (44, + (41) 
regions (40) and (39) will cross and cover the whole of the first octant of plane r, P: 
P>0. (42) 
Conditions (42) represent the necessary and sufficient conditions of stability for system (33) with 
{(e)=co, 00, 
Hence, in the case of (41) the application of both methods provides the widest condition of stability. 
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FILTERING OF A CERTAIN TYPE OF NONSTATIONARY RANDOM VARIABLE 


I, A. Boguslavskii 


This article deals with a method of designing a circuit for generalized 
unshifted—— signal filtering of a random function, whose useful signal is re- 
presented by the solution of a first-order linear differential equation. The 
coefficient of this equation consists of a function of time unknown in advance, 
but continuously determined in the process of filtering, and the right-hand 
side of the equaiion consists of the sum of the products of known functions 
and unknown quantities added to a certain function of time, which is con- 
tinuously determined in the process of filtering. 


Posing the Problem 


In the practical solution of modern automatic control problems there often arises the necessity of using 
filters for picking out a useful signal from a background of noise, Moreover, the unshifted—— signal filtering, 
or as it is often called,filtering without dynamic defects, is of considerable practical interest. In this case the 
signal at the output of the filter after the ending of the transient process must be the same as ai its input, 


Papers [1-3] deal with the problem of unshifted filtering optimum with respect to the value of the mean 
square error, of a nonstationary random process in which the main part of the useful signal consists of a sum of 


the products of predetermined functions of time, and certain quantities which remain constant for each applica- 
tion of the process, 


In a number of cases, however, the main part of the useful signal has a rather more complicated form 
consisting of the solution of a structurally predetermined differential equation whose right-hand side includes 
the same sum of the products of the known function and unknown quantities added to a certain function of time 
unknown in advance but continuously determined in the process of filtering, 


The coefficients of the ‘differential equation are similarly functions of time,unknown in advance but con- 
tinuously determined in the process of filtering. Moreover, for the sake of easy application the filtering circuit 
should be designed in such a way that the solution of the differential equations with its right-hand side consisting 


of the above-mentioned continuously measured function of time is used, rather than that obtained, for instance 
at zero initial conditions. 


The latter condition is expedient because the possible errors in measuring this function will produce in the 
solution of the equation (and hence at the output of the circuit) errors which, with a certain structure of the 
differential equation, will increase with time, 


The problem of filtering of this type arises when designing automatic control systems in which angular 
coordinates (or their derivatives) are used for a moving object. 


Moreover the system of coordinates itself is on a moving base (ship or plane) whose parameters of move~ 
ment are unknown in advance, but can be constantly determined by special transducers and taken into considera- 
tion in the filtering process. 


a 
q 
q 
q 
| 
| 
) q 
2) 
| | 
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Often in these cases it is useful to solve a more generalized problem, which we shall conventionally call 
a generalized filtering problem. Moreover, the output signal should consist of the useful input signal added to 
the sum of the products of the above-mentioned unknown quantities and certain given functions of time. 


We shall examine the design of a circuit for generalized unshifted-signal filtering in the case when the 
differential equation is linear and of the first order. In particular, we shall show that if the above sum contains 
only the product, the circuit can be constructed fairly easily from standard components. 


As an example, the above technique is applied to the problem of smoothing out an angular coordinate of 
an airplane, which is transmitted by a radar station installed on a swaying foundation. 
Generalized Unshifted Filtering Method 
According to the posing of the problem the random function x (t) to be filtered is 


a(t) = y(t) + Ay(t), 
where Ay (t) is the random function of the noise signal and y(t) is the useful signal which satisfies the equation 


y(t) = a(t)y(t) + b(t) + Dd) cefe(t). (1) 


k=1 
In Eq. (1), f,(t) is the unknown function of time, a(t) and b(t) are the functions of time determined in 
the process of filtering and c; are the unknown quantities constant during each application of the process. 


It is required to design a circuit which with an input x(t) would provide an output function whose mean 
effective value after the completion of the transient process in time T corresponds to function 


2(t)=y + (2) 
=1 
where y(t) is the given function of time. 


By T we shall understand the time of the practical (in the engineering sense) attenuation of the transient 
process or the time of its accurate attenuation which presupposes the use of filtezs with a finite memory, con- 
structed for instance on the basis of the discrete counting technique, In the first instance some of the relations 
given below are approximate. 


Thus, the circuit must have the property of producing at its output after the lapse of time T function z(t) 
when function y(t) is fed to its input. Moreover, for reasons given above it is not permissible in designing the 
circuit to use the solution of the equation: 


yr(t) = a (t) y (t) + b(t). 


Let us note that if we assume that ~;(t) =0, we shall obtain the problem of ordinary (and not generalized) 
filtering with a useful signal represented by Eq. (1) with an unknown initial condition. 


The circuit will be designed on the basis of linear, time-sensitive filters. Let W(t, r) be a pulse transfer 
function of a certain filter. 
By immediate checking it is easy to prove this identity: 
t 
= E(t)y (0) + (WE, + Z(t, 
0 


t n t 
| de (t) + fe(*) dz], (3) 
k= 


E(t, =1 E(t)= E(t, 0). (4) 


Let W(t, r) be chosen so that the conditions below hold: 


wie (5) 
t 


From (3) we obtain for t < T 


t 
= E(t) y (0) + \IW(t, ) + E(t, + 
0 


t n T 
+ BU fx (2) de, 
and fort > T 
t t 
2(t)=~ (t, t) + E(t, 2) b (2) de. (1 
0 0 


The technique of designing a circuit for a generalized unshifted (by time T) filtering follows from (7). 


1. A filter is designed to have a pulse transfer function W(t, r) which satisfies condition (5) and to be such 
that function E (i, r) related to function W (t, r) by Eq. (4) satisfies conditions (6), 


2. Function x(t) is fed to the input of the filter with the transfer function W(t, r), and function a(t) x(t)+b(t) 
is fed to the input of the filter with the transfer function E(t, r), This function can always be constructed since 
according to the conditions of the problem functions a(t) and b(t) are con- 


wt) tinuously being determined in the process of filtering. 
Witr, 
= z(t}: 3. The outputs of both filters are added together. In conditions (5) 
aft)yft) *bit) and (6) the lower limit is given as t, (and not zero) in order to make the 
wt) 


circuit provide as from an arbitrary instant t, generalized unshifted-signal 
filtering when function x(t) is fed to its input. 


The circuit connections of a generalized unshifted-signal filtering 
are given in Fig. 1.° At the output of the circuit after the lapse of time T (from the moment when E(t) =0 
when condition (6) begins to be fulfilled) function z(t) with certain random errors will be received, 


Let us denote the dynamic error in building up function z(t) by Ag(t). For any value of t this error is de- 
termined by the relation 


Fig. 1 


t 
a(t) = E(t)y 0) + + EC, 


k=l 0 


* Notations in the circuits (Figs. 1, 2, 3, 4, and 5) are given without taking into account noise in the filtered 
random process. 


where 
t 
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For t > T, the variance o3(t) of random errors in constructing function z(t) can be found from formula 


tt 


where oy(4, ty) is the correlation function of the random function Ay (t). 


If it is known a priori that the initial value y (0) of function y(t) is large, the requirement that the quantity 
E(T)y(0) should in practice be close to zero can be fulfilled only by sharply decreasing the filtering qualities of 
the circuit. Hence, it is sometimes advisable at the cost of a rise in random errors at the beginning of the filtering 
process to eliminate the dynamic error E(t)y(0) for any value of t . For this purpose it is sufficient to store the 
initial value x(0) of the filtering process, multiply it by function E(t) and add the product to the output of the 
system described above, The dynamic error which is proportional to y(0) will disappear but the random error 
will increase by the amount E(t) Ay(0), which is small fort > T. 


It is easy to vary the technique of designing circuits of generalized unshifting-signal filtering, Let us 
transform random process x(t) to another random process x,(t) by means of a certain linear conversion S which 
has the property of converting y(t) into function y,(t) which satisfies an equation of the type: 


Yo (t) = (t) y(t) (t) + (10) 


k=1 


Then having applied to variable x,(t) the technique developed above and having subjected the result of 
filtering to a reciprocal linear transformation S~*, we obtain a generalized unshifting signal filtering of the initial 
random variable. 


Let us now design circuits both for simple and generalized filtering in various particular cases. In designing 
let us assume that we have at our disposal only normal components used in modern computer technique: quadri- 
- poles (filters) with constant and functional parameters, functional potentiometers, wound according to a given law, 
with a slider driven by a special low- power servyo- mechanism. 


Let us examine particular cases in which our problem can be solved by means of above components. 
1, Filtering without dynamic errors a random variable in which the useful signal y(t) satisfies the equation: 


y (t) = a(t) y(t) + b(t). (11) 


The problem is a particular case of the general problem when f(t) = y,(t) =0. Filters with constant para- 
meters should be used. 


Let K,,(p) and Kp(p) be transmission functions of filters with pulse transfer function W(t) and E(t), respectively, 
From condition (5) and from (4) it follows that 


Kw(0)=1, Ke(p)=— ag 


Thus, the filtering circuit without dynamic errors is designed in the following manner (Fig. 2), An arbitrary 
filter is taken which has a transfer function Ky(p) which only possesses the property that when its input is subjected 
to an instantaneous change in voltage its output in time T becomes practically equal to the input. Another filter 
is designed with a transfer function Kp(p) according to formula (12). The input of the first filter is fed with function 
x(t) and that of the second with functions a(t) x(t) + b(t). The outputs of the filters are added. 


It follows from the above that,in time T, function y(t) will appear at the output of the circuit without any 
dynamic errors,and with random errors,in general smaller than at the input. 


ering 


11) 


ectively, 


(12) 


itrary 
»jected 
filter 
function 


any 


The circuit with an unshifted-signal filtering becomes especially simple if we assume that 


Kw (0) = 557" 


Then, for obtaining a filtering with a signal not shifted by time T the input of the delay link with the time 
constant T, should be fed by function 


(1 + T; a (t)) 2 (t) + T, b(t), 


where w T. 
ult) 
Fig. 2 Fig. 3 
2, Filtering of random variables in which the useful signal y(t) satisfies the equation 
y(t) =a(t)y(t) +b) (13) 


here c is an unknown quantity which is constant at each application of the variable, f(t) is the given function 
to be differentiated. 


This problem corresponds to the general problem with n= 1 and p(t) =0. 


Let us assume that x,(t) = oe Then: 


ys (t) = a, (t) +, (t) (14) 
where 


a, (t) b, (t) = 


Let us design a circuit with an unshifted-signal filtering for a variable with a useful signal y,(t). From 
conditions (5) and (6) it follows that 


Kxg(p) = pK (p), Kw(p)=1— (p). 


In order to ensure the filtering properties of links with the transmission functions Kw(p) and Kg(p) it is 
necessary that for |p| -» oo the following conditions hold 


p> K(p)—>1, pK (p)+0. 
Hence, the transmission function K(p) must have the form: 


_ Mp) (15) 


where M(p) and N(p) are arbitrary polynomials with powers smaller that m—2 and m, respectively, where m 
is an arbitrary number not less than 2. 


| 
0) | | 
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Thus, the circuit for an unshifted-signal filtering in the case under consideration is designed in the follow- 
ing manner (Fig. 3): 


a) the filtered variable and function b(t) are divided by f(t) as the result of which variable x,(t) and 
function b,(t) are formed 


b) from the division results signal a,(t) x,(t) + bg(t) is formeds 


c) variable x,(t) is fed to the input of the filter with the transmission function 1 — p*K (p) and the newly 
formed signal to the input of the filter with the transmission function pK(p); 


d) the outputs of the filters are added and the result multiplied by f(t). 


The dynamic errors in constructing function y(t) will disappear from the instant when the function below 
becomes in practice equal to zero: 


t 


Above circuit can be used for unshifted-signal filtering of a variable with a useful component of the form: 


(t) = (t) + 
where go(t),y,(t) are given differentiable functions. 
In this case 


a(t) = 29, = () oun. 


3. Let us now filter in a generalized manner a random variable whose useful signal satisfies Eq. (13). 
According to the posing of the problem the circuit should be chosen so that after the lapse of time T and without 
noticeable dynamic errors, function z(t) = y(t) + cy(t) should be produced in which »(t) is given, 


By dividing the filtered variable x(t) by f(t) we shall arrive at the problem of a generalized filtering of 
variable x,(t) whose useful component y,(t) satisfies the Eq. (14). In this case z(t) becomes: 


(t) = yo (t) + chs(t), 
where ¢, (¢) = Tiny 
For a generalized unshifting-signal filtering of variable x(t) let us choose W (t, r) according to the equality 
W (t, t) = W, (t — t) + Wy (t —*), 
where W,(t) and W,(t) must be found. 
It is obvious that 


E(t, t) = Ey (t — +t) — $.(t) £2 (¢ —*), 


t t 


. and hence 
where 
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Let us denote by Kw, (p), Kw, (p), Kg, (p) and Kg, (p) links with pulse transfer functions W,(t), Wg(t), E(t), 
and E,{t). 


It is easy to show that the required conditions will be fulfilled if 


Kw2(p) = pK2(p), Kxa(p) = K3(p), 


where Ky (p) and K,(p) are transfer functions which have the initial conditions of Ky(0) < @ and K,(0) =1. 


In order to ensure the filtering properties of the circuit the transmission function K, (p) should have the form 


of (15) and the power of the denominator polynomial of function Ky (p) should be at least two powers below that 
of the numerator. 


After the links with the pulse transfer functions W(t, r) and E(t, r) have been selected, the circuit for 
generalized filtration without a signal shift is designed according to the technique described earlier (Fig, 1). 


It will be seen that the circuit is made up of links with constant parameters and functional potentiometers | 


a with a programed displacement of the slides. 


Fig. 4 

Dynamic error in conbessiteig function z(t) will be absent in this circuit from the instant functions Aq, () | 
and Ag2(t) which are determined by the equality given below and in practice equal to zero; 

ality t t 
da (t) {70 Wax] + ear}, | 
0 0 
t t 
0 0 


It has already been pointed out that by storing quantity x(0) and constructing a corresponding additional 
circuit it is possible, at the cost of increasing the random errors, at the initial stage of filtering, to eliminate 
the dynamic error in Ags (t) for any value of t . 


Kwi(p) = 1—pKi(p), = PKs), 
| 
| 
| 
out 
| 
| 


As an example, Fig. 4 shows a concrete circuit for generalized unshifted-signal filtering with: 
1 i 


From the above circuits of filters for ordinary and generalized unshifted-signal filtering it follows that all 
these circuits possess filtering properties inferior to those of the integrating link, since all the enumerated transfer 
functions can be approximated at sufficiently high frequencies by a function of the form k/p, Hence, above 
circuits are only effective for filtering random errors, which 
can be considered as white noise or as the result of the 
transition of white noise through circuits which do not con- 
tain differentiating links, With random errors of the type 
of pure or slightly dynamically distorted derivatives of white 
noise, these filters are ineffective and even sometimes do 
not ensure the production ef the final variance at the output. 


In these cases a more complicated filtering circuit 
Fig. 5 can be used, which consists of series connected simple 
circuits. In fact, let a circuit described by the pulse trans- 
fer functions W, (t, r) and E;(t, r) provide an unshifted-signal filtering y (t)+ Ay(t). Following the transient state 
at the output of this circuit,signal y(t) is formed with decreased random errors. Let us feed this function to the 
input of a circuit with generalized filtering,which is represented by the pulse transfer functions W,(t, r) and 
Edit, r) Fig. 5). 


At the output of the network consisting of the two series connected circuits we shall obtain function z(t) 
with certain random errors which are determined by the filtering properties of the two circuits. If in the place 
of the first circuit the circuit described in section 2 is used and in place of the second that described in section 
3, the entire network at sufficiently high frequencies can be represented approximately by a circuit consisting of 
two series connected integrating links, which will also ensure effective filtering for random errors of the type of 
derivatives of white noise. 


Let us introduce the notation: 
W,(t, t) + = Vit, 
W, (t, t) +a (t) E; (t, t) = (t, t). 


Then 


t 


Den (t) = + + Volts) Bs (6s fe du 
0 


It can be easily shown that both errors tend to zero if the functions W,(t, r) and W(t, r) satisfy condition 
(5) and if for tt, > T the following expressions hold: 


t t 
Ea(t, de — (0), (t, fa(s) dt 0. 
0 


lition 


Optimum Generalized Unshifted Filtering 


Let us examine generalized unshifted (by time T) filtering which would ensure a minimum mean square 
deviation of the errors of constructing function z(t), Filters used for solving optimization problems and possessing 
a finite memory provide the possibility of accurately satisfying conditions (5) and (6). 


The design of an optimum circuit is possible, it would appear, without using any a priori data on function 
a(t). 


Let there be given a most probable application a 9(t) of function a(t) and function V(t, r) chosen which 
satisfies conditions: 


V(t,s)=0 when +t>T, 


t 


(VU, t)exp[— a9 (¢ —u) du] de =4, 


(16) 
(t, t)exp —u) du] de = — h(t), (7) 
where 
F, (t,t) = (¢ — t)exp ay (t —u) du| dx. 
0 
Let us represent function W(t, r) and E(t, r) by expressions 
t—t t 

E(t, +) = exp ay (t — t) dt {1— \ V (2, *) exp[— u) du] (18) 
W (t,t) = (19) 


It can be easily shown that functions W (t, r) and E(t, r) will precisely satisfy conditions (5) and (6) and 
hence circuit 1, made up of links with pulse transfer functions determined from formulas (17) and (18), will 
provide generalized unshifted (by time T) filtering with any form of function a(t), Let us find V(t, r) so that 
the optimization problem be solved for the case when a(t)=ag(t), Since V(t, r)=W(t, t- r)+a,(t-r)E(t,t- 1), 
it follows from (9) that the determination of V(t, r) which reduces to a minimum a (t) should be attained by 


means of the well-known technique, according to which V (¢, t) = do (t)Vo(t, t) + ye (t) Vx (t, *)- 


Here Vg(t, r) and Vj, (t, r) satisfy the integral equations 


exp[— —u) du| = u) ®, («, u) du, 
Fy (t, (t —u) du] (t, u) D, («, u) du, 


and functions ) 9(t) and d,(t) are determined from the condition that V(t, r) satisfies equalities (16) and (17). 
Then 


o8(t) = do(t) — (0. (20 
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The circuit of Fig, 1 made up of links whose W(t, r) and E(t, r) are related to function V(t, r) by 
formulas (18) and (19) provides : unshifted generalized filtering with any continuously determined function a (t) 
and provides a minimum for 03(t) when a(t) =a 


Shifted-Signal Filtering 


The previously examined circuits for generalized and simple filtering were of the unshifted type, i.e., 
ensured that after the lapse of the transient state time interval T, there was either an approximate or precise 
(with elements having a finite memory) elimination of dynamic errors. In a number of cases, however, when 
for instance the probability characteristics of quantities y(0) and cy, are given and it is known that the possible 
values of these quantities are few, it is advisable to allow for some dynamic errors Ag(t) and to reduce by this 
means the mean square deviation of the random errors o,(t). 


Let it be advisable for the solution of this concrete problem that the value of a certain function of Ag(t) 
and o,(t) expressed by F[ Ag(t), a(t] be as small (large) as possible, The value of this function depends on 
the random quantities y (0) and c, which are included in the expression for A g(t) and hence this value will also 
be a random quantity. Therefore, the filtering circuit can be considered good if its parameters ensure for any 
value of t that the mean effective value of random quantity F[Ag(t), o,(t)) reaches a minimum (maximum), 
Sometimes, in the absence of correlation between y (0) and cy the filtering circuit can be selected by striving 
for an extremal of function F[Ag(t), o z(t)] with a Ag(t) corresponding to the mean square deviations of the 
random quantities y(0) and cy. By using above technique, it is possible to solve this problem in the following 
manner. Let us substitute conditions (5) and (6) by conditions 


t 
1) fx (x) dt ~ — dx (t) — ee, (21) 


where quantities €9, €y,..., €,, are mot yet determined, and let us choose corresponding W(t, r) and E(t, r). 
Instead of (7) the following expression will hold: 


t 
2(t) = (t, + E (t,x) y (x) de + 


+) 2) de + yO) 


Hence, the dynamic error of constructing function z(t) my means of the generalized filtering circuit (Fig. 1) 
for t > T will become: 


n 
Ag=y(0)& + 
k= 


The function o,(t) which is determined by formula (9) will depend on the quantities €», €;,...,» €n 


By using the required data of the probability characteristic of the random quantities y(0), cy, . . . , Cp, it 
is possible to find F(t, €9, €4,.+..» €p) which is the mean effective value of function F[Ag(t), oz(t)). 


By finding in the usual manner the extremum of function F, (t, €9, €4,...» €p) we shall thus determine 
the optimum quantities €», €;,..., €p and will obtain, from the adopted point of view, optimum parameters 
for this filtering system. 


Similarly, it is possible to solve the problem of synthesizing an optimum filtering circuit for a (t) = a g(t) 
whose structure is not given. For this purpose conditions (16) and (17) should be replaced by conditions: 


«)exp[— ay (t -- u) = 
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t) y(t, t) exp[— ao (1 u) du] de == hy (t) — ey. 


However, the designing of shifted filtering circuits is complicated by the necessity of using integrating 
links. 
Let us illustrate this point on a simple filtering circuit with f,(t) = p,(t) =0. 


It follows from condition (21) that Kw(0)=1 — 9 # 1, and that the transmission function Kg (p) deter- 
mined by formula (12) is tending to oo with |p| +0, Hence, the link which corresponds to Kg (p) must contain 
integrating elements, Thus, if 


Kw()= 


then 


Ks) = 
The development of the technique of filtering which minimizes a given variance and dynamic error function, 
is also possible when using above circuits of simple and generalized filtering without displacement of t+ @ by 
selecting appropriate values for filter time constants, 


Unshifted Filtering at the Output of a Radar Set 


As a simple example, let us examine the application of this technique to the problem of obtaining an 
angular coordinate of an airplane, tracked by radar, which is located on a swaying base, We shall assume that 
the linear acceleration of the base (for instance of a ship) is negligible. 


The application of a certain statistical dynamics technique for determining an optimum transmission 
function of a radar tracking system or of some of its parameters is described in detail in[4]. There the angular 
coordinate was represented by a polynomial in powers of t and the solution of the problem of optimization was 
carried out with given error coefficients, As the result of it, definite requirements for frequency characteristics 
of the tracking system were found. 


Another posing of the question is also possible which uses to a fuller extent thea priori information about the 
law of variation of the angular coordinate, Let us adopt a requirement that the radar tracking system should only 
have a sufficiently wide frequency band which would ensure in practice the elimination of dynamic errors in 
tracking and had only a small rise in its amplitude-frequency characteristic, Then, we shall obtain at the output 
of the radar set a variable which will differ from the angular coordinate only by some random tracking errors 
which the tracking system will not increase at any of the frequencies as compared with an ideal nonreactive 
system, This variable is connected to the input of a system with unshifted filtering, Let us take advantage of 
the phenomenon that over a fairly wide range of conditions the angular coordinate g(t) with a iinear movement 
of the airplane is sufficiently accurately represented by the differential equation: 


9(t) +254 =e+b(0), 


where c is a quantity constant for a given course of the airplane, b(t) = p(t) + Fo got), Dit), Die) 


is the current distance to the airplane and its rate of change, which are continuously measured by the appropriate 


units of the radar set, gp (t) and g(t) are the angle and the speed of rotation of the base measured by the trans- 
ducers installed on the base. 
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In the case under consideration a(t) = — 2 ae f(t) =4. 


In order to be able to apply the filtering, function b(t) must be continuously constructed by a special 
computer. When, for instance, the filtering circuit shown in Fig. 3 is used, we shall obtain at the output of the 
filter, after the transient state, the angular coordinate p (t) 
approximately without d — errors and with reduced random 


Let us assume that a linear law of change in the distance 
2 to the airplane is most probable and that the random process of 
tracking errors is of the type of a white noise. Then 


= 


Fig. 6 
Tr Day (ti, t2) = 8 (ty — ty). 


It follows from (9) that when Fig. 3 circuit is used for filtering with 


KO) = 


the following relation holds: 


() +75 — t)(T, + 27; —t) x 
f 


In solving the problem of optimum filtering we obtain from (20) 


i 


It can be shown that in the circuit of Fig. 3 dynamic errors are practically eliminated fort > 4Ts, In order 
to compare under the same conditions the filtering properties of the circuit in Fig, 3 and an optimum circuit it 
should be assumed that T=4T,, Figure 6 shows the results of calculations of the ratio oy, (0/oy(t) with Te= 
= 5 sec and T;, = 100 sec. 


— 

where 
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It will be seen that with t > 25 sec the filtering properties of both circuits practically coincide, Hence, 
under the assumed conditions the circuit of Fig, 3 with a chosen form for K (p) is practically equivalent to an 
unshifted optimum filtering circuit. 
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SYNTHESIS OF SERVOSYSTEM CORRECTING DEVICES IN THE PRESENCE 
OF NOISE. 


P. S. Matveev 


(Moscow) 


The results provided in [3, 4] are summarized for the case when signals 
are applied to two elements of a servosystem and for systems of automatic 
stabilization, The suggested method is illustrated by an example. 


In a number of publications (see for instance [1-3]) theoretical and practical results have been obtained 
for determining the optimum characteristic of dynamic systems subjected to random effects. 


These results, however, are applicable to servosystems only when all the effects are received by one 
element of the servosystem, for instance the measuring element, Reference [4] deals with the problem of de- 
termining the optimum characteristic of dynamic systems when the effects are applied to two elements. Hence- 
forth, we shall conventionally call such systems dynamic systems where effects are applied at two points, The 
' defect in posing the problem in [4] consists in the impossibility of determining the basic requirements for the 
dynamic properties of the system to be designed, The book only deals with stationary random functions which 
are fed to the input of the system, and the accuracy of the system is characterized only by its mean square error. 
The object of the present article consists in extending the results obtained in[3, 4] to the case when at the input 
of the system there appear given functions of time in addition to the stationary random function and in applying 
the results to automatic stabilization systems. 


1, A Method of Synthesizing Servosystem Correcting Devices 


Let us assume that one element of the system under consideration (Fig. 1) is subjected to a control action 
y (t) and a disturbing action n(t), The control action is a sum of two components 


y(t) = g (4) + m(t), (1) 


where g(t) is a given function of time, and m(t) is a stationary random function with a given correlation function 

R,, (1) or spectral density S,,(w). The disturbing action n(t) is also assumed to be determined by a stationary 
random function of time with a given correlation 
function R,,(r) or a spectral density S,(w). 


u(t) 
? It is also assumed that the second element of the 
omit) Wels) system is subjected to a stationary random disturbing 
a + action u(t) with a known correlation function R, (r) 
or spectral density S,,(w), moreover, for simplicity we 


shall consider that the functions m(t), n(t) and u(t) 
Fig. 1. Block schematic of the servosystem. are not correlated and have zero mean values. 


The problem can be put in the following manner, From given correlation functions Ry(r), Ry(r), and 
R, (7) a given transient state time T and given error coefficients Cy, Cy,... , Cy find such a pulse transfer 


function k(t) for a closed loop system that the mean 
square error yalues are at a minimum, It is also assumed 
that the pulse transfer function b(t) [transmission function 


t 
Fig. 2 
where 


KWo G)] of the object is given. By means of known 
methods the transmission function of the 
correcting device (see (3]) is determined from k(t). On 
the basis of the aforesaid, it is convenient to use the 
circuit shown in Fig. 2 for the solution of the above 
problem. 


Let the ideal system produce at the output the 
signal: 


zit} 


h(t) = Hg(p)g (t) + H(p) m(t), (2) 


= 


H,(p) = H (p) — pm. (4) 


Assuming that the correlation function of variables m (t) + n(t) has 2q continuous derivatives and hence 
variables* g(t)+ m(t)+n(t) are differentiable g times, let us write down (see [5]) the expression 


r(t)= > (¢—2) + m9 4 — (5) 


j=00 


In expression (5), Wj(r) are functions which have limited variations and characterize the optimum system, 
It is also assumed that they are continuous in the interval (0, T) and may have discontinuities of the first order 
for r =0 and r =T, so that Stieltjes integral exists in (5). 


By integrating by parts with zero initial conditions expression (5) can be reduced toz 


where 


x(t) =| 


q q 
+ (t) + DBO (x —T)] dz, 
j=0 j=o 


2(t—t) = g(t —*) + m(t —t) 4+-n(t —1). 


In expression (6) it is assumed (see [6]) that: 


5 (2) 


q 
= th <<, 


ji=0 


k(t) =0 with Te<tqg ox, 


* It is assumed that g(t) is at least as “smooth” as m(t) and n(t). 
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T 
—2) 80 (e)de = 2651), 
(1) 


7 
(t —t) 8 T) de = (t — 7). 


In further calculations for the sake of simplicity we shall denote the square bracket of expression (6) by 
k(r) and only in the final solution introduce the delta functions at the ends of the interval, 


Since 
= O(s) ig (s) + m(s) + n(s)] + Ka Wa (s) [1 — O(s)} u(s), 
we obtain from (2) the error of reproduction (Fig. 2) in the form: 


e(t) = Hg(p) g(t) + (p) m(t) — \ [g(t —t) + m(t + n(t — 


u(t —e)b(s)do + u(t —t— 0) b(e)dc. (8) 


Assuming that 
= | 4) (9) 
g(t —*) ¥ ae? 


_and in conjunction with (4), we equate the mean value of the reproduction error to zero. As the result we obtain: 


T 
m= 0,4,2,....7), (11) 
0 


In considering Eq. (10) as an identity we can write (r+ 1) limiting values imposed on the pulse transfer 
function by the reproduction operator H (p) and the error coefficient Cj in the form 


T 
0 


T 


The expression for the reproduction error which conforms to conditions (12) and (13) or, which is the same 
thing, to (11) will assume the form (see [2]): 


100 


T 
e(t)= | m(t—+)x(x) —([m(t — 2) nt — — 
— | de. (14) 
0 


By similar reasoning as in [2] it can be easily shown that the necessary and sufficient conditions for a 


minimum in expression (14), when the limiting conditions (11) are observed, consist in the pulse transfer function 
k(t) satisfying the integral equations 


+ Ra(t — de k(x) de b(8)d® x 
xf Ry (t —t—9)b (0) d0 + tad 
Introducing the notation 
Rus(t) = fo (x) b (8) Ry (t — 0) a9, (16) 


let us re-write expression (15) in the following forms 


[Rm (t —t) + Ra(t—*) + Rus = 


= \ (t + Rus (t) + (17) 
i=0 


Let us examine the solution of Eq. (17) for the type of stationary random variables whose correlation 
function 


R(t —*) = Rn(t + Ra(t—*) + Rw (t (18) 


is related in a certain mannér to Green's function of a self-conjugate differential system. This type of system 
also includes, stationary random variables with a rational fractional spectral density of the form: 


bo+ +b, 
S = Sm(o) + Sn(@) + Suv (@) = 


(19) 


By using the relation of the correlation function to Green's function (see [5, 6, 7]) and in conjunction with 
expression (7) we can write the final solution of the integral equation (17) in the form of an optimum pulse 
transfer function 


ak r q q 
k(t) = >) Be’ + Adi + >) (t) 4- DS —7) + 
io 


+ L(p)L*(p) M-*(p) M**(p) [Rus (t) + | Rm (t 
(0<t<T), (20) 
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where gq is determined from expression (19) and is equal to 


q=l—k—1, 


and ), are the roots of the characteristic equation M(Q.)M*(\). The unknown quantities included in (20) are 
determined by substituting k(t) in (17) and in (12) and (13) and solving the algebraic equation (see [2]). Having 
thus determined the pulse transfer function k(t) which satisfies the requirements of dynamic accuracy and 
quality, we find the transmission function of a closed loop system ¢ (s),and next,by means of » (s) and the trans- 
mission function KWo (s) of the object,we find the transmission function KqWgq (s) of the correcting device 

by means of the following formulas 


making use of the methods described in [3, 6, 8]. 


KaWa(s) = 


= (k(t) (21) 
0 
® (s) 
285. (21 a) 
KW (s) (22) 
a (s 


2. Synthesis of Correcting Devices for Automatic Stabilization Systems 


Let us assume that one of the inputs of an automatic stabilization system is subjected to a random stationary 
disturbing action n(t) which has a known correlation function Rp(r) or spectral density S,(w) and the other input 
to a disturbing action which consists of a sum of a time function f(t) given in the form of a polynomial of power 


At) 


Fig. 3. Block schematic of an automatic 


stabilization system. 


At) 


A(t) 


(p), z,(t) 


fit} ou(t) Kg b{t) 1-Pip) 
Fig. 4 


t and a stationary random function u(t) which has a 
known correlation function Ry (r) or spectral density 
Sy(w). 


It is also assumed that the transmission function 
of the object is known and is equal to KgWo(s) (Fig. 3). 
Let us assume for the sake of simplicity that the correla- 
tion between n(t) and u(t) does not exist and that their 
mean values are equal to zero, 


The problem consists in finding by means of 
correlation functions Ry(r) and Ry(r), the given co- 
efficients of system Kj and the given transient state 
time T,the pulse transfer function k(t) which would 
provide a minimum value for the mean square error. 


Let an ideal system compensate the signal f(t) 


according to the equation 
h(t) = Hy (p) f(t), (23) 
where 
H (24) 
= 


Then, the error at the output of the system will again be of the form (Fig. 4) 


(t) = Hy (p) f(t) —\ n(t—) + de 


T 


kde @—t—0) + u(t —t—0)] b(6)do. 


—0o 


—oo 


By equating the mean value of the error to zero we obtain 


T 


Hy(p) f(t) =§ f(t —2) de + do = 
0 —co 


T r 


r T 


| = (t)— (28) 
n=0 
where 
&= (i =0,1,2,..,r), (27) 
ary 
en = (n =0,1,2,...,r). (28) 
Considering (24) and that f(t) is a polynomial of the power of r it is possible to obtain the equality 
la- 
” Considering Eq. (29) as an identity, it is possible to write r+1 limiting conditions imposed on the pulse 
transfer function by operator H fP)- These conditions can be determined in the forms 
K; = Yow (30) 
For the error, upon fulfillment of conditions (30), expression (25) takes the form: 
3) T 
k(e)dt (.u(t n(t k(n) de — 
0 —co 0 (31) 
By using above reasoning, it is easy to show that the necessary and sufficient conditions for obtaining the 
minimum of expression (31) and fulfilling the limiting conditions (30) consist in making the pulse transfer function ' 
k(t) satisfy the integral equation 
T T 
| Ra(t— 1) k()de | —0)b (0) do = 
@ . (32) 
= | Ru(t (0) dd + 
—oo —oo 
123 


oo 
r 
a) r r 


Considering notation (16) it is possible to re-write the integral equation (32) in the form: 


r 
(Ralt— +) + Res (t— de = Ru (t) + Det’. (33) 
0 


i=0 


Further, if it is again assumed that the spectral density which corresponds to the correlation functions 
Ry(r) and r) has the forms 


bot +,..+5,0" (@) M*(o) 


S () = Sn (w) + Sy = = (a) L* 


(34) 


it is again possible to write the solution of the integral equation (33) in the form: 


2k r q 
k(t) = Be™' + + (t) + 
j=o 


q 
+ 2 (t — T) + L(p) L*(p) (p) (p) [Raw 
0 


(0<t<T). (35) 


For the synthesis of the correcting devices formulas (21) and (22) can be used again. 


Example, It is given that 

Ry, (t)=0, (t) = Ae 8(t)= gt, Ry (t) = N% (x), H(p) =H, (p) =1. 
, For this example, the pulse transfer function which corresponds to the transmission function KgWo(s) has 
the form: 


b(t) = ae~* (a= 


Ry, = b(t) de 6) (6) d0 =Ntae~** > 0). 
The spectral densities which correspond to R,,(r) and Ryp(r) have the form: 


2A N*a* 
Hence 


2A Nig? _(2ABa* + + (248 + 


Whence 


M (p) M* (p) = (2ABa* + — (248 + pt 
L (p) L* (p) = pt — + p* + 


On the basis of expression (20) it is possible to write 


| - 
| 
| 


k(t) = + Bye 4 Ag + Art + Cod (t) + Dod (t —T)4 
N2 
+L (p) L* (p) (p) (p) | + 
It can easily be shown that a term of the form 


L* (p) (p) (p) [ 4 0 


and hence 


k(t) = Bye™! + Bye! + Ay + Ait + Cod (t) + Dod (t —T), (36) 
O<t<T, 
where 


By substituting the expression for the pulse transfer function in condition (12), (13), and the integral equation 


(17), we obtain six algebraic equations for determining By, Bg, Ag, Ay, Cy, and Dy. By solving these algebraic equations 
and inserting the results obtained into (36) we obtain the final expression for the pulse transfer function. 


By using formulas (21), (21a), and (22), it is possible to determine the transmission function of the correcting 
device. 
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ON THE SYNTHESIS OF PULSED CORRECTING DEVICES FOR SERVOSYSTEMS 


A. A. Krasovskii 


(Moscow) 


The optimal distributions are determined for the weighted coefficients 
of closed servosystems for the case when the signal is slowly varying and the 
noise is of an arbitrary stationary random nature, and also for the case when 
the signal is random and stationary, Formulas are derived for the values of 
the coefficients of a linear pulsed correcting device corresponding to the given 
weighted coefficients of a closed system. 


The use of digital controlling machines as links in automatic control systems has led to the recent appear- 
ance of a number of works devoted to the analysis and synthesis of some pulsed systems [1-7]. In works (2, 3, 5] 


is solved the problem of synthesizing a pulsed correcting device under the condition that some standard controlling 
actions be reproduced, 


The paper of V, P. Perov [5] contains the solution to the problem of determining the optimal weight function 
for a pulsed (sampled-data) system when the system input is made up of stationary random noise and a regular 
component of the controlling action. 


The present paper considers an analogous problem, but with somewhat different conditions imposed on the 
closed system's weight function and on the controlling actions, A simple method is proposed for deducing the 
equations which define the optimal weight function of a closed system, as well as relationships on the coefficients 
of the transfer function of the correcting device to be synthesized. 


Distribution of the Weighted Coefficients of a Closed System . 


The servosystem to be considered consists of a pulse element PE (Fig. 1) in series with a discrete (sampled- 
data) correcting device (DCD) and a continuous linear portion (CP), all shunted by a rigid feedback path. 


The pulsed (sampled-data) element, one of the so-called 


cp type 1 elements [8], transforms the continuous error signal, € = 
= Xin ~ Xout, to a series of rectangular pulses, Modulated in 
amplitude and with length equal to the repetition period. The 
Fig. 1. Schematic of a servosystem with discrete correcting device, in conjunction with the pulsed element, 
a discrete correcting device, implements an operation of the form 


where ¢, is the value of the input quantity (the error signal) at the beginning of the i'th repetition interval (r 
is the repetition time), 1; is the value of the output quantity during the i*th repetition interval and bj and a; 
are constant coefficients, to be determined by the synthesis of the correcting device, 
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Apart from digital computers, the operation given in (1) may be realized by means of various memory 
devices without the input and output being expressed in digital form, One of the possible designs for such a 
device is the magnetic drum scheme (Fig. 2) [7]. 


4 K 
7 
o- 


Fig. 2, Logical block diagram for a magnetic-drum pulsed correcting 
device for realizing the operation in (1), 1) A writing head, 2) a 
reading head, and 3) an amplifier. 


The system's continuous portion, CP, has the transfer function W, (in operator form) and xoy; =Wan. The 
input to the servosystem (the controlling action), xin, consists of a signal component, x, and random stationary 
noise § x, 


If we carry out the discrete Laplace transform on the variables 


bed 
k=0 k=0 
@o 
Zin = > Zin ete, 
k=0 
where z we obtain 
= Wis’ = Wi (in — ou 
Whence® 
d 
. 
out 1+ . 
Here, 
bi , bs 
+ Oy tt... + 
2” 


* In the case when the pulsed element's input function has discontinuities at times t= kr, one should use the 
more exact expression for the transfer function of the closed system which is given in (5). 


ing 
107 


is the discrete transfer function of the correcting device, W} is the discrete transfer function, determined in 
accordance with well-known rules (8), of the system's continuous portion. 


If the closed system's transfer function is expanded in powers of z~', the coefficients, c,, in that expansion 
will equal the values of the system's reaction to a unit pulse since, if 


WwW, -1 -2 


then the system's impressed motions are defined by the expression 


The sequence of coefficients, cy, cy, Cy, . . . , is called the distribution of the weighted coefficients, or 
the pulsed discrete transient response of the closed system. 


The desirable distribution of the weighted coefficients of a servosystem will be determined below from 
the conditions for the best possible reproduction of the input quantity's signal component for noise of given 
characteristics, 


It is first necessary, however, to formulate certain general requirements which must be met by the transient 
response of a closed system. 


1, For cg 0, a jump variation in the input would produce a corresponding jump in the output quantity, 
which is only possible if the executive organ has unbounded power. 


Thus, energy limitations require that one, or several, of the first weight coefficients equal zero, 


Further consideration is given to the case when cy =0 and cy # 0, It may be easily generalized to the case 
when an arbitrary number of the first weighted coefficients equal zero, 


2. It is well known [1, 2] that, in systems with sampled-data correcting devices, a finite duration of the 
transient response may easily be provided, i.e., a response in which all the weight coefficients, starting with the 
(M + 1)-st, are zero: 

CM4i = Cuya=...= 0. (4) 


With this condition, sometimes called the condition of an infinite degree of stability, the regulation time 
is rigidly limited, and equals 


T, = Mr. (5) 


3. If we set 6x =0 and 


= (to — t) = 2 (to) — (to) (t — to) + (to) (t 


— 


; i. 
Lin (h—p)= — + — + Sees 


= (kt). 


The 
| 
we obtain 
where 


ne 


If we substitute this expression in Relationship (3), we find that 


k = dot, + yx, + + ..., 
where 


dy = Cot Cy + + Cy, 
dy = — + + Bey... + Mey), 


dy= — (4+ 2+ 3°cg+... + (7) 


The quantities dg— 1, dy, dy, dg,... are the error coefficients of the closed servosystem. 
The condition for first-order astatism has the forms 


d,—1=0 


or 


Cot Cy + ey = 1. 
The condition for second-order astatism is 


d,—i=0, d,=0. 
The condition for q’th-order astatism is 


4. In addition to the conditions already stated, the condition must also be satisfied that the closed system 
is stable, 


It would seem, at first glance, that the holding of the condition of exact limitation of regulation time (4) 
would guarantee stability of the closed system. However, this is not the case, The necessary and sufficient con- 
dition that the closed system be stable is that the roots of the characteristic equation 

1+ =0 
lie in the unit circle, 


Optimal Distribution of the Weighted Coefficients 


With the above limitations taken into account, we may present the servosystem's transfer function in the 
form: 


(2) = + +... (8) 


where, for a system with first-order astatism, the weighted coefficients, cy, Cg, . . . , Cyy are connected by just 
the one relationship 


Cy + eg +... ey= 
Correspondingly, expression (3) takes the form: 


nt 
| 
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That distribution of the weighted coefficients will be called optimal for which the servo system's mean- 


square error, 
V Cour — 2) 


is a minimum, 


We shall determine the optimal distribution of the weighted coefficients for two cases: a) the signal 
component x of the input is a slowly varying function of times b) the signal component of the input is a 
stationary random function of time. 


1, If we set x;,, =x + 65x, where the signal component, x, is given by the power series in (6), we find that 


If the system has first-order astatism thep, for a slow change in x, we may limit consideration only to the 
velocity error, dyxy_ in reproducing the signal. In correspondence with this, 


Lout — Ze = + + +... + 


The dispersion of the error equals 
If we denote by R(Q)=Ry, R(r)=Ry, R(2r)=Rg,..., the values of the autocorrelation function of the 


stationary noise for the argument values 0, r, 2r,... , then we may express the dispersion of the error’s random 
component by the formula 


+ + = (C2? + 2+... + Ro + 
+ 2 + Coty +... + Ry + 2 + Cary 
+ Ry +... + 


When cy 1 =0, the conditions that a function of the variables cy, cg, ... , Cyy be 
minimum have the form: 


where ) is the Lagrangian multiplier, 
By carrying out the indicated differentiation, we obtain 


+ Cops) Rit + + + eyRy_, = — veds (2), 
dh, — + 204+... + Men, ay 


This system of linear algebraic equations defines the optimal distribution of weighted coefficients for a 
given correlation function of the noise and given average values of the squared derivatives of the slowly varying 


signal component of the controlling action. 
There is no particular difficulty in solving System (11). 


A, Thus, if the correlation time is less than the repetition period, more accurately, if R} =Rg=...=0, 
then [7] the optimal distribution of the weighted coefficients is the linear one: 


(z’) 
’ 


a 
— 
where 


(M + 1)* 

Ry 6a + 

2R, 
Om 

As an example corresponding to this case, the form of the correlation function, the optimal distribution 

of the weighted coefficients and the “envelope” of the transient response are shown in Fig. 3, a. 


Noise correlation Optimal weighted Optimal transient 
function coefficient distribution response 
a 
“sx av 
Rix} 
a 
7 at 
eT 
Ze, 
- 
av 


Fig. 3. Optimal weighted coefficient distributions and optimal transient responses, 
corresponding to various forms of the noise correlation function. 


With the given optimal distribution of the weighted coefficienrs, the dispersion of the error of the servo- 


system equals 
— 2) = Poly +35 


If the noise level is very high, such that a » M®, the given optimal distribution may be expressed as the 


uniform ones c, = = and 


| 


If the number of intervals, M, is large and a «< M® we obtain another uniform distribution: c, = ~ 


ARy 
and o= 4 


B. If the correlation time is so large that the correlation function can be considered constant during the 
regulation time, T,= Mr, 


then, in accordance with Eqs, (11), the optimal! distribution has the form: 
Cy = Cg = 


_ _M(M—1)—a 


The corresponding graphs are shown in Fig. 3, b. 

In this case, o = Rg, i.e., no “smoothing” noise is present. 

C. If the correlation function of the noise has a rectangular form, wherein the correlation time Ty, is 
greater than the repetition time 1, but less than the regulation time T,; > T, =Nr > r+ (N an integer), then 


the optimal distribution of the weighted coefficients consists of two sequences, one of which has the form of a 
step- function: 


a 2d 2d 
eu = — + = = = Cun = — 
ra 2d 4d 
Cu—n-2 = — + (M—1)=, CmM—N—3 == CM—N—4 = = => 
a 2d 6d 
CmM—2nN—3 = — + CmM—2N—4 = = CmM—3n~—3 = 


The corresponding graphs are shown in Fig, 3, c. 


2, We now consider the case when the signal component of the controlling action, just as the noise, is a 
stationary random function of time. 


In accordance with (9), the error of the servosystem may be given in the form 
Let — = + + + CM Le + + CoB + 


If it be assumed that the random functions x and §x are statistically independent then, if their auto- 
correlation functions are denoted, respectively, by Q and R, one obtains the expression for the dispersion of the 
error 


— 2) = (1 t+... + Qo t+ 2(— + + 
2 + Crem) — 2emQy + (2 + 02+... +04) 

2 + +... + Ry + 2 + Coa 
+ +... + 


With cy + Cg+...+Cyy=1, the condition for minimum mean-square error will have the form: 


=0 t,2,...,M) 
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or, in expanded forms 


Cy (Qo + Ro) + (Cv, + Cv41) R,) + (Qs + R,) + ees 
+ eu (Qy_, + Ru) = —+ + @,, 


CQ +eg+...+¢ey = 1 (v=1,2,..., M), (12) 


This system of equations defines the optimal distribution of the weighted coefficients in the case considered, 


For given autocorrelation functions Q and R, the solution of System (12) is not particuiarly difficult, 


Thus, if the correlation functions Q and R are rectangular, the correlation time of the noise is less than the 


repetition interval r, but the correlation time of the signal is comprised of N repetition intervals: 


Ryo+9, R,=R,=...=09, 
+0, Aner = =% 


then 


The optimal distribution of the weighted coefficients for the case M=10, N = 2, Q, = Rg, is shown in Fig. 4. 


The corresponding optimal transient response of the servosystem has an oscillatory character (Fig. 4). 


Fig. 4. The correlation functions of the noise and signal 
components of the input quantity (a), the corresponding 
optimal distribution of the weighted coefficients (b) 
and the servosystems transient response (c). 


Determination of the Coefficients of the Transfer Function of a 


Discrete (Sampled-Data) Correcting Device 


After one is given the optimal distribution of the weighted coefficients and, consequently, the transfer 
function of the closed system 


(2) = + -|- + ™ 
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the transfer function of the series- connected correcting device is determined from formula (2): 


Daly 
Whence, 
The transfer functions W{ and W$ are rational functions of the variable z~*: 
By by + bis +... + 
ata 
A, *+...+ 
The very simplest solution of Eq. (14) 
4 — — — —e,s 
is 
By = Aj +- +... + = ALO’, 
A; = By (A — — — . —cyz-™) = — ©). (15) 
The characteristic equation of the closed system has the form: 
A;A, + BiB; = 0. (16) 


With the coefficients of the correcting device determined in accordance with (15), this characteristic equation 


takes the form 


A\A, + BiB, = By(1 — ©’) A, + = = 0. 


It follows from this that the closed system to be synthesized by the method given will be stable if the 
discrete transfer function of the system's continuous portion has no zeroes or poles (with respect to z~+) inside of, 


or on the boundary of, the unit circle, 


In those cases when the continuous portion of the system does not satisfy this condition (is, for example, 
unstable or *neutral*), it is possible to use continuous correcting circuits which will provide the desired property 


to each individual continuous portion of the system, 


By equating coefficients of equal powers of z~!, we obtain the expanded form of Eqs. (15): 


m=p+M, 

by = 0, 

by = Aer, 

bg = Aol, + 


bm = Op Cu, 


n=q+M, 

Ay = Bo, 

a, = By — Bots, 

= Bs — Bic, — Pots, an 
= Bs — Bac, — — Bots, 


Gn = — 


| 


These relationships determine the coefficients of the discrete correcting device to be synthesized from a 
given distribution of the weighted coefficients of the closed system, and given transfer functions of its continuous 
portions, 


Example, As an example, we shall consider a system for the automatic control of an aircraft's dynamic load, 


It is known that, in longitudinal motion, the aircraft dynamic loading (without account being taken of 
gravitational force or of winds) is related to the angle of deflection of the elevator by an equation with an 
oscillatory term 


Y + 2Ewoy + of y = kod, 


where y is the dynamic load normal to the velocity vector, 6 is the angle of deflection of the elevator and 

we and 2£€ we are parameters related to the velocity, the aircraft altitude and to several other factors. It will 
be assumed that the control of the elevator angle is an electro-hydraulic executive mechanism with a rigid feed- 
back loop, described by an equation with an aperiodic term 


Tt kw, 
where u is the signal at the output of the correcting device, 
Figure 5 gives the block schematic of the control system for aircraft dynamic loading. 


PE Ke y 


Fig. 5 


The controlling action, consisting of the signal x and the noise § x, is added to the signal from the device 
which measures the normal dynamic load (accelerometer), The measurer of normal dynamic loading is taken to be 
inertialess, and is represented in Fig. 5 by the rigid feedback path. 


The error signal is applied to the pulsed element, to the discrete correcting device and then to the control 
of the elevator angle.The signal component of the controlling action is a slowly varying function of time whose 


dispersion, o (x*), equals 2sec™*, The experimentally determined correlation function of the stationary noise 
is shown in Fig. 6, 


a 
‘5 cy 
as} 
02+ 
as 
0 
Fig. 6. Noise correlation Fig. 7. Distribution of the weighted coefficients 
function. (a) and the transient response for the servosystem (b). 


On the basis of the design specifications, and the reaction of the system to step-function controlling actions, 
we may take the regulation time to be T, = Mr =1 second and the repetition interval to be r =0.1 second, By 
substituting these values in Eqs. (11) we find the optimal distribution of the weighted coefficients of the closed 
systems 


0.139, Cs = 0.240, C3 = 0.146, = 0.115, Cy => 0,200, 
= — 0,03; = 0.220, cg = —0.100, co = 0.160, = — 0,083. 
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The corresponding graph is given in Fig. 7. 
The discrete transfer function for the system's continuous portion equals 


where 


i 
— (—E + wo, = (— E— VEF— 1) woo, 


C1 = koi Tay Ca = Ca = 


A = (1 + Ts) + + 5*). 


In the case being considered, Relationships (17), which define the coefficients of the discrete correcting 
device, have the form: 


bo = 0. a = 0, 

= 0,139 Bi, 

be ~ 0.240 + 0.139 B, — 0.139 

bs -- 0.146 0.240 a, + 0.139 ag, a3 = Bs — 0.139 — 0.240 
big = — 0.085 as, ~~ 0.085 By. 


In particular, if w=5sec™*, =0.1, T=0,1sec and k=10, then 


by = 0.139, by = —0.044, by = — 0.134, ... , dig = 0,029, 
= 0.187, — 0.023, a3 = —0.01, ..., = — 0.0028. 


In principal, there is no difficulty in realizing these values of the coefficients of the discrete correcting 


device. The dispersion of the error in reproducing a given dynamic load, if one assumes an optimal pulsed correcting 
device is, according to formula (10), equal to 


— Z) = 0.21 (x’) + 0.59. 


If o(x")=2sec* then Volxg,, ~) w1. Thus, the mean-square error in reproducing a given dynamic 
loading in the given system is of the order of magnitude of 1. This best result is possible under the conditions 
given (a regulation time of 1 second, a mean-square velocity of change of the controlling action of 1.4 sec “1 
and a mean-square value of the random noise of 1.9, with the correlation function as given). 
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ON THE SIGNIFICANCE OF STABILITY CRITERIA FOR AUTOMATIC EXCITATION 
CONTROL OF FREE-~-RUNNING SYNCHRONOUS MACHINES 


L. V. Tsukernik 


Kiev 


The stability of automatically-controlled excitation of free-running 
synchronous machines, running in parallel and not connected to the rest of 
the energy system, is considered as a necessary (but not sufficient) condition 
for system stability, It is shown that excitation stability in this mode of 
operation may be used as an effective preliminary criterion of the per- 
missibility of coupling (for the analysis of the complex energy system's 
stability) of multi-stand machines, working in parallel, into one equivalent 
machine, 


In the development and investigation of automatic excitation control (AEC) for the Volga GES (State 
Energy System) there arose the question of the permissibility of replacing, for purposes of analyzing system 


stability, the proposed symmetrically-connected group of synchronous machines by one equivalent machine 
with the identical parameters and AEC. 


Tests on electro-dynamic model MEI [1] showed that the permissibility of such a replacement depends on 
the type of AEC, Then, in the Electrical Engineering Institute of the AN USSR[2, 3] and in MEI [4] there were 
developed methods for analyzing the stability of systems with AEC on several machines. The calculated results 
were in good agreement with the experimental data from the electro-dynamic models. 


However, the results provided by these workers did not give the necessary methodological recommendations 
for setting up designs of complex energy systems with AEC which would permit analysis of static stability accord- 
ing to complete criteria which would take into account the possibility of self-starting. 


The paper of V. M, Matyukhin [5] is a step forward in this direction, but it limited its consideration to 
ideal symmetric designs and multi-stand modes of operation, 


In actual conditions, the case is unavoidable when an unloaded machine operates in parallel with otheral- 
most completely loaded machines (immediately after synchronization). This spoils the assumption of ideal 
symmetry. Also possible are the cases when the circuit for switching the stands is non-symmetric, or when the 
parameters of the machine and the AEC are significantly different— and both the one and the other frequently 
occur. Under these conditions, the recommendations of paper [5] cannot be used, 


_ The problem of providing an admissibie equivalently connected machine stand for the analysis of an energy 
system's static stability may be significantly eased if it is realized that, as an initial mode of operation for such 
an equivalent, one may choose, as a preliminary, free-running machines working in parallel and unconnected 
to the external system, 


We consider now the equations for the disturbed regimen of a group of free-running machines connected 
to the machine-stand bus. For this purpose, we shall use Eqs, (1.16) and (1.25) of [2]*, in which the relative 
angles are calculated from the magnetic axis of the rotor of machine 1: 


* Or the identical Eqs. (4) and (5) of [3]. 
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<1 OP, 
(M mp* + + Ady + A£q = 0, (1) 
Fam Fam S8m-+ (RY? + p) + 
i=] 
(2) 
+d 104 + AEq — AEqnm =0. 
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We denote the machine stands by the numbers s = 2, 3,... , n+ 1 and, for machine 1, we assume a constant 
inertia My = oo and voltage V =const, Thus, machine 1 provides infinite power to the busses, the relative voltage 
vectors of which are all computed at the angle 6, We might note that the basic Eqs, (3) and (14) of (5) differ, 
for the particular case given (M, = 00), from the equations of [2] only in notation. 


Since loads and lines are not replaced in the design, it is completely permissible to ignore the active 
impedances of the machines and the transformers. Consequently, in the formulas for computing the partial de- 
rivatives p, 10), we shall have 54,70, ay, (k=s-1=1,2,...,m k j). 


If we assume that the AEC maintains an identical voltage on all the machines, we may set Ex, Eo) = 
= Vj. 


A number of terms in Eqs. (1) and (2) will be zero since, in the given case, 


oP, OV ay av, ax 0 
= = = 


If the automatic excitation control is implemented by varying current, voltage and their derivatives (AEC 
by varying angles and their derivatives is not employed in energy systems), it is then obvious that Eqs, (1) and 
(2) become mutually independent. 


Equation (1) defines only the disturbed mechanical motion of the rotorss* 


OP 
[Mup* + Dap + = 0. 
Equation (2) determines only the electromagnetic processes in the excitation circuits 


(Of + TY p) AEq — SE = 0. (2") 


For a free-running machine, the stability condition for Eq. (1*) is always satisfied. 


Thus, for the analysis of stability, it suffices to set the characteristic n'th-order determinant of Eq. (2") 
equal to zero, 


Ag — 0, (3) 


‘ Here, and in the sequel, we shall omit the unit in the subscript of 5 i, 80 that angles are computed with respect 
to the bus with infinite power: 64; = 5;. 
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= + TYP p— 1, ax; = Qs” + Tp (kj). 


As was stated, the equation A; =0 allows one to obtain a necessary, but insufficient, condition for system 
stability (2), p. 8). We now give this definition on actual physical meaning, We shall consider a necessary 
and sufficient condition for the stability of a free-running machine unconnected to the system as a necessary 
but not sufficient condition for stability of the system in the general case of all possible energy system modes 
of operation. 


Indeed, the free-running mode can occur due to an accident in the line, in particular, it all the line, or 
some portion of it, is single-circuited. The loss of stability in parallel operation of the machines of a multi-stand 
GES is inadmissible in this mode. 


It is of interest to note that the stability conditions for the closed-loop automatic excitation control of 
each individual machine are obtained from the characteristic equation 


+ TY p—1 =0. (4) 


It is obvious that these conditions must also be considered as necessary, but not sufficient. 


Very significant is the fact, disclosed both by calculation and by energy-system experience, that if ex- 
citation control is implemented by varying voltage, current and their derivatives, then the free-running mode is 
generally the most difficult, both as regards stability of excitation of individually-working machines and as re- 
gards stability of motion of machine-stands working in parallel, In the latter case, the assumption of the absence 
of connection with the external system also burdens the stability conditions with respect to machine operation. 
This may be explained as follows. If it is assumed that all the stand busses are infinite power busses then the 
reciprocal influence of the machines on each other's operation is completely removed and, consequently, there 
is no longer any question of the stability of relative motion. As the machine busses get further from infinite power 
(i.e., as the connection with the external system is weakened), the mutual influence of the machines on each other 
becomes less constrained by the system's reaction. Mathematically, in the denominator of the expression of the 
reciprocal moments of the machines there enters the reciprocal of the impedance, which increases monotonically 
as the resistance of the connection between the equivalent system stands and the busses of the grouped machines 
decreases, The same conclusion may be drawn with respect to the influence of the magnitude of the static load 
connected to junction points supplied by the group of machines, The free-running group of machines, unconnected 
with other groups (stands) must be considered as one of the initial modes of operation in the analysis of system 
stability. Therefore, in a number ot cases, there will be indicated variants which will be unsatisfactory from 
the point of view of machine motion stability in a given group of machines and with other initial modes of 
operation, With this the problem is eased, in that, for the free-running mode of operation, stability can be 
analyzed only with respect to electromagnetic coordinates (stability of excitation of synchronous machines), 

For this purpose, the domain of stability in the space of the chosen parameters is determined by means of Eq. (3). 


Thus, the necessary and sufficient condition for stability of an energy system can be decomposed into 
necessary, but not sufficient, conditions under particular conditions, which are best analyzed in a definite order, 


1. The stability is analyzed for a fully calculated design for a replacement energy system. In this design, 
each group of machines or stands is first connected to its equivalent machine on the basis of general considera- 
tions of design symmetry, degree of connection between the machines, type of AEC, etc, In the analysis, the 
system's initial mode of operation is varied, and the parameters both of the design and of the AEC are investigated. 


2. For optimal variants, in the case of new types of AEC, the stability is analyzed of automatic excitation 
control of free-running machines unconnected to the energy system for those groups (stands) earlier given in the 
form of equivalent machines, with respect to which nothing is given as regards stability of the mode of operation 
to be considered. For example, it is not necessary to analyze the stability of free-running stands of the receiving 
system, in which the proposed stand is connected to the new-type AEC being investigated. _ 


3. Stability of excitation is verified for the individual free-running machines using a new type of AEC, 
or for the corresponding groups of machines for grouped AEC, 


If the AEC of a machine is checked with respect to the natural mode parameters of each machine (for 
example, with respect to current or its derivatives), then a final check is made with respect to the motion of 
those machines which, in the computations of number 1 above, were joined to equivalent machines, With this, 
each symmetric group of machines in the design scheme is given as two identical machines (with identical AEC). 
The initial modes of operation correspond to the basic variations computed in number 1. 


Equations (3) and (4) used in the computations in numbers 2 and 3 above, were obtained without any 
assumptions as to the identity of the parameters of the machines and AEC, or as to the design's symmetry, The 
fact that the analysis of stability is carried out for the free-running mode of operation, without any connection 
to the system, significantly simplifies the computation. It becomes quite possible to carry out the design for 
systems of three or four machines without the use of complex computing devices, Consideration of the parallel 
operation of groups of identical machines in symmetric designs allows several important conclusions to be reached, 


We first substitute in the elements of the determinant in (3) the values of fh) and ri” and the corre- 
sponding partial derivatives which, for a free-running mode of operation under the assumptions made and for 
V;, =1, equal 


OE OE on Tek GE 9; Qi ) 


We obtain, 4s the result, 


z 
s [p. ar 4 — "ae (6) 
kj ak tk Ik dok P 


For a symmetric system of n identical machines in parallel and with AEC by internal current and by 
common voltage, it is easy to find the values for internal and mutual impedance of each machines 
n 
= Taj = (1) 


For a symmetric design, all the diagonal elements in the determinant in (3) are identical (aj, =a,) and 
all the remaining elements are equal to one another (ayy = a), After substituting Eq, ('7) into (6), we obtain 


a 


Expansion of the determinant with elements a, and a leads to the well-known expression [6] 


+ (mn — 1) a] (ag — ~. 0. (9) 
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After some transformation, we obtain from (8) and (9) the two characteristic equations: 


(10) 


—Taap —1 = + (n—1)a] = 0, 
FF; — 2a =: — a) = 0. (11) 


Equation (10) coincides completely with the characteristic equation of one real or equivalent machine 
working with an open stator circuit in a free-running mode, This equation depends on the parameters of the 
AEC with respect to total voltage and its derivatives, and does not depend on the parameters of the AEC with 
respect to current and its derivatives, nor on the reactive impedance of the machine. Conversely, Eq. (11) does 
not depend on the parameters of AEC with respect to voltage, but does depend on the parameters of AEC with 
respect to the current of each of the machines working in parallel (and also on their impedances, xq and xj). 
This equation corresponds to the operating mode of a three-phase short-circuited machine, 


Let us continue the analysis of Eqs. (10) and (11), If the transfer functions of the exciter and the AEC 


equal 


then Eqs. (10) and (11) take, respectively, the following forms: 
(1 + Tao p) (1 + Tsp) (1 + Top) pL =0, 
+ (1 + Top) (1 +71 p)— (kK pth’ p)=0. 


Here, T.= T 
d 


We now consider in more detail Eq. (14), which determines the admissibility of the equivalent machine 
connection. For regulation by current we obtain, from the condition that the free term in Eq, (14) be positive, 
the limitation of the coefficient k, known from the theory of compounding synchronous machines, for the con- 
dition of aperiodic stability 


In addition, we obtain from Eq, (14), by the Hurwitz criterion, a second limitation, defining the condition 
for the absence of self-starting: 


In the case when k* =k’ =0, Inequality (16) always holds. 
For positive values of k* and k* we obtain from Eq. (14) the following limitations, from the condition that 
the coefficients of the characteristic equation be positive: 


< + TT a’ +7 1T a] 


As was stated at the beginning of this paper, it was noted in running tests on the MEI electrodynamic 
model that, for AEC based on the first and second derivatives of the generator current, it was not possible to 
obtain stable parallel operation of the two machines which modeiled the generators ‘of the Volga és (1), p. 23). 


In this case, if we take into account the rigid negative feedback path for the exciter voltage, a path shunting 
the exciter and the regulator amplifier output and decreasing by a factor of five the values of Tp, T;, k", k’ and 
k (2), p. 13), we obtain from Eq. (17) 


kimax< Ts + T14+Ta') ta- an 


kmax 4800 andk 350. 


| : 

max < 74: (15) 


3) 


5) 


ion 


16) 


hat 


17) 


The limiting conditions for a loaded mode of operation required that k" be of the order of 10° and k’ of 
the order of 10°, Thus, the analysis of the free-running mode which we have been considering corresponds to 
experimental results, and allows the necessary conclusions to be drawn on the basis of uncomplicated computa- 
tions. 


It was also noted in the experiments that the unstable mode of parallel operation was observed particularly 
with respect to the reactive currents and, consequently, the reactive excitation of the machines, This is also in 
good correspondence with the relationships presented above. 


Instead of an AEC with respect to the current of each generator, one might use AEC with respect to average 
current (or total current on the outgoing lines), in which case the stability of parallel machine operation will be 
determined only by Eq, (10) or, respectively, by Eq. (13) since, in Eqs. (11) and (14), the terms depending on the 
AEC parameters will not appear. This implies the permissibility of the equivalent machine connection. 


If the automatic excitation control for each generator is implemented, not with respect to the voltage at 
a common bus, but with respect to the voltage at a point separated from the common bus by a reactive im- 
pedance x (for example, by a transformer impedance) then the right members of Eqs. (7) will contain, not Xq 
but q+ x). With this, Eq. (11) takes the following form: 


+ Fr) — (xa + 2) Ta, p — (2a + 2) = 0. (18) 


Equation (14) also changes correspondingly. 
Instead of the limiting inequality of (15), we obtain the followings 


kmax< + 2(1 +0). (19) 


Comparison makes it clear that, in the case considered, the condition for aperiodic stability is improved, 
since the allowable value of k,,,, is increased. The conditions for self-starting are simultaneously changea. 
However, with AEC on the derivatives, an increase in the value of kya, is particularly cogent, since this gets 
relatively large even for small values of x, In practice, the place to connect AEC with respect to voltage does 
not limit the admissibility of the equivalent connection of synchronous machine stands, 


SUMMARY 


1. The stability of an automatic excitation control system for free-running, parallel-operating, syn- 
chronous machines, unconnected to the rest of the energy system, can be considered as a necessary (but not 
sufficient) condition for stability of the energy system. 


2. The mode of operation referred to in number 1 above is defined in terms of a grouping, previously 
judged to be admissible, of connected synchronous machines (with small mutual impedances) into equivalent 
machines, 


3. The type of AEC, not the parameters and design of the machines connected, is of decisive influence 
on the stability of the excitation of parallel-operating synchronous machine stands, If the AEC is implemented 
with respect to total voltage and current of a group of machines with identical parameters, the replacement of 
this group by one equivalent machine is admissible, 


If the AEC is implemented with respect to the individual mode parameters of parallel-operating machines, 


the admissibility of connecting an equivalent machine must be verified, 
APPENDIX 
For each synchronous machine, we introduce the following notation for Eqs, (1) and (2): M and D are the 
constant inertia and the damping coefficient, P is the power (taken to be numerically equgl to the moment), 


Eg is the emf corresponding to the transverse synchronous impedance xqgj Fp, Fy, Fj, andF, are the transfer 
functions of, respectively, the exciter, and AEC with respect to voltage, to current and to angle, 
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The remaining nomenclature is in common use. 
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ON THE ERROR OF A LINEAR INTERPOLATOR FOR A 
PROGRAM-CONTROLLED DIGITAL SYSTEM 


Vv. V. Karibskii 


(Moscow) 


The principles of operation of a digital linear interpolator are described, 
and a theoretical investigation of interpolation error is made, A general ex- 
pression for this error is derived, and an estimate is made of the maximum 
possible absolute value of this error. 


Today, in many branches of technology, much attention is given to program-controlled systems which 
cause the controlled object to follow a given trajectory with high accuracy. The program for such a system is 
ordinarily given either in the form of a series of successive values of the trajectory coordinates, x;, Xg, ... » Xt 
Yi» Ya - - - » Yps OF, more frequently, in the form of a series of successive coordinate increments, Ax, AX, 

- » AXns Ay. . » AYp, Where Axk = — AYk ~ Yk- The control of the motion between 
the bracketing points is implemented by means of an interpolator in such fashion that the difference between the 
actual trajectory and the given one lies within the admissible limits, defined by a given accuracy, Moving 
successively from segment to segment, the object describes the entire trajectory. The program also gives the 
speed of translating the object along the trajectory, An example of the use of such a system is a program-con- 
trolled metal- working machine, in which the controlled object is a cutting tool. In the present work we shall 
consider linear interpolators which provide rectilinear motion of the object, with a given speed, between two 
points. For each k*th segment there is given the magnitude of the displacement, As, = ¥ Ax, + Ay; and the 
time of displacement whence the speed is V, = V3, + Vyk » Where Vy, = Ax /Tk and =Ayy/ Ty. 


The digital interpolator supplies signals in the form of control pulses; its problem consists of supplying, 
during time interval T;, the number oi pulses determined by the program for each coordinate. The distribution 
of pulses supplied to each channel during time T;, must be uniform in time. With this condition, the object will 
be translated along the trajectory in the way best approximating to the given straight line if a single control 
pulse causes the identical translation in each coordinate direction (identical scale of measurement). 


To carry out a uniform distribution of an arbitrary, previously given, number of pulses in the given time 
Ty, for different channels, is a difficult task and, in practice, an unsolvable one, Therefore, a device is used 
which does not provide a completely uniform pulse distribution during the time Ty. There is no technological 


difficulty in implementing such a device. Its schematic, based on the use of digital differential analyzer elements 
[1], is given in Fig. 1. 


The pulses with the fundamental frequency f max ate supplied from the main generator to the input of the 
first flip-flop, T1, of the binary counter. Before operation commences, all counter flip-flops are cleared to zero, 
In this counter, a differentiated pulse at the left anode of the flip-flop for each place (bit) is used to trigger the 
following flip-flop, while a differentiated pulse at the right anode of a counter flip-flop is applied to a coincidence 
gate, corresponding to the given place, and controlled by the memory cell register. 
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Fig. 1. Block schematic of a five-place interpolator. T,, Tz, ..., Tg are counter 
flip-flops, and G2ce are cycle-end gates; Gly, G5x, Gly,..., G5y are gates 
for, respectively, the X axis and the Y axis; MClx,..., MC5x, MCly, noe 
MC5,, are memory cells for, respectively, the X axis and the Y axis registerss 

f max is the repetition rate of the input pulses and 6 is the unit translation (value 


of a pulse division). 
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Fig. 2. a) The pulses applied to the gates 
from the counters b) the trajectory of the 
instrument, inertia being left out of account. 


As a result of the frequency division, pulses of 
frequency fe are applied to gates Gl, and Gly, of 


frequency fas to gates G2, and G2y, etc. These 


pulses are distributed in time in the fashion shown in Fig. 
2, a. The time of motion between two neighboring points, 
Ty» is determined by the number of counter stages and by 
the repetition rate, fax, of the input pulses. The fre- 
quency fmax is given by the generator, but the number 
of stages is varied by a code, giving the time Ty. If, in 
the circuit of Fig. 1, the code is given by a counter length 
corresponding to four stages, then only gate Glog” is 
opened, With this, when four stages of the counter have 
returned to the zero state, a termination pulse is applied 
to gate Gl.¢, this pulse passing through Gl., and inhibiting 
the application of input pulses to the counter, If gate G2c¢ 
is opened, but gate Glce is closed, the termination pulse 
is taken from the fifth stage when all five stages have re- 
turned to zero, and the time of motion is thus doubled. By 
means of an example, we shall now consider in more de- 
tail the operation of all these devices during the course of 
one cycle, 


* The subscript “ce* denotes the end of a cycle, while Glce is a cycle-end gate. 
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For simplicity of presentation, we shall take the magnitude of translation for one pulse to be 6 =1 cm, 
while Ax and Ay are measured in centimeters, Before motion begins, let the numbers Ax/5 =6 and Ay/é= 
= 3 be introduced into, respectively the X and Y registers. In the binary code (binary number system), these 
are written as 0110 and 0011. Let the cycle length be given as corresponding to four stages, This means that 
units ("1"s) will be found in memory cells MC2, and MC3x or the X register and MC3y and MC4y of the Y 
register, while gates G2,, G3, G3, and will be open, Asa result, in one cycle six pulses will go to the X 
axis and three to the Y axis, nonuniformly distributed in time, as shown by Fig. 2, b. 


For the circuit described, the magnitude of the translation, Ax, and the magnitude of the average velocity, 
Vx ay» are defined by the relationships 


Az >) 
i=o 
Vay = —max >) 2'Kix 
i=0 


where Kj, is 0 or 1, in correspondence with the binary code which gives Ax, n in the number of places (bits) in 
the binary counter, 6 is the amount of translation per one pulse and fj, is the repetition frequency of the 
pulses applied to the counter input, 


The magnitudes of Ay and Vy ay are defined 


" pulse No, of pulses 
per vole analogously. 
ee eee Teel ae It is clear from the description of the scheme shown 
Fig. 1 that the interpolation method used, by virtue of 
pulses) its discreteness and of the nonuniform distribution in time 
jet td 3 pulses) of the output pulses, cannot provide an output motion 
pulses) which strictly follows the straight line whose slope is given 
by At; As is clear from Figs, 2 and 3, the trajectory of 
4 motion is given by a zig-zag line, symmetrically oriented 
with respect to the desired line and joining the same two 
7 af end-points of the spatial interval given for the particular 
b 7 | Ay cycle, We now seek a general expression for the magnitude 
AS... \y of the deviation A of the zig-zag line from the required 
7” straight line (Fig, 3) at any discrete point of the trajectory. 
4 For this, we shall suppose that tan a o AL. = 1, since 
. ae si any zig-zag line which lies above the line for which 
tan a =1 has a line symmetric to it which lies in the 
Fig. 3. Trajectory of motion for Ax = 4 and region below this straight line (Fig. 2, b). 
Ay =3. 


In order to find the desired expression, we first de- 
termine how many pulses of the i'th frequency are distinguished in the N pulses applied to the counter, We shall 
number the frequencies as shown in Fig, 3, i.e., the frequency with which the counter puts out one pulse during 
the complete cycle we denote by 0, the frequency with which it puts out two pulses, we denote by 1, etc, 


It is clear from Fig. 3, that for the four-stage counter (n = 4), there are isolated, from the N pulses of fre- 


quency fmax applied to the counter's input, (+2) ip pulses of zero frequency, (Bed ip pulses of 
frequency one, (M22) ‘ pulses of frequency two, etc., where the notation “ip” denotes the integral part 


of the expression in brackets, Consideration of Fig. 3 and of the expressions given above allow one to note that 
the denominator within each bracket may be given in the form of the ratio 2";2!, where 2" is the number of 
pulses applied to the n‘th order of the counter during the complete cycle, and 2! is the number of pulses of the 
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i-th frequency isolated during the complete cycle, Thus, the number of pulses of the i-th frequency isolated by 
the n-stage counter from N input pulses may be expressed as follows: 


The corresponding total number of pulses given by the code may be expressed as 


Not = + Shp Ki, (2) 


i=o 


where the numbers K; take the values 0 or 1 and are the coefficients of the corresponding powers of 2 in the given 
number, expressed in the binary number system, i.e., Ky is the coefficient of 2, K, is the coefficient of 2, Ky 
the coefficient of 2%, etc. 


We now find the expression for A. On the basis of Fig. 3, the expression for A will have the following form: 


(3) 


4 = (y; — ys) cosa. 


We now determine the quantities which enter into Formula (3): 


(4) 


A 
y= 


We express Ay and Ax in the binary system, taking 6 =1, and we then find the expressions for the current 
coordinate x in terms of the number of pulses, N, impressed on the counter in the given interval of time: 


n—1 
Ay= 
i=0 


n—i 


At = >} 2! Kis, 
i=o 


Each pulse supplied to the X axis, or to the Y axis, gives rise to a translation of the executive organ by 
unit length. Consequently, the running value x equals the total number of pulses, Nx, isolated in correspondence 
with the code number for Ax, from the N pulses applied to the counter: 


n—1 
z=N,= [wae + Th Kix. (6) 


By substituting Expressions (5) and (6) into Formula (4), we get the result 


n—l 

> 

Ki. i=o ip 


The expression for y, can be written analogously to the expression for x : 
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We now express cos a in terms of Ax and Ay: 


n—1 


2 2' K,.. 


By substituting (7), (8), and (9) in Formula (3), we obtain the expression for the error in this method of 
interpolation: 


2K nai" +5] K,.— [ + | K 
2 w 2 Jip ix 2 iz ty 


aw (10) 


where V = 1,2,...,2"—1. 


The formula just obtained permits the computation of the error for any value of the ratio at and for 
any discrete point of the trajectory for any number of places, n. 


For the selection of the value of the pulse division * in constructing a digital control system with the aim 
of achieving a given accuracy, it is important to know the maximum error that may be expected from the 
method of interpolation, as defined by the given scheme for a given number, n, of places. For this it is desirable, 


for fixed n, to determine, not only the magnitude of Amaxe but also the value of ~L and the number N for 


which the error is maximal. Moreover, it is of interest to determine the dependence of Ang, on n, and the 
limit of Amax a8 0 — ©, if this exists. An exact solution of these problems seems very complicated, so for a 
beginning we shall try to give at least a rough estimate of the magnitude of Amax.- 


of 4 Ke an reducing by uy each tem of [ Ky where 


we start with the inequality, N2—" + 5 —1<[ <N Then, Formula (10) 


takes the forms: 


= au) 


—i 


If we replace the letteri by j 2 KyandY (N an — 5) Ky and 


carry out algebraic transformations in the numerator, we obtain 


* The magnitude of the translation for one pulse. 
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Ben.) + (S24) 
finally, 
Consequently, 
Ky) cosa + (3 Kin) sin al, (a3) 


The quantities > Kx and Kin are the sums of the digits in the binary expressions for Ax and 
imo 


Ay, and take their maximum value, equal to n, when all the Kj, = Kjy =1. For a rough upper estimate for 


n—1 n—1 
Amax: We can assume that the sums >) Ki, and >) Kj, retain their maximum values as a varies, 


Then, Formula (13) may be written in the form 


nex max (cos + sina). | (14) 
The sum cos a +sin a has its maximum value for a = 45°. If we bear in mind that A may also have a 
negative value (if, in Formula (9), the subtrahend is increased, while the minuend is decreased by unity), we 
obtain, finally, 
| 4 max! <7 (15) 


It may be shown that this estimate is very rough, since the dropping of the restriction to integral parts of - 
numbers in Formula (9) leads to a substantial coarsening of the inequality. 


In reality, this is not the case, It is possible to cite several unfavorable cases for deriving n, when the 
error differs from the maximum by a factor of 2,5 to 3, and increases with increasing n, With this, it should be 
remembered that, in actual devices, n is determined by starting from the required values of Axra, and Aymaxe 
since AxXmax * AYmax = 1. 


As an example, we shall consider in more detail the case when n is even, Ax is given by the code , 
101010, i.e., Kj, =0 foreveni and Kj, =1 for odd i, and Ay is given by the complementary code ,... 
010101, i.e., Kiy =1 fori even and Kiy =0 fori odd. In computing the error A for different values of n al 
this case, we may note that A will have its maximum value for N=Ay. 


Bearing this fact in mind, let us determine the dependence of 4 on n in the particular case when N= Ay, 
using Formula (10) for this, To do this, we consider the dependence of each individual sum entering into (10) 


on ns: 
n—1 n/a 
i=o 


(2° + 23) + 2° + 2%) 4 (29-4 2" 4 28 4 25) + (294 27 4 294 284 


Here, = 24 + is a geometric progression with initial term a» = 2 and ratio q= 44 


the sum of this expression is given by 


— — 


a—l 


i—n i = 
+3), = OF 14+ 5424485 4 


(29 2) $ (204 294 28) (PH dy (7) 


k=l 


n—1 


P 2+ 3 Fhe Ky = 


i=o 


n—1 
Ay = Ky = 2+ P+ (19) 


15) 


(20) 


731 
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If we substitute the expressions just obtained into Formula (10), and carry out the computations for arbitrary 
n, we obtain, for the given case, 


+A 


3V5(2"—1) 3\ 2 3 3V5(2"—1) 3\ 3 2 


A= 


n n A n 


In analogous fashion, we may show that, for the case when Ax is given by the code... 111111 and Ay 
by the code , . . 1010, with n even, the error has its maximum value for N = Ay, and is approximately 0.835 nA, 
For the case when Ax is given by the code ...10101 and Ay by the code ...01010, then at the points N= 
= Ay, for odd n, we shall have 


bin n n 4 


For 00 


SUMMARY 
1, With an increase in Ax;ax = AYmax = 2"— 1, there occurs an increase in n, the number of counter 
stages, which leads to an increase in the absolute error for some values Ax and Ay. 


2. The absolute error cannot exceed 7° 


3. As shown by a number of particular examples, the absolute error may increase proportionally ton and, 
n n 
in certain cases, reaches a value from ove @ Ova" 
4, Formula (10) allows one to compute the error for an arbitrary case and to estimate its maximum possible 
value, but it does not answer the question as to which case is the very worst. 


5. It should be remembered that a number of combinations of Ax and Ay are favorable. With these com- 
binations, the error does not depend on n, the number of counter stages, and always remains less than unity. 
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INVESTIGATION OF ALGEBRAIC EQUATIONS USING 
ANALOG COMPUTERS 


Vv. M. £l'yasberg 


(Leningrad) 


A simple method is presented for the investigation of algebraic equations 
on analog computers, the method being based on the reproduction of the given 
polynomial as the solution of the differential equation defining it, The method 
permits the immediate determination of the real and imaginary roots of al- 
gebraic equations, the approximate estimation of the value of their complex 
roots, and also the investigation of dynamic systems for stability, using the 
Mikhailov criterion. 


The methods given in the literature for solving algebraic equations on analog computers are based on the 
investigation of dynamic systems, the characteristic equations of which are the given algebraic equations (see, 
for example, [1, 2]). 


The present work proposes a method for investigating algebraic equations which is based on the use of 
analog computers to obtain graphs of the given polynomials, 


Basic Method 


Let us consider a method for obtaining the graph of a given multinomial on an analog computer, 
Let Yy(p) be an n-th order polynomial in p : 


Yn(p) = + + p™*+ +++ + + ao. 


It is simple to show that the polynomial in (1) is the solution of the differential equation 


(p) _ nt (2 
dp” 
with the initial conditions 
Y (0) 
Y’ (0) = 1! a,, 
Y"(0) = 2! ay, (3) 


Y (0) = kl a. 


The solution of Eq. (2) can be obtained by means of an n-fold integration with the initial conditions in (3), 
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Since time is the independent variable in an analog computer, solution of Eq, (2) on such a computer 
necessitates the substitution: 


p=mt (4) 
p=—mi. (4a) 


or 


The solving circuit will consist of n integrators in series, 


If substitution (4) is used, the graph of the multinomial in (1) will be obtained for some interval of varia- 
tion of the independent variable, 0 to + po, while if substitution (4a) is used, the interval will be from 0 to — po. 
By varying the scale of the transformation, one may vary the size of the interval of solution within which the 
size of the polynomial, Y,,(p), and its derivatives do not go beyond the limits permitted by the given concrete 
computer, 


Determination of an Equation's Real Roots 


If the multinomial in (1) is the left member of an algebraic equation 


Yn (p) 0, (5) 


then, in obtaining the graph of this multinomial, one may determine the real roots of the equation as the abscissas 
of the points of intersection of the curve Y,(p) with the axis of abscissas. 


In the solution of the problems of automatic control, one has to do, in the majority of cases, with algebraic 
equations with positive coefficients. As is well known, such equations cannot have positive real roots. Therefore, 
for the determination of the real roots of the characteristic equation of an automatic control system, it suffices 
to use the substitution in (4a), 


The accuracy of the solution is determined by the technical capacity of the model (the analog computer). 
Thus, in solutions effected by the author, of seventh-order equationsoncomputers MN-7 and IPT-5, the error in 
determining the real roots did not exceed 5%, 


As an example, let us take the solution of the equation 


p® + 44.62 p* + 733.5 p® + 5590 p* + 19670 p + 23120 = 0. 


The roots of this equation, found analytically by an iterative method, are 


hs = 17.75. 


— 2.25; — 6,504 — 11.62; 


To obtain the graph of the multinomial Y,(p) in the left member of Eq. (6), we carry out the substitution: 
p= -— 2. Then, in correspondence with (3), the initial conditions are defined as follows: 


| 
| 
Fig. 1 


a) 


Pos 


er). 


6) 


ons 


conditions 
Y (0) = 23420 46.2 
Y’ (0) = 1! (—2)-19 670 = — 39 340 78.7 
ad (0) = 21 (—2). 5 590 44 720 89.4 
m0} == 3! (—2)*-733,5 = — 35 210 70.7 
= 17 130 34.3 
— 5! —2\8 = — 3840 7.68 


The scale of ¥ (p) was sy = 500 units/volt, 


The equipment used for the problem on analog computer MN-7 is shown in Fig. 1, while the graph of the 
multinomial Yg(p) is: given in Fig. 2. 


Fig. 2 Fig. 3 


Figure 3 gives the graph, obtained from the computer, of the multinomial Yg(p) in the left member of the 
equation 


p® + 47.7 p® +- 1105 p* + 14370 p® + 97590 p* + 302900 p + 302500 = 0. (1) 
The roots, found analytically, are 


— 1.8; = —6.0; = —7.2; Y= —13.6; — 9.64/14, 


As is clear from a consideration of the oscillogram provided, the values of the real roots obtained from the 
computer are sufficiently accurate for many applications. 


The given method is also applicable in the case of multiple real roots. As is well known, if A, is a root 
of a polynomial with multiplicity k, then it is also a root of all 
derivatives of this polynomial up to order k~ 1 inclusive, The 
solving circuit permits the graph of all derivatives of a polynomial 

Pw) to be found simultaneously with the graph of the polynomial itself, 
By using this property of multiple roots, one may easily determine 
the multiplicity of a real root in every case, 


a By using the principles mentioned one may, in addition to 
finding directly the real roots of an algebraic equation, also de- 
termine the imaginary roots, find approximate values of the com- 
plex roots and investigate the stability of a system wing the 
Mikhailov criterion. 


Fig. 4 
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Obtaining the Mikhailov Hodographs and Determining Imaginary Roots 


In the multinomial in (1), the left member of a system's characteristic equation, let us substitute jw for p- 
We may then write 


(jw) = P(w) + (0). (s) 
With this, the multinomials P(w) and Q (w) are alternating polynomials of the following form: 
P (w) = — + — + — ages, + (9) 
Q(w) = — age? + — + — (10) 


If we make the substitution w = m;t, and use the method described above, we may obtain on the computer 
the graphs, simultaneously, of the polynomials P(w) and Q(w) or the characteristic curve, P= f (Q). Thus, the 
given method allows us to use both forms of the Mikhailov criterion. 


The graphs of the polynomials P(w) and Q(w) may be used in determining the imaginary roots of the 
equation. 


When there are pure imaginary roots, of the form A, yy, = + jwy, the multinomials P(w) and Q(w) 
reduce to zero simultaneously, From the magnitude of w for which these multinomials go simultaneously to 
zero, one may judge the magnitude of the imaginary root. 


Figure 4 shows the graphs of P(w) and Q(w) for the equation 


p® + 18.7 p* + 474 p® + 3.710-10%p? + 5,22-10*p + 7.82-10* = 0, (11) 
The roots, found analytically, are 


hy = —1.65; = —8.53 + 14.1. 
It is clear from this oscillogram that the imaginary root is determined with sufficient accuracy. 


Approximate Determination of Complex Roots 


If the imaginary part of a complex root is less, in absolute value (modulus) than its real part, or is some- 
what larger than it, the graph of the multinomial in (1) has a minimum absolute value, The abscissa of such an 
extremum correspondsapproximately to the real part of a pair of complex conjugate roots and, from the mag- . 
nitude of the ordinate, one may frequently judge approximately the magnitude of the imaginary part [see Eq. 
(6) and Fig. 2). 


Fig. 5 Fig. 6 


If, on the other hand, the imaginary part of a pair of complex conjugate roots is much greater in absolute 
value than the real part, then the magnitude of the imaginary part may be obtained, approximately, from the 
graphs of the polynomials P(w) and Q(w). To illustrate this assertion, Fig. 5 gives the graphs of P(w) and 
Q(w) for the equation 


p® + 17.6 p* + 453 p® + 3.38-10%p* + 5.17-104p + 7,82-108 = 0, (11a) 
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(11a) 


with roots 
hy = — 1.65; dos = + 0.554 713.2; 


and Fig. 6 shows them for the equation 


p® + 19.6 p* +. 495 p® + 4.04.10%p? + 5.27+104p + 782-104 = 0, (11b) 
with roots } 


hy = — 1.65; hes = — 0.55 +- 13.2; hes = — 8.53 + 714.1. 


These equations differ from Eq. (11) only in that the pair of roots, \g s, in addition to the imaginary part, 
has a relatively small real part, positive in the first case and negative in the second, 


As is obvious from the oscillograms, the imaginary part of the roots may be approximately determined as 
the mean arithmetic yalue of w for which the graphs of P(w) and Q (w) intersect the axis of abscissas (of course, 
only in the case when the points of intersection of P (w) and Q(w) with the axis of abscissas are sufficiently close 


to one another). Moreover, it is possible by this means to judge roughly of the size of the real part of the com- 
plex conjugate roots. 


For a more exact determination of the complex roots, one may use the method based on the construction 
of a family of Mikhailov curves, and described in [3]. 


Special Characteristics of the Scale Transformations 


The scale transformations must be such that in the desired interval of solution, the values of the polynomial 
itself, and of its derivatives, do not exceed the limits permitted in the analog computer i At the same 
time, the requirements on the accuracy of the solution must be fulfilled. 


The domain of variation of the values of the polynomial and its derivatives in the required interval of 
argument variation (and, in the majority of cases, this interval as well) are not known beforehand, Therefore, 
it is impossible to give a general dictum on the choice of the scales, However, as practice has shown, the choice 
of the scale factor in any given concrete case poses no particular difficulty. We give now some recommendations 


for choosing the scale factors when using standard analog computers, of the types MN-7, IPT-5, MPT-9, etc., in 
a number of characteristic cases. 


1, If the boundary of the region of root (real) distribution is previously known (the magnitude of the 
largest root is approximately known) then the time scale must be so chosen that the duration of the time interval 
of the solution is not greater than 10 to 12 seconds, Experiments have shown that, in this case, the absolute 
values of the initial conditions may be taken sufficiently large to provide satisfactory accuracy. Moreover, with 
such intervals of solution time, the error due to the zero drift of the dc amplifiers is small. 


2. If the region of root distribution is not known beforehand (the usual case), itis then necessary to choose, 
as a first guess, a time scale such that the initial conditions, for all orders of the derivatives, are approximately 
of the same order of magnitude, The solution which then follows will show whether or not the time interval 


within which none of the polynomial derivatives exceeds the magnitude imposed by the computer does contain 
the entire region of root distribution. 


This is easily determined by a comparison of the graphs of all the polynomial's derivatives. If the given 
interval does not contain this region or, conversely, if it covers the region a little too well, then the time scale 
should be changed correspondingly. 


3. If the roots of the equation are widely separated from each other then, to obtain the graph of the poly- 
nomial over the entire range of its roots, several scale transformations must be made, each of which allows one 
root, or a group of roots (real), to be determined, However, in this case, the most efficient procedure, after 
having found one root, or several roots, is to reduce the degree of the polynomial by a division, thus easing the 
choice of the proper scale factor for determining the remaining roots. 
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4. The Y scale factor must be so chosen that the largest initial condition corresponds approximately to 
the largest possible value of voltage in the given computer, If, with this value, some derivative exceeds the 
permissible limits, then the Y scale factor may be changed, necessitating a proportional change (in this case, a* 
decrease) in the initial conditions for all the derivatives. With this, the value of the smallest of the initial 
conditions must constitute no less than 1% to 2% of the maximum value allowable in the computer. 


Similar recommendations for the choice of scale factor may also be given in other particular cases. We 
might mention that the change from one scale factor to another is very simple, since it is not required to re- 
calibrate the equipment, but only to vary the values of the initial conditions, 


The establishment of a coefficient of an equation on the computer is done only by establishing given 
initial conditions on the corresponding integrator, and may be easily varied, Therefore, the use of the given 
method is convenient for the investigation of equations where one coefficient, or several coefficients, are varied 
(for example, for the study of questions of root migration), 
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A DEVICE FOR SOLVING HIGH-ORDER ALGEBRAIC EQUATIONS 


Jifi KryZe 


(Prague) 


The paper describes a high-speed, completely automatic, device for 
simulating algebraic polynomials and solving high-order equations. The 
device is constructed on a new principle, and does not contain electronic 


vacuum tubes in the simulating (analog) portion, thus removing the basic 
source of errors, 


Simulating a Complex Variable 
Any device designed to solve algebraic equations must simulate, in a definite way, polynomials of the 


form 


W=X4+j¥ ae, 
i=0 


where z is a complex variable and the aj are the polynomiai's coefficients. In the general case, these co- 
efficients are complex, 


The device must determine those values, z;, (k=1, 2,..., 0), for which W =0, Le., 


X=0, Y#=0. (2) 


In order to simulate the polynomial in (1), use was made of damped oscillations of tuned RLC circuits, 


Fig. 1. The connection of oscillatory loops to obtain independent 
damped oscillations, 


Figure 1 shows several such circuits, all connected in the anode circuit of a vacuum tube, If the tube is 
conducting, a current I flows in its anode circuit. After a large negative voltage is suddenly impressed on the 
tubes grid, the tube is cut off and, in the loops connected in its anode circuit, there arise, independently of one 
another, free damped oscillations, described by the equations 


| 
Uy, 


up, w,t. (3) 


where UL is the instantaneous value of voltage in the i-th loop, Lj and Cj are the values of inductance and 


capacity in the i-th loop, wj = 7, is the frequency to which the i-th loop is tuned and 9; is the damping 
coefficient, defined by the expression 


(4) 


(Q; is the Q-factor of the i-th loop). 
By means of the phase-shifting loops, one may obtain from the voltage UL the voltages 


u, = uy = sin wit, (5) 


where U is some initial amplitude. 


If these voltages are applied, respectively, to the hori- 
zontal and vertical deflection plates of a cathode-ray tube, the 
z face of the tube will show the logarithmic spiral, represented 
in Fig. 2. With this, the quantities u, and uy may be considered 
as modelling the complex variable z: 


motion of point z in the complex plane, if z varies in accord- 


z=2-+ jy, (6) 
; y 4 m3 and the spiral of Fig. 2 may be considered as the image of the 
ance with the law 


z= = (coswt + jsinet). (1) 


Thus, the voltages ux and uy will be the images of x 
Fig. 2, The logarithmic spiral on the and y if9=9; and w=wj. 

face of a cathode-ray tube, correspond- 
ing to the motion of the point z in 

the complex plane. 


Consequently, by means of the damped oscillations one 
may simulate (model) a complex variable which varies in 
accordance with the law of the logarithmic spiral, With this, 
the point z “searches” the area of the complex plane swept out by its motion. The area of this "search" is 
defined by the relative distance between two neighboring arms of the spiral: 


z(t) 


i= z(t) 


3 


2a, (8) 
| 

It is possible to observe the entire spiral of z on the face of the cathode-ray tube. But the electron beam 
may be passed only at those moments of time when z is a root of Eq. (1), i.e., when the voltages modelling the 


real and imaginary components of the polynomial in n (1) are, in accordance with (2), equal to zero. Then, 
phosphorescent points will appear on the face of the tube, these points defining the distribution of the roots, z = 


= in the complex plane* (points zg, Zs on Fig. 2). 
* This principle for showing roots was also used in [2-4]. ’ 
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In order to determine the roots of Eq. (1), it is necessary to model, not only z, but also the quantities z!, 
If z is defined by means of the time parameter t , in accordance with Eq, (7), then the i-th power of z can be 
expressed in the following way: 


at = = Zi iwt sin iwt). (9) 


It is easy to convince oneself, by a comparison of Eqs. (9), (7), and (5), that the i-th power of z may be 
modelled by means of the i-th loop, if this loop is so tuned that the following conditions are mets 


=io, = ih, (10) 


where w=uyj and $=9; are, respectively, the frequency and damping coefficient of the first loop, modelling 
the first power of z. 


Consequently, to model an n-th degree equation it is necessary to connect n loops in the circuit shown in 
Fig. 1. These loops mtst be tuned to the frequencies w, 2w, 3w,..., nw. The Q-factors of these loops, 
according to Eqs. (10) and (4), turn out to be identical: 


@, 
This is certainly very convenient from the builder's point of view, 


Formation of the Terms a;z! 


If z is defined by means of the parametric representations in (7) and (9), for simplicity we take Z;,; =1, 
and the coefficients of the polynomial are expressed in polar form 


a; = A,ei*, (12) 
then the polynomial in (1) takes the form 
ian ian 
W = >} [cos (iw t + + j sin (iwt + (13) 
i=o 


Consequently, to model the individual terms of the polynomial, it is necessary to have voltages proportional 
to the real and complex components of these terms: 


Re (a;z') = Aye cos(iwt'+- Im (ajz*) = sin (iw t + 9). (14) 


As is clear from (14), these voltages must be shifted in phase with respect to the voltages of the tuned loops, 
and must have initial amplitudes which are proportional to the absolute values of the coefficients, 


In order to obtain such voltages with controlled phase and amplitude, a small take-off coil with three inter- 
laced windings, 1, 2, and 3 (Fig. 3), is placed in the coil of each oscillatory loop. This take-off coil can rotate 
about two mutually perpendicular axes. The field in which the take-off coil is placed may be considered homo- 
geneous. However, the relationships given below remain valid even for inhomogeneous fields if the take-off 
coils are wound in a certain way, for example, as shown in Fig. 4, The number of turns of such a winding, per 
unit length along 0, is constant, 


Axis A, about which the bearing supporting axis B rotates, is perpendicular to the direction of magnetic 
induction. The planes in which take-off coil windings 1 and 2 lie are mutually perpendicular, and intersect in 
axis B. Windings 1 and 2 are each divided into two identical portions, 1a, 1b, 2a, and 2b. Winding 3 is per- 
pendicular to windings 1 and 2 and to axis B. Axis A lies in the plane of winding 3. 
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We may assume, for simplicity, that the number of turns in the individual windings, 1a, 1b, 2a, 2b, and 
3, are such that voltages identical in magnitude are induced in each of these windings, with the condition that 
each winding's plane is turned to a position perpendicular to the direction of magnetic induction, However, the 
assumptions made are not necessary, and the results which will be obtained below may be easily stated also for 
windings with different numbers of turns, different forms and different diameters. 


Fig. 3. a) A schematic representation of the winding distribution 
of the take-off coil with a coiled spherical form; b) the mounting 
of the take-off coil with three windings (represented by one turn) 

in bearings. The state is that for which aj =0. 


Fig. 4. Completed take-off coil 
winding (schematic) for which all 
relationships hold, even in a non- 


homogeneous field. 


ul 


We introduce the following notation for the take-off coil in the 
i-th oscillatory loop, whose frequency is wj=iw: y; is the angle by 
which axis B is turned from the position for which the plane of winding 
2 is parallel to the direction of magnetic induction (for oj = 0); this 
angle equals the angle formed by the plane of winding 1 with axis A; 
oy is the angle by which axis A is turned from the position in which 
winding 3 is perpendicular to the direction of magnetic induction; 
Ugas Ub» Uga, Ugh, Ug are the instantaneous values of voltage in the 
individual windings; Uj, Uj, and uj are, respectively, the vector, 
amplitude and instantaneous value of voltage in any winding of the 
take-off coil in the i-th oscillatory loop in the position in which this 
winding is perpendicular to the direction of magnetic induction. 


The voltage u; varies according to the law 


u; = U; sin wt, = (15) 
where Uj inj is the amplitude at the initial moment of time. 


The individual windings are joined to the phase-shifting loop in accordance with Fig. 5, a such that the 
polarity of their voltages satisfies the expressions uy, = — Uyb, Ugg = — Ugh (the measurements are made with 
respect to the terminal common to all four windings). The resistances and capacitances in the phase-shifting 


loops satisfy the relationship 


oR,C; = 1. (16) 


The instantaneous values of voltage on the take-off coil's windings are: 


Ujo= — = Uy COS Ugg = — Ugs = Uy SiN SiN aj, Us = Uy COS (17) 
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After some elementary computations, we obtain the formulas for the instantaneous values of voltage at 
the terminals + Re and 4 Im of the phase-shifting loops: 


= ini Ky cos a; sin (wit + + doi), 
Une = Ui jim Ki cos a es! cos (wt + di + doi), (18) 
= —Uim, U—Re = —URe, 
where 


Ro(1 +i 


By comparing Eqs. (18) and (14), we find that 


Ure = Re (aj2*), = — Re (a;2'), = Im(ajz*), = — Im (aj2') 
for 


ini Ki cosa, = Ay, $i + Gor = (19) 


Consequently, by turning the take-off coil about axis B, one can establish the argument , of the coefficient 
aj; by turning the same take-off coil about axis A, one can establish the amplitude, Aj of coefficient aj. 


~ n-ne 
? 
4G, 
c d 


Fig. 5. a) The se: for obtaining the voltage correspond- 
ing to the term ajz* (Z is the total load impedance); b) 
is the equivalent total impedance; c and d) the loads on 
winding 3. 


However, the circuit as given would have one essential disadvantage. It is clear from the diagram of the 
phase-shifting loop in Fig. 5, a that the take-off coil, with the loop connected to it, draws active and reactive 
power from the oscillatory loop. This power depends on the coil's position. A rotation of the coil would have an 
effect on the magnitude of uw; and Q;. This would result in an essential distortion of the results. 


The apparent power required of the phase-shifting loop is 


= Ut a4 (7 — jw C1) = Ut cos* = Ut costa, . (20) 
(%— 720; 
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This expression is obtained by the symbolic method from well-known formulas, It is a consequence of 
Eq. (20) that the apparent power depends not only on ~; but also on Zj. The same power would be required by 
connecting the voltage Uj cos aj to a total impedance Z,j (Fig. 5, b) consisting, for example, of a resistance 


R= a in parallel with a capacitor C = 2C; or a resistance R= 


in series with the capacitor C = 4C; 
(U, is not the effective, but the maximum value, of voltage). 


According to Expression (16), there is a voltage us = uj sin aj on winding 3. If winding 3 is loaded with 
total impedance Zo; (Fig. 5, c ord), then the apparent power drawn by the third winding will equal 


Sy Ut — jw C,) = Ut = U? sina, (21) 


i 
C; 
The total apparent power required by the take-off coil will equal 


S=Sy+5,= (22) 


and will not depend on Z; or aj. This has the consequence that a rotation of the take-off coil will not influence 
the load on the oscillatory loop and, consequently, will not affect the quantities wj and Qj. 


This last circumstance is very important, since it allows one to draw a significant amount of apparent power 
from the oscillatory loop (more than 1 watt) and to obtain a large voltage from the take-off coil (more than 100 
volts), 


In order to obtain the real, + X=+ Re(W), and imaginary, + Y =+ Im(W), components of the polynomial 
W, it is necessary to add the voltages on all terminals, + Re and + Im, corresponding, respectively, to the real 
and imaginary components of the individual terms of the polynomial, The polynomial's free term, ay , is 
modelled by adding in dc voltages proportional to Re (ay) and Im (a 9). 


For the summation of the voltages from the individual phase-shifting loops, it is possible to use any of the 
well-known adding circuits built of linear passive elements, The input total impedances to the adding circuit 
will be comprised of the total load impedances, Z;, shown in Fig. 5, a. 


For determining the instant when the quantity z passed through the value of a root, a discriminator is 
necessary. This must develop a pulse which gates the electron beam only when the quantities X and Y become 
zero simultaneously or, what is completely equivalent, when the largest positive (or negative) of the quantities 
+X, — X, +Y and — Y equals zero, Since such electronic circuits are familiar, they will not be described 
here. 


Over-all Block Schematic 


Figure 6 gives a block schematic of the device. The circles denoted by the aj represent the take-off 
coils, The block labelled £ is the phase-shifting and adding loops, while the discriminator is denoted by D. 
Deflection of the electron beam along the spiral z is implemented by means of a voltage induced in a sub- 
sidiary winding mounted on the coil of the first oscillatory loop, and shifted in phase by means of RC elements. 


Advantages of the Device 


The device operates at high frequencies. The basic frequency, — , is best chosen so that its value lies 
between the limits of 20 to 50 kilocycles. 


The repetition frequency of solutions is about 40 cycles, With this, the spiral z will have rolled itself 
up to a fraction of a percent of its original diameter, None of the devices hitherto known can implement such @ 
rapid repetition rate of solutions, In the highest-speed device [3], solutions are repeated at the rate of one per 
second, 
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The high frequency of solution repetition is very useful, since it permits observation on the cathode ray 
tube's screen of the inanner in which the roots of an equation move in the complex plane when there is a smooth 
variation in the equation's coefficients. If a cathode-ray tube with long persistence is chosen, one can observe 


the entire curve describing the roots sought. This latter property is particularly important in the analysis of 
stability of automatic control systems. 


= 
7 


Fig. 6. Block schematic of the device, 


The high-speed automatic device considered does not have vacuum tubes which form the voltages re- 
presenting the polynomial values. Therefore, variation of the parameters of vacuum tubes does not affect the 
accuracy of the device, except for the zero drift of the discriminator, But, thanks to the high values of the 
input voltages, even the requirements on smal! discriminator zero drift are not particular stringent, Thus, the 


number of vacuum tubes is relatively small, and does not depend on the degree of the equations for whose solution 
the device is designed, . 


Analysis of Some Sources of Error of the Device 


Cathode-ray tube, The first source of possible error in the device is possible nonlinearity in the cathode- 
ray tube, and inaccuracy of readings on its screen. This error may be easily eliminated, A pulse is impressed on 
the cathode-ray tube network both for the roots of the equation W =0 and also for the root of the equation 
Zg ~~ a=0. Then the representation of the one point, z, =a, is obtained on the screen. 


By varying the value of a, one can assure that the point corresponding to a recuts successively all the points 


corresponding to the roots of the equation W=0, The corresponding values of a_ are also the roots, zy, of this 
equation. 


In view of the fact that the coefficient for z»_ is constant, it suffices, in forming the subsidiary equation, to 
make direct use of the voltage used for deflecting the electron beam in describing the spiral z, The formation 
of the term a is implemented by means of two dc voltages analogously to the formation of the term a, in the 


polynomial W. For the solution both of the equation W =0 and the subsidiary equation, one discriminator (with 
switching) can be used, 
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Frequencies and Q-factors of the loops. A second source of possible error in the device is a nonconstancy 
of frequency or Q-factor in one or more of the individual loops. In the process of solution the vacuum tube V 
is cut off and, therefore, the inconstancy of its parameters has no effect of the device's operation. Thus is 
avoided a source of instability which is fundamental in electronic generators. 


A tuned loop is, in itself, a stable element, and the stability of its parameters can be provided by good 
mechanical construction, 


Moreover, it is clear from the equations provided earlier that, for accurate operation of the device, there 
is required, not so much the maintenance of absolute constancy of the frequency and Q-factors of the individual 
oscillatory loops, but rather the maintenance of the accurate ratio between the individual frequencies and the 
equality of the Q-factors of all the loops, The absolute value of the frequency affects only the correct operation 
of the summing R,Cj loops and, therefore, it is completely satisfactory if the frequency oscillations do not exceed 
+ 0.1%, 


The error due to faulty maintenance of the mutual relationships of frequencies and Q- factors of the in- 
dividual loops accumulates from period to period. Therefore, the accuracy of solution may be increased by 
making the substitution z'= kz, The coefficient k is so chosen that, if possible, the roots of the polynomial, 
after the substitution has been made, all lie on the outermost arm of the spiral. 


Computation shows that if the mutual relationships of the frequencies are maintained with an accuracy of 
5*10°5 and the relationship of the Q-factor with an accuracy of 3- 107°, then the error in determining a simple 
root will not exceed 1% of the spiral'’s maximum diameter. 


Fig. 7. Image on the Fig. 8, Image on the 
cathode-ray tube screen cathode-ray tube screen 
obtained in testing fre- obtained in checking the 
quency. Q- factors. 


For verifying the accuracy with which the frequency and Q-factor relationships are maintained in the 
device, a simple control is provided. For example, to verify the accuracy of correspondence of the frequencies 
in the first and the i-th loops, one solves the equation 


a,z' = 0. 


With this, the inputs + Y to the discriminator are switched off, The polarity at the discriminator output 
will be positive (the cathode ray appears) if the following relationship holds: 


X = Re(W) = Re (ajz*) = Re = 0, 
i.e., for 


cos (iw t + = 0, 
iot+q= (k=0,4,...,2i), 
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A star is obtained on the screen (Fig. 7), made up of 2i straight line with identical angles between them, 
If the frequency, w;, of the i-th loop is not exactly i times the frequency of the first loop, spirals will appear 
on the screen instead of straight lines, This is explained by the circumstance that the passage through zero of 
the sinusoid sin wt (which controls the moment at which the ray appears) will be translated by an ever-increasing 
distance in one direction or another relative to the passage through zero of the sinusoid sin wt (which controls the 
position of the ray), By rotating the star by means of a change in the argument, gj, of the coefficient aj, one 
may verify the correctness of the scale of y;. By interchanging the inputs X and Y, it can be verified if the 


voltages at these inputs are shifted with respect to one another by = » etc, 


Verification of the equality of all the Q- factors is implemented by solving the equation 


t 

(z— = 0. 

The solution is carried out with either + X or + Y switched off from the discriminator input, The screen | 

| 


will show the same star as in the previous case, but shifted by the amount b with respect to the first (Fig. 8). 
Incorrect tuning of the Q-factor of the i-th loop leads to spirals, instead of rays, on the star. 


EXPERIMENTAL RESULTS 


A model of a device for solving third-order equations, based on the principles described above, was con- 
structed in the Laboratory for Automation and Remote Control of the Czechoslovakian Academy of Science, 


Fig. 9. Circuit for the phase-shifting and summing networks. wRC =1, 
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Testing of this device verified the advantages of these principles; in particular, the stability of the relative 
frequencies and Q-factors of the oscillatory loops was verified, as well as the possibility of drawing significant 
voltages from the take-off coils. 


The phase-shifting loops were connected in one circuit (Fig. 9). This circuit, for each harmonic in- 
dividually, was identical to the loop shown in Fig. 5, a. Since the coils of the oscillatory loops were implemented 
in the form of simple cylindrical solenoids without shielding, the interaction between loops gave rise to noise. 

It was, therefore, only possible to obtain results with an accuracy of about 10%, It appeared, with this, that it is 
necessary to devote a good deal of attention to the discriminator construction, so that this latter might be able 
to produce pulses no longer than 0.1 microsecond, 


The device operated with a solution repetition rate of 50 cycles, The frequency of the fundamental 
harmonic was 50 kilocycles, The peak reactive power of the oscillatory loops was approximately equal to 1 
kilovolt-ampere, while the output voltage of the circuit in Fig, 9 was about 100 volts. The device contained 
six vacuum tubes, An ordinary oscilloscope was used to indicate the points corresponding to the roots of the 
equation to be solved, On the basis of these experimental results, work is presently being conducted to develop 
a more highly perfected device, ; 
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INVESTIGATION OF THE INFLUENCE OF VARIATIONS IN MEMBRANE 
CURVATURE AND EFFECTIVE SURFACE ON THE OPERATION OF 
PNEUMATIC DEVICES CONSTRUCTED ON THE PRINCIPLE 

OF FORCE COMPENSATION 


Chou Ching-lieng 


(Moscow) 


An investigation is made of the effect of deviations in effective surface 
and curvature of rubber-tissue membranes on the quality of operation of 
pneumatic regulating devices constructed on the principle of force compensa- 
tion, It is shown that deviations of membranes’ effective surfaces from their 
nominal values can lead to a disruption of regulator adjustment and to the 
appearance, in isodrome regulators, of residual irregularities. 


One of the basic elements in the construction of modern pneumatic devices designed for inspecting and 
regulating productive processes is the rubber-tissue membrane, Thus, for example, all the blocks of pneumatic 


unified plant systems (UPS) are constructed with membranes (see [1-3]). The quality of the performance of such 
devices depends heavily on the membranes and their characteristics. 


There are works in the literature which are devoted to investigations of membranes (see [4, 5]), but no 
work in which questions having to do with the effect of various individual membrane parameters on the quality 
of performance of pneumatic devices is considered, In particular, in works dealing with methods for designing 
devices constructed on the principle of force compensation, the important question of how to determine the effect 


of variations in the membrane's effective area on the course of these devices’ basic characteristics is left un- 
touched, 


Observation of the operation of control devices which contain membranes shows that variations in the 
membrane characteristics, related particularly to changes in their effective areas, have an essential effect on 
the course of these devices’ characteristics, From the point of view of the operating conditions for a regulator, 


the estimation of admissible deviations in membrane effective area from the nominal value has an important 
practical significance. 


Deviations of the effective area from the nominal value of membranes prepared from the same specifica- 
tions and from the same material may arise due to the following reasons, 


1. Deviations in the geometric dimensions of the membranes which are within the tolerances permitted 
in their preparation. 


2. Differences in material properties of membranes of one and the same brand ("trade-mark") arising 
from inaccuracies in the preparation of their bases and irregularities in their rubber films, 


3. Installation conditions, on which depend the initial positioning of the rigid membrane center with 


respect to its sealing surface and with respect to the nozzle, i.e,, the initial membrane flexure and, hence, its 
effective area, 
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4. Variations in the pressure drop supplied to the membrane, 


Variations in a membrane's effective area arising from the first three causes are constant for each con- 
crete membrane, i,e., they do not change during operation of the membrane. 


Deviations due to pressure drops vary in magnitude since the drop varies continually during operation of 
the membrane. 


Vaph a The present paper presents a portion of the work 
Ym) OPER being carried out on membrane investigation at the 
= \N Pneumo-Hydro- Automata Laboratory of the IAT AN 
SSS SV 5 SSSR, Its aim is to estimate, by computation, the effect 
Priv Wy of variations in membrane effective area on the course 
P arom ; 2 S of characteristics of regulators constructed according to 
WANN 
zag? inZ the principle of force compensation, Here also are 
rie treated questions dealing with the effect of membrane 
GY ‘out J curvature on the operation of certain pneumo- automatic 
4 Y devices. 
Y 2 The investigation of the effect of variations in 
membrane effective area was carried out using the 
= —* 4 example of the regulating block for a UPS, as shown 
Supply in Fig, 1, which was chosen as being one of the most 
Fig. 1. Sketch of a UPS's control block. characteristic designs for devices of this type. The 


investigation was carried out for the cases when the 
deviations of the effective area were constant in magnitude, and when they varied during membrane operation. 
We determined the influence of variations in membrane effective area on the magnitude of the regulator's gain 
and on the magnitude of the control point's * displacement or, what amounts to the same thing, on the magnitude 
of the regulator's residual irregularities, 


Le Investigation of the Effect cf a Constant Deviation of the 


Membrane's Effective Area 


Derivation of the basic formulas, We shall use the following notation: Pgiv, Pms» Pout and Pp, are 
the excess pressures in the chambers, in accordance with the nomenclature used in Fig. 1, i.e., the pressures 
corresponding to the given and measured values of the controlled parameter, the pressure at the regulator output 
and the pressure in the positive feedback chambers; AP=Pms — Pgiv is the error pressure at the regulator in- 
puts Fy, F 9, Fs, F4, and F, are the nominal values of the membranes’ effective areas (Fig. 1); for the calculations 
we take F, =F, =F, =F, =Fy and F;=mFy (for the control block of a UPS system, m=0.6)3 
5F, and 5F, are the absolute values of the deviations in effective area of the corresponding membranes from 
their nominal values. 


For the computations, we assume that the delivery characteristics to be regulated, and the constant values, 
are linear. Then, for a static regulator, Gy=c (Pp, — Pg) and G2=cz (Pour — Pp), where G, and G, are air 
delivery (mass-second), respectively, through valves 7 and 6, cy and cy are the delivery coefficients for these 
valves, and Pg is the pressure in the isodrome chamber. The pressure Ps, constant in magnitude, can be established 
in those cases when, for an isodrome regulator working as a static regulator (i.e., for the isodrome switched out), 


* The term “control point,” frequently used in describing regulators, is understood to mean the given value of 
a controlled parameter which must be maintained by the regulator. An isodrome regulator must always 


maintain a parameter at its control point, so that a constant pressure is established at its output, corresponding 
to the object. 


The control point displacement is the magnitude of the error, + AP, between the measured and given 
values of the regulated parameter, for which the constant pressure is established at the regulator output. 
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it is required to establish some definite pressure Po, yt the regulator output for zero value of error, AP, at the 
input (control point).* The magnitude of this pressure is established arbitrarily by completely covering the iso- 
drome valve at the moment when the pressure at the regulator output and in the isodrome chamber equal the 
required value of pressure, Pout, (with this, the hermetic sealing of the isodrome chamber must be maintained). 


From the condition of equality of the discharges, G, = Gg, in steady-state operation, we determine the 
value of the pressure, Pp, in the positive feedback chamber 


Po = Py, + A(Poy, — Ps), (1) 
where 
= 
Amada’ 


We now set up the equation for the forces acting on membrane block 1, and the regulator's static equation, 


a) With the nominal values of effective membrane areas, the equation for the forces acting on block 1 has 
the forms: 


Pry Paint — Ang) + — Pour = (2 


By substituting Eq, (1) in Eq. (2), we obtain the expression for determining the value of pressure at the 
regulator output: 


Pout = (Pms— Poiy) + Po- (3) 
The expression for determining the regulator's gain has the following form: 


Pourdn_ 1—™ 


where Poyr, is the output pressure for zero error pressure ‘at the regulator input. 


In Formula (4), the subscript *N" on K and Poy;, indicates that these values are chosen for the nominal 
values of membrane areas. 


b) When account is taken of the deviations of the effective membrane areas from their nominal values, 
the equation of the forces acting on membrane block 1 has the form: 


[Ps + A (Pour — (Fn + 1) — + 
+P Pg) (mF y+ 5) + + — Pour (Py + = (5) 


It follows from this equation that 


Pout = (P_.—P..)+ 
ad 
(Fe Poiyt (1 + (6) 
F, + 


* Such an adjustment of a static regulator, ordinarily called control point tuning, has the aim, for nominal load 
on the regulator, of maintaining the parameter exactly at the given value, and the regulating organ is thereby 
found at the necessary position (usually at the center). 


| 
il 
t 
| 
1 
| 
ons 
lished 
ut), 
ig 
751 


We find from Eq. (6) that the gain equals 


N 
K= (7) 


We now write the expression for determining the relative deviation of the gain: 


Figure 2 contains graphs illustrating the variation in the magnitude of AK as a function of the various 


magnitudes of relative deviations, and , of the effective membrane areas from 
their nominal values, The design characteristic for various adjustments of the controller is defined by the quantity 
A. The relative deviations, for computational purposes, were chosen between zero and 20%, and varied both on 


the side of increasing, and on the side of decreasing, area. 
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In computing the graphs of Fig. 2, as with all following computations, it was assumed that a deviation of 
effective area occurs only in the given membrane, the area deviations of the other membranes being taken equal 
to zero. 


It is clear from the graphs of Fig. 2 that variations in the effective areas of membranes 1 and 5 have the 
greatest effect on variations of the relative gain, AK, this effect being the greater, the larger the magnitude of 
A, i,e., the larger Ky. 


In view of the fact that, in the case considered, the values of 5F;, SF, 5Fs, 5Fg, and 6F, are constant, 
the deviations of AK engendered by them will also be constant and may therefore be compensated for by 
corresponding adjustments of the valve range so that they do not affect regulator operation; however, with this 
the scale of valve range will be confused. 


The regulator shown in Fig, 1 can operate either as a static regulator (whereby the isodrome valve is 
completely closed) or as an isodrome regulator, it being more frequently used in the latter capacity, The re- 
gulator is so constructed that, when the isodrome is switched in, the regulator possesses no residual irregularities 
and the static regulation errors equal zero. 


However, as will be shown below, such errors can arise in the regulator due to variations in the effective 
membrane areas from their nominal values. * 


This may be explained in the following manner. At the end of the regulation process implemented by an 
isodrome regulator, the pressure in the upper chamber of the isodrome element is compared with the pressure in 
the regulator's output line, P,,,;. With this, the pressure in the lower chamber of the isodrome element and the 
pressure in the positive feedback chamber also become equal to Poy. If there were in fact an equality of the 
effective membrane areas, F; =F, =F, =F, =F, then the regulator, brought to an equilibrium state by the 
pressure equalizations cited above, would exhibit the condition, Pojy=P ms. If these effective areas were not 
equal among themselves, equilibrium could only be established n Pgiv #~ Pimms . In this case, the regulator 
operates as a static regulator with a small astatism. 


In order to determine the possible magnitude of the residual irregularities, we substitute the pressures 
= Pour = Ps in Eq. (5). We then obtain an equation of the following form: 


After transforming this equation, we obtain an expression for determining the magnitude of the residual 
irregularity: 


giv 

AP = (Pros— Paid = Pom — 7 — (10) 
(1—m— + (1—m— *) 


It is clear from Eq, (10) that if the deviations, §F, and 6F,, of membranes 2 and 4 do not equal one another, 
but the deviations, 6F, and 6F,;, of membranes 1 and 5 are equal, or do not exist, then the regulator possesses a 
degree of irregularity, AP, which is constant in magnitude, If there do exist deviations, unequal in magnitude, 
of the effective areas of membranes 1 and 5, then the residual irregularity, in this case, will be a function of the 
regulator's output pressure, i.e,, of the load, Since the deviation, 5Fs, of the area of the flat central membrane 3 
enters only in the denominator of Eq, (10), it may not lead to the appearance of residual irregularity. 


* Apparently, this is precisely the reason why mass-produced isodrome regulators are given a certain tolerance 
on the displacement of the control point, 
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2. Investigation of the Effect om Regulator Operation of Variations in 


The effect of a variation in effective area of some membrane or other on the magnitude of the residual 
irregularity is shown graphically in Fig. 3. It is clear from these graphs that a very insignificant deviation of a 
membrane's effective area may lead to the appearance of a large residual irregularity. 


Since the deviation of effective membrane area which is permissible from the point of view of residual 
irregularity is significantly smaller that the deviation giving rise to inadmissible variations in the gain, the gain 
will automatically be kept within permissible limits if the minimum amount of residual irregularity is held to, 


So that the residual irregularity not exceed the established value of 1% to 1.5% for pneumatic isodrome 
regulators, the deviation in effective area of one membrane, when no deviations exist in the other membranes, 
must not exceed 0.5 to 0.6% of the nominal value. Since the deviation in area of one membrane may be com- 
pensated by the deviations of other membranes, the actual permissible deviation may be a bit larger than 0.5% 
but, in any case, it must be quite small. This bespeaks the particular importance, for the accurate operation of 
pneumatic isodrome regulators constructed on the principle of force compensation, of carefully choosing mem- 
branes, of equal effective areas, and carefully installing them. 


Fig. 3 


Effective Membrane Area Due to Variations in Supplied Drops 


In this case, the deviations of the effective membrane areas from their nominal values are functions of 
the pressure drops acting on some men.brane or another, and may be described in the following manner: 


0 Of 06 06 kg/cm 
~Q0/ 
~Q03 
-004 
-005 
1\ 
oF, 
N 
Fig. 4 


Fy = /1(Pout), 


IN 


Fs 
(11) 


= f,(AP) for Pyiy = const and = fs (Pour )- 


An experimental investigation of membranes* provided 
characteristics expressing the dependence of the membrane's 
effective area on the pressure drop, both for flat and for crimped 
membranes, Figure 4 shows one of these characteristics, 


Sel = f (AP), taken off for a flat membrane with a caprone 
base. As is clear from the figure, this characteristic is nonlinear. 
If, instead of the relative deviations in the effective 
membrane areas, a » we substitute their values, f; (AP), in 

Eq. (6), we obtain a following equation: 


* This investigation was carried out by the author in conjunction with G, T. Berezovets at the Pneumo-Hydro- 


Automata Laboratory of the IAT AN SSSR, 
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Pout (1 —A— (Pour) + fs (Pour)) = AP [1 — m — fy (QP) + + 


+ fa + Ps (4 — A) + fi (Pour (12) 


It follows from (12) that the regulator gain, and the relative deviation of the gain, are expressed, respectively 
by the formulas: 


1 —m + AP [f, (AP) — f, (AP)] + fe(AP) — fa (AP) + (AP) 
1 A+ Pout [fg(Pout) Af, (Pour + fo(Pout) — 


(13) 


1—A 
AK = —— K —1, (13a) 
1—m 
where f{ (AP) is the derivative of the function f; (AP). 
If these equations are used with the assumption that all the membranes are flat, and the computation of 


AK is made by substituting the values of f(AP) and f*(AP) which are taken from the graphs of Fig. 4 for the 


corresponding pressure drop AP, then the curves of Fig. 5, showing the variation of AK as a function of AP, may 
be obtained, 
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Fig, 5, Curves of the variation of the relative deviation of the gain. 
a) For a deviation of effective area of membrane 1; b) the same 


for membrane 5; c) the same for membrane 4; and d) the same 
for membrane 3. 


Before considering these graphs, we should note that the various membranes may have drops of differing 
magnitude acting on them, Since there is atmospheric pressure on one side of membranes 1, 2, 4, and 5, these 
membranes are always subject to a one-sided pressure drop which, due to the actual range of pressures in the 


regulator, cannot be less than 0.2 atmospheres. The center membrane 3 can be acted on by very small drops, 
acting on either side, 


We now consider the individual graphs. The dependence of AK on A, i.e., on the adjusted range of re- 
gulator valving, exists only when there are deviations of effective area of membranes 1 and 5. As before, the 
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5 effect of these deviations is the greater, the larger the magnitude of A, i,e,, the larger is Ky. For membrane 1, 
X with the output pressure Poy; varying only between the limits of 0.2 to 1 atmosphere, the magnitude of AK 
cannot exceed 5% for each of the values of A given on the graphs. But for membrane 5, with the same limits 
4 of variation of P,,,,,, the size of AK increases rapidly with increasing A, For large A, the variation of AK can 

4 exceed 50%, The variation of AK brought about by a deviation of the effective area of central membrane 3 

H cannot exceed 14% for error signals in the range AP =0 to 0.4 atmospheres. 


The deviation of AK brought about by a deviation in the effective area of membrane 4 can, for small 
values of Pyiye reach a large magnitude (greater than 20%). 


1 We now investigate how a deviation of effective area which depends on pressure drop effects the magnitude 
} of the residual regulator irregularity, By substituting Expression (11) in Eq. (10), we get 


[fs (Pout) — (Pout)! (fe (AP) — fa (Pay 
AP = Fins — Fyiv = + Ja — fo (AP) Tm BP) 


4 In the equation just obtained, a deviation of effective area of the individual membranes shows the same 
} influence on the magnitude of residual irregularity AP as the analogous deviation of effective membrane area 
which entered into Eq. (10). 


Thus, it is clear from the results obtained that, for accurate operation of isodrome regulators, it is necessary 
that the deviations in effective membrane area, both those which are constant in magnitude and those which 


depend on pressure drop, be quite small. 


3. Estimate of the Effect of Membrane Curvature on the Accuracy 


of Operation of Pneumatic Devices 


The nozzle and the flapper are two of the most important elements of pneumatic devices. The flapper, 
as a rule, is bound either to the rigid center of one membrane or to the block which connects several membranes 
and, therefore, the flapper movement is determined by the movement of membranes, 


Since, during operation, a membrane can move with curvature, 
the correct operation of a nozzle- flapper element may be disturbed, 


A Wy HZ. since the flapper will not be established parallel to the nozzle face. 
z 

ss (Ss A skewing of a flapper with respect to a pete due to the 
VIEL curvature of a membrane, is particularly undesirable for certain 
toe pneumatic devices, since the accuracy of operation of the device 

2, depends on this. 

As an example, we might cite the multiplier-divider scheme, 

A, shown in Fig. 6, for which a detailed description was given in work 


Fig. 6 [6]. In this device, curvature can have a particularly adverse 
. 8: effect on the quality of operation. We shall show how this is 
indeed the case. 


The basic equation of the multiplier- divider is 


P,P 


. where P, is the output, and P;, P,, and Ps are input pressures, this equation being derived under the assumption 
- that the flapper which, in this case, simultaneously covers two nozzles, moves with respect to each of these by 
one and the same magnitude, h. 


This condition, on which depends the accuracy of the device's operation, will hold if the flapper face and 
the two nozzles are all parallel to each other, and if the motion of the flapper is strictly perpendicular to these 
faces. Curvature of the membranes can violate this requirement. Indeed, if the membranes are so canted that 
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the plane of the faces of the flapper is established at some angle to the plane in which the nozzle faces lie, then 


the distances between the flapper and each of the nozzles will no longer be identical, but will be some hy and hg 
(Fig. 7). 


We now consider in more detail how a possible curvature of the 


membranes can affect the accuracy of the operation of multiplication, 
defined by Formula (15). 


YY Yy, In the multiplier-divider, both the constant and the variable valves 
Ve must have linear characteristics. If this condition holds, we may write the 
following equilibrium equations for relay chambers B and C: 


Fig. 7. Positioning of the 


fla with th 
K;(P3— K,Pz, = KyP,, (15) 


is curved. where Ky, Kg, Kg, and Ky are coefficients of proportionality of the delivery 


write 


drops at the individual valves (corresponding to Fig. 6). 
Let us assume that Kg=K,=K. Bearing in mind that the pressure P, accurately follows pressure Py, we may 


K K 
Po= 37K, Pe (16) 
For the operation in (15) to be executed accurately, it is necessary that K; = Ky, 1.e., equality of the pro- 


portionality coefficients of the two valves which, as was already stated, cannot always be attained due to 
membrane curvature, 


If we denote by h the distance from the flapper to the point E, lying midway between the nozzle centers, 


then the quantities hy and hy will equal: 


h=h+atga, h—h—atga, 


where a equals half the distance between the nozzles or, alternatively, 


h,=h—atga, h=h+atga 


depending on which side the membrane is canted, 


We now express the coefficients K, and Ky as products of the discharge coefficients by the cross-sectional 


area of the corresponding valves: 


K, = pyrd,hy = + atga), Kg = = pnd, (h —a tga), (17) 


where dy and dy are the-diameters of valves 1 and 4, 


If we substitute Eqs, (17) in (16) and carry out the corresponding transformations, we obtain 


In order to determine the relative error, we put Expression (18) in the form 


Py = +9). 


The relative error, 6, can be found from the formula 
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Let us assume that pgd, = 14d). In this case, the formula for the relative error has the form: 


B=: (1— tga (20) 


We now estimate the order of magnitude of 6. We take for the nozzle, Ky = j1y7 dyhy, letting a = 10° 
(tan « = 0.0029), a = 2 mm and hy =0.05 mm. We consider the mode of operation corresponding to the mean 
point of the device's scale, for which P; =0.5 Ps. In accordance with Eqs. (16), in order that the equation P, = 
=0.5 Ps hold, it is necessary that Ky = Ks = K (we ignore the effect of air compressibility), i.e., that pgm dghy = K, 


With these conditions, 


= —1 =13.1%° 


—i= 


2a 2 
tga K 1 —2 0.0029.0.5 


We chose h, = 0.05 mm as a rough, order-of- magnitude, figure. 


The investigation provided above shows that curvature of membranes with flappers attached can show a 
significant affect on the accuracy of pneumatic devices similar to the multiplier-divider considered above. For 
membranes used in such devices, even small curvature is inadmissible. 
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* In the experiments conducted by the author with blocks which included three membranes, careful placing of 
the membranes led to minimum angles of cant of the membranes" rigid center, of the order of 10°, The experi- 
ments, described in [6], on an assembled multiplier-divider gave high accuracy in performing the operations 
considered, This indicates that, in the model used for the experiments, the flapper slant was less than given 
above, Moreover, it should be remembered that it was assumed for the computations here that curvature occurred 
in the plane of the nozzle axis, the most unfavorable case. 
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TRANSIENT RESPONSE AND FREQUENCY CHARACTERISTIC OF 
DIFFERENTIAL PHASE-SENSITIVE AMPLIFIER-RECTIFIERS 
WITH RHEOSTAT-CAPACITIVE LOADS 


V. 1. Anisimoy 


(Leningrad) 


A method is considered for computing the transient response and frequency 
characteristics of differential phase-sensitive amplifier-rectifiers with rheostat- 
capacitive loads. An equivalent circuit to the amplifier-rectifier is given and 
the pulsation coefficient of the output voltage is computed, 


Phase-sensitive amplifier-rectifiers are widely used in various automatic regulation and control circuits 
for transforming ac voltages into dc voltages, the magnitude and polarity of which depend on the phase shift 
between the signal applied to the circuit's input and sorne standard voltage. Most frequently, phase-sensitive 
amplifier-rectifiers operate in a mode where the phase shift between the input and standard voltages is either 


zero or 180°. In this case, a variation of the phase shift of the input signal by 180° corresponds to a reversal of 
polarity of the output voltage. 


Besides an ac voltage, a dc voltage (usually a feedback voltage) is frequently applied to the input of the 
phase-sensitive amplifier-rectifier. With this, the amplitude and polarity of the output voltage are determined 
by the amplitude and polarity of the dc input signal, 


In order to decrease pulsation of the output voltage, the load of the phase-sensitive amplifier-rectifier is 
ordinarily shunted by a capacitor, which leads to an increase in the time required to establish the output voltage 
after a signal is applied to the circuit's input. If the phase-sensitive amplifier-rectifier has applied to it a 
modulating voltage, the amplitude of which varies sinusoidally, then the presence of a capacitor leads to a 
frequency distortion of the signal envelope. The steady-state behavior of phase-sensitive amplifier-rectifiers is 
considered in[1]. The present paper is devoted to the analysis of the transient responses in differential circuits 
for phase-sensitive amplifier-rectifiers with rheostat-capacitive loads, The purpose of the investigation is to 


obtain calculated relationships and equivalent circuits which will be useful for engineering designing of phase 
sensitive amplifier-rectifiers, 


1. Basic Equations and an Equivalent Circuit for Differential 


Amplifier-Rectifiers 


The circuit for a differential phase-sensitive amplifier- rectifier is shown in Fig. 1. The circuit given 
corresponds to the most usual case, when an ac signal voltage, wp, and a de signal voltage, ujn., are im- 
pressed upon the input of the amplifier- rectifier, 


In the circuit shown, the grid bias on the tube's grids is provided by connecting, in the grid circuit, 
sources of dc voltage, U,,, and ac voltage, Upj;, and also connecting them to resistor r,, for purposes of auto- 
matic biasing. Naturally, only one of these methods of biasing is used in any given actual circuit but, for 


generality of the exposition, we shall assume that the biasing is implemented by simultaneous usage of all three 
of the methods given, 


5 =. 
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The load on the phase-sensitive amplifier-rectifier can be either the resistances rg (divided load) or the 
resistor r; (non-divided load). If the rectifier operates on a divided load, the resistor ry, is, obviously, not pre- 
sent. 


Fig. 1 


Thus, the most common case is when the rectifier operates on a nondivided load, This case is also used 
as the basis of the analysis which follows, all the computed formulas being valid for divided loads if the sub- 
stitution = is made in them. 


For setting up the equations describing the circuit of the phase-sensitive amplifier-rectifier, we made the 
following assumptions. 


‘a 


in 


Fig. 2, Tube characteristics of the phase-sensitive amplifier- 
rectifier, 


_ 1, The tube characteristics are linear (Fig. 2) and, consequently, the equation of these characteristics, 
for negative voltage on the grids, has the form: 


"% 4 


2. At the circuit’s output is connected a capacitor of sufficiently large capacity that, in the steady-state, 
the voltage on the capacitor may be considered to be constant. 
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3. In the transient response, the amplitudes of the ac voltages and currents, and also the dc components 


of the rectified voltages and currents, vary insignificantly during one period of carrier frequency, so that these 
quantities may be considered as being slowly-varying functions of time. 


4. The phase shift between the ac voltages applied to the tubes is either zero or 180°. 


5, The tubes operate without grid current, and the average value of the anode currents during one period 
of the carrier frequency does not equal zero. 


If we bear the foregoing assumptions in mind and, for definiteness, if we assume that the signals up 
and yj. are applied in such fashion that the anode current in tube L, increases while the current in tube Ly 
decreases, we may write the following equations for the anode currents: 


where 

Ue, = Emecos wt — Uy — (iay baa) Tk» 

= Ema cos wt — U eg — + ing) Tee 

= Ugo + Umpi —2U min).c08 ot — U (iar + tas) 

= Ugo + (Um bi + WU min) 008 wt + + (ian + bas) ree 

Vins 
n= is the transformation coefficient of the input transformer and = 
By solving these equations for the currents igy and igs, we get 
ie = lEmae m in (1 + cos wt + (1) 


+ (1 + 24) — goe— (1 + 1) + 
iat = — min (1 + 27)} cos wt — 


— pin = (1 + 21) —pU goe—Uee (1 +1) + Vert, @ 
where 
4 
Emae = Ema— PU mpi» =Ugp 


The cut-off angles, 6, and 63, can be determined by equating the instantaneous values of the currents 
ig, and igs to zero, whence 


— (1 + 27) + (4 + 1) 


cos 8, = Ema et m (+ 21) 


BY gee Wig (1+ 2 + 


cos 8, = moe PU min (i + 2y) (4) 


The average value of the anode currents may ve determined by integrating Expressions (1) and (2). 


By taking into account that, in accordance with our assumptions, Urpin, Ujne Ucy and Ucg are 
slowly varying functions of time, we obtain 
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1 
le = \ dwt = ((Emae min (1 + sin 8, — 


— goe— = (1 + 24) + (1 1) — %), (6) 
+6, 
1 
lea = | = (mee (1 + sin 8, — 
— goet+ Win =(1 + 24) + -+ 1) — 9). (6) 


If the voltage pulsations on the capacitors are small then, for the most general case of a non-divided load, 
the phase-sensitive amplifier- rectifier may be described by the following equations: 


dU U U.,—U dU U 


By substituting the values of the currents, 1g, and Igy, in these equations, and by taking account of 
Expressions (3) and (4), we obtain, after some transformations, 


aU. U. 
27) [Ema -+ min (1 + 27)] (sin — 9, cos 6,), 
dU U U..—U 
24+ 


= |Emae—wnU min (4 + (sin 8, — 6, cos 64). (8) 


The simultaneous solution of Eqs. (7) and (8) determines the variation of the output voltage when a signal 
is connected to the input of the phase-sensitive amplifier-rectifier. However, the solution of these equations 
may only be obtained by a numerical or graphical integration for given values of the circuit parameters. 
Naturally, such a method is ineffectual for practical usage and it is therefore more efficacious to simplify some- 
what the expressions defining the mean value of anode currents so that one might obtain a solution which, even 
if somewhat less accurate, would be of a more general form. For this purpose, it is convenient to use the common 
piecewise linear approximation function for » = sin @ — @cos@, choosing as the independent variable cos 6 
[2, 3]. 


Ordinarily, for approximating the function », one uses either two straight lines, »; = A— Bcos 6, where 
A= 0,98, B = 1.22 for 0< cos6< 0.7, 
A=0.41, B=0.41 for 0.7 cos6< 1, 
or approximates it by one line, where 
A=0.8, B=0.8 for 0<cos6< 1. 
Approximating the function » by two lines leads to the necessity of connecting the solutions at the point 
cos @ =0.7, which significantly complicates the computations for the nonsteady-state processes, Therefore, 


it is more efficient to use one line to approximate the function ». The decrease in accuracy of the solution 
in this case is completely warranted by the simplicity of obtaining the computational relationships, 


If the approximation cited is used, and Expressions (1) and (2) are taken into account, the system of 
equations describing the circuit takes the following form: 


+6 

} 

rL 
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U, 


+ {[Emae + m in(1 + 27)) A — 
— [20 — WU (1 + 24) + Ver (1 +4) — BY, (9) 
i 
— min (1 + 27)] A— 
— [WU -+ in = (1 + 24) + Vea (1 1) — Vert) BY. (10) 


~ 


The system of equations just obtained allows, in case of necessity, the determination in general form of 
the temporal variation of the voltages Ucy and Uc3. 


Ordinarily, the output coordinate of the circuit is the voltage Upy;. Subtracting Eq. (10) from Eq. (9), 
and remembering that | 


Vout = Ue, — 


we get 
where 


r 2r 


If we use the single-line approximation to the function y, the expression for the coefficients a, b and 
Te will have the form: 


0,84 


Equation (11) allows one to obtain, in general form, the law of temporal variation of the output voltage 
for any form of signal envelope applied to the input of the phase-sensitive amplifier- rectifier, 


It should be borne in mind, however, that this equation is only valid if the condition holds that the mean 
value of the anode currents is not zero, i.e., if the inequality holds 


(U + 2U min) (1 + 21) < Emae — PU gee — (1 + 1) + Vert: (12) 


If Inequality (12) holds, it is always the case that cos @, < 1 and cos @, < 1 and, consequently, the mean 
value of anode current during one period is not equal to zero, 


If Inequality (12) does not hold, one of the tubes is cut off during several perio is of carrier frequency and 
the capacitor discharges through resistor rq until such time as the tube again conducts, i.e., Inequality (12) again 
becomes valid, With this, the variation of the voltage on the capacitor is determined by the equation 


aU, 


Since ordinarily T > T,, the time to establish the output voltage is significantly increased in this case, 
If, at the circuit input, there is impressed a modulating voltage, the amplitude of which varies sinusoidally, the 
inapplicability of Inequality (12) leads to the appearance of a strongly nonlinear distortion of the output voltage, 
since the form of the voltage at the circuit's output is, in this case, nonsinusoidal, 


e- 
mon a 
| 


Inequality (12) does not allow a direct determination of the magnitude of the input signal for which the 
average value of anode current does not equal zero, since this inequality contains the values of Uc, and Ucg 
which are determined by solving the system of Eqs. (9) and (10), 


By solving Eqs. (9) and (10) for a definite form of input signal envelope, 


rt and substituting the solution in Inequality (12), one can obtain an expression 
for immediately estimating the voltage of the signal for which the average 
| er u value of anode current does not reduce to zero, 


g Thus, for example, for the most difficult case, when there are si- 
multaneously impressed signals of an ac voltage Um in, of invariant ampli- 
Fig. 3 tude, and a dc voltage Uj,., Inequality (12) takes the form: 


Ema e (1 — — goe (1 — 


Uin= + in< (13) 


Aa 


Ba 
1 #(i+27)+ Ba’ x(i+2y)+ Ba 


If the inequalities given above hold, one can construct an equivalent circuit to the differential phase- 
sensitive amplifier- rectifier, 


To this end, we write Eq, (11) in the following form: 


dU 
qT. U out bEe, (14) 


where 
= min + 2pUin 


i 
lie "le 

Be rae "le 
+ Ba rae" 
Tact "Le 


2nr 
» "Le = TL, Ce = 0.5C. 


In correspondence with Eq, (14), the equivalent circuit of the phase-sensitive amplifier-rectifier has the 
form shown in Fig. 3. It has already been stated that the best choice for the approximation coefficients A and B 
is 0.8. 


For this case, 


2p + nU min), 
2.5nri, Tee = = "Ls Ce=0.5C. 


By using Eqs, (11), (14), or the equivalent circuit, one can determine the complex transmission coefficient 
of the phase-sensitive amplifier-rectifier, 
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For an ac signal voltage, the complex transmission coefficient is 
Tout 
where 
4 PP 
2un 
Ta 
i- "ie + ra r Le 
is the circuit's static gain for an ac signal voltage, 
The complex transmission coefficient for a dc voltage signal is 
) 
where 
Vin. Tae? 
"Le 
is the circuit's static gain for a de voltage signal, 
) 2. Pulsation Coefficient of the Output Voltage 
The equivalent circuit and the basic equations of the phase-sensitive amplifier-rectifier do not allow de- 
finitive recommendations to be made respecting the choice of parameters and the mode of operation of the 
amplifier- rectifier, since in designing such a circuit one is ordinarily given the pulsation coefficient of the out- 
put voltage, Therefore, for a more complete estimate of the qualitative indicators of the phase-sensitive ampli- 
fier-rectifier, it is necessary to find the relationship connecting the circuit's time constant with the pulsation co- 
efficient. 
The pulsation coefficient of a rectified voltage is usually defined as the ratio of the amplitude of the 
fundamental harmonic to the dc component of the voltage. 
Since the impedance of condenser C at the carrier frequency w is much less that the impedancest, and 
t,, the amplitude of the output voltage's fundamental harmonic can be determined approximately by the ex- 
pression 
The amplitudes of the first harmonics of the anode currents, I}, and If, will be 
i | 
nt With a sufficient degree of accuracy, the function g =@~— sin @ cos 6 can be approximated by the 


straight line gy; =M~—N cos @ with the values of the coefficients being M = N = 1,57 [3]. 


{ 
{ 
| 
j 
i 
165 


By using this approximation, and bearing in mind Expressions (11) and (12), we obtain 


= Ema e+ unm + 21) — + Win = (1 + 24) — 
— Ue (1 +- 4) + Vesti, 
= — min (1 + 24) — goe— 
—pU (1 + 24) (1 +1) + 


By substituting the values of I}, and Ip, in Eq. (15) and by remembering that in the steady state 


Uo our = m in + Uin-) (for A = B), 
we get | 


Um out = Seats, Vo out 


Thus, the pulsation coefficient for the output voltage is 


i—a 1,96 2 
where 
2a 
By recalling that 


Te . 


"ae ‘Le 


one can express the circuit's time constant, T,, in terms of the pulsation coefficient of the output voltage: 


(for A= B=0,8). (17) 


140.25 755 


For operation of the rectifier on a divided load (tr, = a), we have 


Te = (18) 


The expressions obtained allow one to draw certain conclusions regarding the choice of parameters for the 
circuit of a differential phase-sensitive amplifier- rectifier. 


As follows from Expression (17), the circuit's time constant can be lowered due to an increase in the | 
carrier frequency w, of the pulsation coefficient p, of the coefficient of load division a, and by means of a 
decrease in the undivided load coefficient 8, The first two methods of decreasing the circuit's time constant 
are inapplicable in practice, since ordinarily the carrier frequency and pulsation coefficient are given in the 
design. By varying the coefficients « and 6, i.e., by selecting the necessary relationships between the tube 
impedances and the divided and undivided loads, one can vary the circuit's time constant within broad limits. 


Thus, for example, for rectifier operation on a divided load (6 =0) and for p=0.1 and f =50 cycles, we 
have 


- 
i 
1 
| 
fi 
| 
| 


Te =0.05 seconds fora=1, TT, =0.01 seconds for a = 20, 


Thus, the circuit's time constant was decreased by a factor of 5 by varying a from 1 to 20, 


Consequently, to decrease the time for establishing the output voltage, and decreasing the frequency dis- 
tortions for planned differential phase-sensitive rectifiers with rheostat-capacitive loads, it is necessary to select 
tubes with low internal impedances, r;, but resistances of divided, rg, and undivided, r; , loads should be chosen 


as high as possible. A given pulsation coefficient is thereby provided if the capacitance C is properly chosén, its 
magnitude being determined starting from Formula (16). 


SUMMARY 


1. The approximate analysis of the circuit for a differential phase-sensitive rectifier with rheostat- 
capacitive load allowed the construction of an equivalent rectifier circuit, described by a first-order linear 
differential equation, If the input signal does not exceed a definite magnitude, then the phase-sensitive amplifier- 


rectifier is an aperiodic first-order link with respect to the input signal envelope. With this, the time constant 
and static gain of the link are easily determined from the circuit parameters. 


2. The pulsation coefficient of the output voltage depends on the magnitude of the output capacitance 


and on the impedance of the divided and undivided loads but, to a first approximation, does not depend on the 
tubes’ internal impedance, 


3. The time constant of the phase-sensitive amplifier-rectifier with rheostat-capacitive load decreases 
with an increase in carrier frequency, pulsation coefficient, and ratio of the impedance of the divided load to 
the tubes internal impedance, The inclusion of an undivided load impedance, r,, increases the circuit's time 


constant (with the condition that the pulsation coefficient of the output voltage remains invariant due to the 
corresponding variation of the capacitance C). 


4. Experimental verification of the formula herein obtained showed that the error in determining the 
static gain by the formulas given above did not exceed 10% to 20%, and the error in determining the time 
for establishing the output voltage varied from 20% to 30%. 
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GRAPHICO-ANALYTIC METHOD FOR COMPUTING THROTTLE CONTROL 


CHARACTERISTICS OF ASYNCHRONOUS MOTORS WITH 
LARGE SCHENFER ROTORS 


O. B. Rosenbauli and R. N, Rodin 


(Moscow) 


A graphico-analytic method is presented for the calculation of throttle 
control characteristics of asynchronous motors with large Schenfer rotors, and 
an example is given of the construction of such characteristics. 


In highly automated electric drive systems containing asynchronous motors, one of the most sensible 
methods for speed control is the method of control via choke (throttle) saturation, Motors with large steel 
Schenfer rotors are frequently employed in such systems. This is explained by the fact that such motors, when 
used in systems of speed control, possess a number of significant advantages over motors with rotors of the 
“squirrel-cage” type [1]. 


Despite the widespread use of systems with throttle-controlled speed of asynchronous motors with Schenfer 
rotors, methods of computing the characteristics of such drives have, until now, not received their due attention 
in the literature. The present paper presents a graphico-analytic method of constructing the characteristics of 
throttle-controlled asynchronous motors with large steel rotors. 


This method is useful for various forms of drive schemes, using various methods of compounding by current 
and voltage. 


. Computing the Mechanical Characteristics of a Drive 


In computing the mechanical characteristics of a drive, the initial data are the throttle saturation parameter, 
the feedback coefficients of current and voltage, and the motor's natural mechanical characteristics. 


It is assumed that in all modes of drive operation, the currents and voltages in the chokes ac circuit are 
close to a sinusoidal form, Such an assumption is justified by the good agreement of experimental results with 
computed data, and also by the possibility it offers of using methods with complex variables, which significantly 
simplify the computations, It is also assumed that the leakage of magnetic flux and losses in the throttle's 
magnetodrive are negligibly small. To decrease the throttle's weight, it is desirable to use, for core preparation, 
high-quality magnetic materials with large values of magnetic permeability and inductance, such as types 50NP, 
65NP, etc. The active loss in such cores is extremely small, Moreover, the presence of active losses affects the 
position of the simultaneous magnetization curves, and the losses themselves are taken into account in the com- 
putations. Therefore, the assumption of negligibly small losses in the throttle cores is completely justifiable. 


For illustration, Fig. 1 shows a widely known electrical circuit for one form of throttled drive with current 
and voltage compounding. 


The computational method presented below is illustrated in terms of a drive using a power throttle. How- 
ever, the method can be extended for use with drives where the throttle is realized in other forms. 
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We shall assume that the throttle for each phase consists of two three-leg cores with ac windings connected 
in parallel, Then, the equation for the magnetizing force for the choke's ac circuit will have the form: 


Ww. 
Aan = =, 
0.4nV 2 


where H,,, is the field reversal force in the throttle's core, 
in oersteds, W,, is the number of turns in the ac winding 


~U per throttle core, 7,, is the average length of the mag- 
‘ — A 98 netic lines of force and 1_ is the drive's phase current. 
Cc The control and feedback windings are common to 


all chokes, With this arrangement of the control windings, 
the magnetizing force for the control circuit may be 
. CT described in the form 


Hylgy=0.4e (IW, 


Equation (2) shows that, even for 1, =0, the action 
of the feedback results inH, # 0. 


rv. 
¥ 


M We take the direction of the current vector to be 
along the real axis, Using the vector diagram of the 
Fig. 1, Electrical circuit of the throttled drive drive's voltage and current, shown in Fig. 2, we write 
with current and voltage compounding. M) the equation of the ac circuit in complex form 
The motors CT) the current transformer; : 
VT) the voltage transformer; SC) the satu- Uy = (U~ cos +] 
rated chokes; D, and D,) rectifiers (diodes) (3) 
and Ry and R,,) adjusting resistors. sin Px): 
Here, U,, is the network voltage, U_ is the phase 
¥ U. voltage on the motor, R_, is the active impedance of the 
choke's ac winding, is the choke's emf and 
Vi, sing is the angular shift between the motor's voltage and 
current. 
In terms of moduli, we have 
U. m 
U} = (U~ 008 + w 
Fig. 2 + (E Unsin 


The ampere-turns of the feedback can be determined from the formula 
Ip, Wa, = kU Wy— k, Wil _, (5) 


where k and ky are the feedback coefficients of, respectively, the current and voltage, The values of the co- 
efficients ky and kg are determined in computing the parameters of the feedback circuit (see [1]) and will there- 
fore be considered as initial data in calculating the drive characteristics. 


Simultaneous solution of Eqs, (1), (2), and (5) gives the following expression for the voltage on the motor: 


4 lay 
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The choke's emf in volts can be found from the expression 


E = 4.44{W By S-10-* (1?) 
or 
E y= (8) 
where f is the network frequency, B,,, is the amplitude of inductance in the choke's cores and S is the cross 
f section at the center of the magnetic circuit bar. 
By solving Eqs. (1), (4), (6), and (8), simultaneously, we obtain the equation of an ellipse 
| 
| 
Here, 
LW H 1.W, H 
Bs = Bm + helm + hy = By + — (10) 
LW H 
(12) 
U 
Hy = (13) 
Ww. 
sin 
” a@ 
sin 
k, = (15) 
0.2nV2W _ cos m 
Ke = — (6 
u 


Expression (9) defines an ellipse with semi-axes By and H, on the By, H, plane. On the Bi,, Hyp plane, 
the same expression corresponds to an ellipse which, as control current I’, varies, moves along the straight line 
which makes the angle a with the H,, axis. With this, the ellipse’s semi-axes, B, and Hy, remain parallel to 
the axes B,, andH,,. To obtain the geometric locus, on the Bm, H,, plane, of the points defined by Eq. (9), 
it is necessary to subtract the ordinates of the line B= kgH,, from the ordinates of the ellipse constructed on the 
Bin» Hyp axes. 

For a motor with aSchenfer rotor, working with a controlled drive over a broad load range, the angle 
%m_ Varies insignificantly [1]. Therefore, the method for constructing the mechanical characteristics of the 
drive can be reduced to the following. 


1. We prepare the master e!lipse with semi-axes By and Hy, 
2. We find the subsidiary functions: 


ke 
Hm= He an 
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and 


Hm = (18) 
av 


3. On the B,H plane we construct the line at the angle a, whose tangent, in the scales of B and H, equals 
the ratio kg 


4, For a number of values of H,, for which there are curves of simultaneous magnetization, we determine, 
by the use of the function in (17), the abscissas of the centers of the ellipses with semi-axes By and Hy and super- 
impose the master (template) at these points, as shown in Fig, 3, a. With this, the semi-axes of the template 
must be parallel to the axes B,, and Hp. 


a b 
6,, | Template 6m| Template “0 
Yi, 
WV 


Fig. 3 


5. On the template we find the points of intersection of the ellipse with the corresponding curves of si- 
multaneous magnetization, 


6. By subtracting the ordinates of the line By, = kppH,, from the ordinates of the curve constructed from 
these points, we obtain the curve which is the geometric locus of the working points of the choked drive in the 
Br» Hyp plane when the choke's control current equals zero, 


7. To determine the geometric locus of the working 
Sn . points when |, # 0, it is necessary to shift the center of the 
25.U 2 ellipse, corresponding to a definite value of H,, to the left, 
ai as shown in Fig, 3, b, the magnitude of the shift being de- 
fined by Expression (18), after which one proceeds with the 
\ construction as described in numbers 5 and 6 above. 


8. After having thus obtained, on the B,,, Hy, plane, 
the family of working point curves for various values of I,, 
we construct the subsidiary lines in the fourth quadrant 


0.28 2W_ Hm = (Hm), 


and, in the third quadrant, 


LR. 


Fig. 4 U_=I_R_=f (I). 


In the second quadrant we draw an arc of the circle, centered at the origin of coordinates, with radius 
equal to the network voltage Up. Being given a series of values of Hm, we then, by means of the supplementary 
construction shown in Fig. 4, find the geometric locus of the points 0 which are the end-points of the vectors 
Ech Corresponding to the series of values of control current, Here, the scale of Ech is determined by Expression 
(7). Then, from these points, at the angle gp, we draw straight lines up to their intersections with the circle 
Un. As is easily seen on the vector diagram in Fig. 2, the segments OA obtained as a result of this latter 
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construction have the scale of the voltage on the motor, Um. Using the values of U,, thus obtained, we con- 
Mi struct the relationship, in the third quadrant, 


Um (I~). 


9. Having the calculated or experimentally determined family of the motor's mechanical characteristics, 
n= f(M), for a number of values of current and voltage, and using the relationship, U,, = f (I), obtained in 
number 8 above, for various values of I,, we construct the drive's mechanical characteristics for various control 
currents or ampere- turns, 


By, gauss 
47] 
160 120 80 2 
cae vai Oprsted 
| a4 


Um 1. 
Fig. 5 

For calculations of drives with other forms of power chokes than the type taken here, Expressions (13) to 
(16) will differ in dependence on the form of Eqs. (1) and (2), written for the specific form of choke. In this 


case, in analogy with the deductions from Eqs. (13) - (16) used in this paper, one should find new expressions 
for H, and the constant coefficients and then construct the drive's mechanical characteristics by following the 


sequence of steps given here. 
We now illustrate the method of constructing a drive's characteristics by determining the relationships 


Um = fi.) and n= f (M) for various control currents for the drive with the following parameters: 
U = 220y, 0,82, P=0.2 kw, W_=1700 turns, W,, = 158 turns, 

W,=265 tums, W,=33tums, 4, = 33cm, R_ = 11.5 ohm, 

| f = 50cps, = 2.1-10-* 1 = 3.8, cos 0.58, 

) the choke material being £4 A steel, 

1. We determine the values of the coefficients: 


laysin V2W_ sin? m 


0.34-108,. 


ky W 
u 


7112 


- 


2. We find the ellipse’s semi-axes: 


U,0.2nV2W 
== 54 oersted, 
ky Wi 


tay (Re + 


3. We find the expressions for subsidiary Relationships (17) and (18): 


U, 
B= 8.34-10° gauss, 1, = 


ke 2 keW, 
Hm, = Hy, = I, =1A4lg. 


4. In correspondence with steps 3-9, we construct the relationship U,, = f(.) (ig. 5). 


5. Using the experimentally determined mechanical characteristics of the motor, shown in Fig. 6, and 
also using the relationship U,,= f (l_), obtained in step 8, we determine the characteristics for various 
values of control ampere-turns, The constructed mechanical characteristics are given in Fig. 6. 


Um = 
66v 58 "0 132 154 


1, 0528 a4) as7 a656 082 


Fig. 6. Motor's mechanical characteristics, Along the axis of abscissas 
is plotted the relative moment on the motor shaft, while the motors 
relative speed is plotted along the axis of ordinates, 


The example shows that the method presented allows a relatively simple computation of the mechanical 
characteristics of a throttle-controlled asynchronous motor with a large Schenfer rotor. 
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THE USE OF HALL-ELEMENTS AS PHASE-SENSITIVE RECTIFIERS 


Vv. N. Bogomolov and V. A, Myasnikov 


(Leningrad) 


The operation of a phase-sensitive detector, based on the Hall effect, is 
considered, and certain basic relationships are provided. A simple graphical 
method is presented for choosing elements to compensate for the temperature 
dependence of the Hall-element parameters, 


A table of experimental data for three different phase-sensitive detectors 
is given. 


Fundamental Relationships 


As is well known, the Hall effect in semi-conductors is characterized by the magnitudeof the increasing 
Hall emf 


Uy = RZ 10-, | 


where R is the Hall constant in cm*/coulomb, | is the unit's current in amperes, H is the magnetic field force 
in oersteds, d is the thickness in centimeters of the semi-conductor disk and Uy is the Hall voltage in volts. 


Frequently, the coefficient of voltage transmission is used, this being the ratio of the Hall emf, Up, to the 
voltage U, impressed on the disk and generating the current I: 


(2) 


By using the formula for the electrical conductivity of the semiconducting disk 


Il 
° = iad (3) 
and substituting (1) in (2), we obtain 
K 


where y » Ro is the current carrier mobility in the semiconductor in cm” /volt-second, and a and 1 are, 
respectively, the width and length of the unit in centimeters, In the formulas for the Hall emf and the voltage 
transmission coefficient, it is assumed that the current flowing through the unit, and the magnetic flux through 
the unit, constitute sources of constant emf and mmf. 


For a phase-sensitive ac Hall rectifier, the voltage at the unit's output contains two components: a 
constant (dc) output and an ac output which changes at twice the supply frequency. In this case, two transforma- 
tion coefficients should be introduced: a transformation coefficient for the dc voltage and a transformation co- 
efficient for the voltage at the doubled frequency. 
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The transformation coefficient for the dc voltage is 


= (5) 
where K is the transmission coefficient and Uj; is the dc component of the Hall emf. 
In Formula (5), the peak values are chosen for the magnetic field and the applied voltage. 
The transformation coefficient for the voltage at doubled frequency is 


U, 


where U;; _ is the component at doubled frequency of the Hall emf, 


In Formula (6), the magnetic field can be taken either at peak value or at virtual value. 

Single-crystal semiconducting units, in the majority of cases, are prepared in the form of parallelepipeds 
with the ratio —+ =1.5-2, In its construction, a phase-sensitive rectifier which uses a Hall-unit is a choke in 
whose air-gap a semiconducting lamina, a Hall-unit, has been placed, 


Phase-Sensitive Rectification 


Phase detectors, based on the Hall effect, are effectively used in those cases where high power is not 
required, for example, for operation with the control winding of a magnetic amplifier, with a relay winding, etc. 
Such detectors are linear, reliable and long-lived. 


Let us apply a voltage uy =U, sin (wt + ¢ ) to the 


J —o unit's current electrodes, atid a voltage ug= U, sin wt to 
Ty the choke winding, where is the phase-angle shift between 
— a U, and U3, The magnetizing current in the choke will equal 
Yin Dy 
o— 
Tr, 4 is = sin (or 
Ps and the magnetic field force in the choke's air-gap is de- 
termined by the expression, H = where K, is the pro- 
& portionality coefficient. The current through the model 


Fig. 1. Circuit for compensating for the equals i; =Ip, sin (wt+ g ), By using Formula (1) we 


second harmonic component at the Hall 
output. 


Uy = [sin 9 + sin (2ot + (1) 
where 
K. R K,-10-* 


It follows from the expression obtained that, at the output of the Hall detector, there is a dc component, 
varying with the sine of the phase-angle shift, and an ac component with doubled frequency. It is easy to 
convince oneself, that, at the Hall electrodes, a dc component will be produced only in case that field frequency 
and current frequency through the mode] are absolute identical, Since, it is assumed that the phase rectifier 
operates on a fixed frequency then, by shifting the current in the choke or the unit by 90°, one can obtain a 
dependency of the dc component on the cosine of the phase-angle shift between the voltages on the choke and 
on the Hall unit. 


115 


4 
K. = =;K, (6) 
| 
4) 
h | 
rma- 
| 
| 


If the presence of the second harmonic of the carrier frequency at the detector output is undesirable, one 
may use a circuit which allows compensation for the ac component at the detector output (Fig. 1). With this, 
the power supplied to the load by the dc component is doubled. The transformers, Tr 1 and Tr2, in the circuit of 
Fig. 1 are matched. The capacitor, C, is used for tuning and resonance, We assume that through D, there 
passes a current i, =1, sin («'t+y). The current inducing the magnetic field is ij =Ij cos wt. 


At the Hall output of unit D, there will be a voltage 
i 
Ou, = 1, K, [sin (20t + 9) + sing). 
Analogously, one may write for the second unit, D,: 
ig =I, cos(wt+ 9), i,=J;sin ot, 


Uy .= + K, [sin (2t + — sin 9], 


where Ky and K, are proportionality coefficients, 


TABLE 
Permall 
Units parameters unit °y Ferrite unit KhVP unit 
Material of specimen InSb InSb Ge 
Type of magnetic circuit Sh 12 Sh 12 Sh 42 
Bat cross-sectional area, 
cm 0.72 0.72 0.72 
Air-gap size, mm 0.4 0.4 0.4 
No, of turns of the choke's 1600 192 1600 
core PEL 0.22 PEL 0.54 PEL 0.22 
Semi-conductor dimensions, mm 8x4x0.3 6x4x0.3 8x4x0.3 
Unit's output impedance 0.62 (B=0) | 1.34(B=0). 1000* 
Ty, ohms 0.83 1,44 
B = 13 280 ‘¢ = 3360 
Unit's input impedance, 0.4(B =0), 47 (B = 0), 800 
ohms 1.55 | 6 
(B = 13 280 gs)| (B = 3360 gs) 
Maximum current throngh unit, amp 1,2 0.4 0.025 
Maximum impressed voltage, v 0.6 0.64 20 
(B = 13 280 gs)| (B = 3360 gs) 
Freq. of current and mag. field, cps. 2250 50 
4 mw 1.4 2 
Maximum output power, mw (B= ey (B= meen (B = 15000 
a a a 
Max, transformation 2 
coefficient vs 14 13 10 
Total power required by the 
unit, w 2 i 2 
Wt. of the unit with matched trans- 
. formers 250 220 250 
Current camtler mobility 7000 22000 2500 


* Depending on the germanium frequency, the output impedance can vary from a hundred 
to several thousand ohms. 

** By using InSb with a current carrier mobility u = 8000 to 100,000, it is possible to 
obtain K = 30 % with a corresponding power of about 20 milliwatts for the same dimen- 
sions of the unit. 


If we connect the outputs of units D, and D, in opposition, and set 1, = 15, Ij = 1; , Ky = Ky, we obtain, 
finally, Un, Up, = K; sin y. With a matched load, ry, we obtain a power output twice that which 
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would be obtained from one unit. In particular, such a circuit can be successfully used as a power relay in relay 
protection systems. The reference voltage, usually constant in magnitude, is fed to the choke core and the signal 
to be detected is applied through the model. Such connections allow the circuit to have a small time constant 
since, in practice, the semiconducting unit is inertialess. For the Hall effect, the time constant lies between 
10°" and 10° ™ seconds. 


The essential difference between linear phase-sensitive detectors based on the Hall effect and linear de- 
tectors using vacuum tubes, crystal diodes or crystal triodes is the following. While the output of a Hall effect 
detector is a voltage with only sum and difference frequencies, the output of any other type of detector contains 
a large number of voltages of different frequencies, part of which may be compensated for by nonlinear elements 
or by other methods, This leads to an instability in time of the detector's null as the temperature varies, In 
theory, the Hall-effect detector is free of this disadvantage. 


ma 


Fig. 2 Fig. 3 


We now consider the experimental! results of investigations in the operation of rectifiers built with units of 
indium stibide and germanium (InSb, Ge). The investigation was carried out at a frequency of 50 cycles with 
choke cores of Permalloy and KhVP-steel] and at a frequency of 2250 cycles with choke coils of ferrites. The 
parameters of these devices are given in the table, 


In correspondence with Formula (7), Fig. 2 gives the dependence of load current on the phase-angle shift 
between the current flowing through the unit and the current in the choke, for the ferrite unit with r; = 5 ohms, 
1; = 0.2 amperes and B = 2500 gauss. 


a) Input and output unit impedances, A Hall-effect unit may be considered as a generator whose internal 
impedance depends on the magnetic field strength. The effect of varying the impedance of a semiconductor 


in a magnetic field is estimated by the formula, a w A(wH)*, where ry is the semiconductor’s impedance 


for magnetic field swength H=0, Ar is the increment of impedance in a magnetic field of strength H and A is 
a coefficient embodying the semiconductor's properties. 


The graphs of Fig. 3 show the dependence of the semiconductor's impedance (its constant component) on 
the magnitude of the magnetic field, for a Permalloy unit. Here, rj is the impedance between the current 
electrodes and ny, is the impedance between the Hall electrodes. It follows from the curves for my that the 


impedance of the Hall output varies very significantly, increasing by 34% from the magnitude at which the 
induction B = 0, 


b) Power of the constant component of the unit output voltage. The mode of operation which gives 


maximum power in the load, as is well known, occurs when the load and generator impedances are equal, an 
assertion which is also correct in the case of Hall-elements, The power delivered to the load is determined by 
the expression 
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By letting yy =r, one can write P; = Uj,/4t,,. The power delivered to the model is Pj = U?_/2r. If 
Formulas (1) and (3) are used, the power transmission coefficient may be written as 


2 


Consequently, to increase the power transmission coefficient, one should choose semiconductors with high 
current carrier mobility, and one should increase the inductance. It is also possible to increase the power drawn 
from the unit by increasing the current through it; with this, the upper limit on current is determined by heat 
transfer conditions. In steady-state operation, the heat generated in the unit due to passage of current will equal 
the heat dissipated by the unit in unit time: 


_ Yin dao 


The quantity of heat dissipated from unit surface of the mode! in unit time is: 
= w/cm* +sec 


Without special measures having been taken for dissipating heat from the unit, i.e., for operation directly 
in air, q< 1. The magnitude of q can be increased by taking special measures to improve the heat transfer 
properties, for example, by immersing the unit in a tank with oil or glycerine. Thus, for “standard” units whose 
dimensions are no greater than 12x 6 x 0.3 mm, with no particular measures for heat transfer having been 
taken, the power drawn does not exceed 3 to 5 milliwatts at audio frequencies, 


6. gauss gauss @ 
Fig. 4 


c) Dependence of output voltage and current on the magnetic flux through the unit, Figures 4, a, b and 
© show, respectively, the dc component of the Hall voltage, the power transmission coefficient and the trans- 
formation coefficient for the dc voltage as functions of the magnetic field strength for the Permalloy unit. 
Beginning with a value of inductance approximately equal to 8000 gauss, the dependence of output voltage on 
inductance ceases to be linear. The deviation from linearity is the more noticeable the closer the load im- 
pedance is to the unit's impedance. 


Under conditions of signal detection, the current through the unit varies and, consequently, it is important 
to have a linear relationship between the current, ly, and the voltage, Uj, of the Hall-element output and the 
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variable current, lj, through the element, As Fig. 5 clearly shows, a Hall detector possesses linear character- 
istics until that value of element current is reached for which the magnitude of q becomes higher than the 
admissible value for the given type of unit. The dashed lines on Fig. 5 correspond to the maximum values of q. 


ma a Uy, mv b 


Fig. 5, Output current and voltage as functions of the current 
through the unit. a) Fora ferrite elements b) is for a Permalloy 
unit. 


If, through an InSb ferrite unit, there is passed a current I; = 0.8 amperes, twice as large as for q,,,,, then the 
current in the load will drop by 11% in 20 minutes as a consequence of the overheating, 


Dependence of the Parameters and Characteristics on 


Temperature and Thermocompensation 


One of the disadvantages of Hall-elements is the dependence, already noted, of their parameters, and 
characteristics on the temperature of the ambient medium. Thus, the transformation coefficient for dc voltage 
in the experimental units fell by approximately 20% at temperatures of 0° and 70° C from their value at 20° C. 
There are two possible ways to increase temperature stability. The first way consists in using semiconductors 
whose current carrier mobility and Hall constant have a small temperature dependency, For example, it is 
advantageous to use indium arsenite, InAs, [4] or HgSe [3]. 


The second way to increase temperature stability is by artificial temperature stabilization. Here two cases 
must be considered: when the input quantities are voltages or currents. 


Let us consider the first case, when the input quantity is voltage. 


a) Temperature compensation by parallel resistances. Let us turn to Fig, 6, where E is the voltage source, 
R is the impedance of the voltage source, p is the element impedance between the current electrodes and r 
is the compensating impedance, R and E being temperature-independent. In order that the Hall voltage does not 
depend on temperature, we should, in the expression 


(8) 
set 
'K 
(9) 
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where K‘ is the transmission coefficient. The new transmission coefficient, K' = BK, does not depend on tempera- 
ture, Since B < 1, the transmission coefficient of the device is lowered. It follows from Eq. (9) that 


= i (10) 


From this we find the condition for the existence of the compensating resistance r, namely, R < o(2 -1) 


In Expression (8), K and p depend on temperature, Knowing the functions K= f (t*) and p = f (t*), one can, 
by using Eq. (10), construct the function r = f(t*), If a material is selected whose dependence on temperature is 
analogous to r= f(t*) then, in the temperature range considered, the Hall voltage U,; will not vary. 


AL. 


Fig. 6. Circuit for tem- Fig. 7. Circuit for tem- Fig. 8. Circuit for tem- 


perature compensation by perature compensation by perature compensation 
parallel impedances. series impedances. when the input quantity 


is current. 


b) Temperature compensation by series resistances (Fig. 7). If we set the source's internal impedance, 
R=0, then the Hall voltage is 


KEp KE 


+! 


KU, = 


In order that Uj does not depend on temperature we should set 


= BK, 
r 
>t! 


where B < 1. In this case, the desired expression for the compensating impedance will be 


If R% 0, then the following expressions will be valid: 


r=£—R—p;, R<p(7—1). 
P 


When the compensating impedances are connected in series, one can use for them thermistors or impedances 
made of semiconducting materials (Ge, HgTe or InSb), i.e., materials with negative temperature coefficients. 


We now consider the second case, when the input quantity is current (Fig. 8). In this case, it is only 
meaningful to consider thermocompensation by compensating impedances in parallel. The Hall voltage is 
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Uy= KU, =K (11) 


The quantities which do not depend on temperature are the current I, the emf source E and the source 
impedance R: 


E 


In order that the Hall voltage does not depend on temperature, within a given range, it is necessary that, 
in Eq, (11), the current coefficient be made equal to the intermediate parameter 


pr 
R Ke 


this parameter having the dimensionality of impedance and not depending on temperature, Setting R' = KR", we 
define the compensating impedance 


the condition of existence of which isp > R". This determines the magnitude of R*. Knowing the functions 4 
K= f(t) and p = f (t°), one can find the function r= f (t*), In this case, one must use material with 
positive temperature coefficients as the compensating impedances. 


SUMMARY 


1, Phase-sensitive rectifiers made of Hall elements can be used for controlling the operation of magnetic 
amplifiers, relays, as output stages for dc amplifiers with transformers, etc, 


2. The dynamic range of detectors built of Hall elements is greater, on the weak signal side, than that 
of any other detectors. With this, the characteristics for the output current (voltage) are linear, starting from 
zero, 


3. In contradistinction to all other detector types, whose outputs contain a large number of voltages of 
different frequencies, the Hall detectors’ outputs contain only sum and difference frequencies, which provides 
a stable null, 


4. Phase detectors using the Hall effect have virtually infinite life spans, since the volume properties 
of semiconductors are used, 


5. The disadvantages of such phase detectors are low efficiency and low temperature stability. However, 
the latter may be increased by special choice of materials or by the use of temperature compensation. 


6. A circuit element with characteristics similar to those described above, in addition to its use as a phase- 
sensitive detector, can be used for solving a large number of practical problems. 
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LIMITING DIMENSIONS OF ELECTRIC ELEMENT WINDINGS 


S. P. Kolosov 


(Moscow) 


The possibility is considered of decreasing the dimensions of electric 
element windings at the expense of increasing their overheating. It is 
established that, due to the effect of the resistance temperature coefficient, 
these dimensions can be decreased only down to a definite minimum. 


In connection with the wide usage of temperature insulators and the even greater potentialities in this 
domain, the question arises of determining the effective limits of increased working temperatures of the windings 
of various electrical elements. 


The designer of electrical devices which are to operate at high temperatures must give particular attention 
to the technological problems, taking account of the significant changes in the magnetic properties of ferro- 
magnets when they are overheated, etc. These questions may turn out to be decisive in the choice of the upper 
limits of operating temperatures, but the possibility of decreasing winding dimensions, even with an increase in 
their overheating, is very tempting, 


Let us consider the relationship between winding temperatures and dimensions, 


For this, we shall start with the assumptions that the necessary useful life of the insulators can be guaranteed 
for any choice of operating temperatures, and that these latter are chosen only with respect to an advantageous 
decrease in winding dimensions, Today, such a posing of the problem is realistic only in certain domains of 
technology, where the equipment has an abridged life (for example, in aviation, particularly in pilotless aircraft). 
However, when the rapid development of heat-resistant materials is taken into account, it may be expected that 
in the near future, the problem may be thus posed for devices designed for general industrial usage. 


So that the discussion may have a concrete character, we shall consider initially the connection between 
winding dimensions and overheating for a dc electromagnetic relay. 


When the electromagnetic relay operates, the electromagnet develops a certain force F in the maximum 
gap equal to the opposing force of this gap. This force can be found from Maxwell's formula 


F = 
where @ is the magnetic flux in the air-gap and S, is the area of the effective cross section of the air- gap. 
With this, 
wo 
(2) 


where w is the number of.turns in the winding, R is the impedance of the magnetic circuit, r is the winding 
impedance for a given temperature of the ambient medium, Samp and Uc, is the voltage for operating the 
telay under these same conditions, 


4 


The number of turns of a cyclindrical winding is related to its dimensions by the following relationship (1): 


ray" (3) 


where p is the specific impedance of the conductor material (copper) at the temperature 9am, lay is the 
average length of a turn, I, and he are the winding dimensions (length and height of the core) and f, is the 
core's charging coefficient, 


Finally, based on the heat-balance equation, in a simplified form, one may establish the relationship 
between the winding impedance and its steady-state (operating) overheating: 


r(i+ a8) = 1530, (4) 


Umax 


where Uy is the nominal voltage of the network, K, = is the coefficient of possible increase in the 


network voltage in use, a is the resistance temperature coefficient of the winding material (copper) at the 
temperature S,pp. u is the equivalent heat-transfer coefficient and S, is the cooling surface of the winding, 


If we ignore heat dissipation from the end surfaces, we shall have 
S,=n(d,+d; (5) 


where dy and dj are the winding's inner and outer diameters. 


By eliminating the quantities ¢, w, r , and S, from the expressions given above, we obtain the direct 
relationship between the operating voltage, U,;, on the one hand and the design parameters of the relay and the 
overheating of its winding, on the other: ° 


Ua = Unk, Fly (6) 


my ) (20 (1+ a8) 


We recall that the quantity Uc, is the operating voltage at winding temperature 9,,,, that the relay of 
given design will actually have if the steady-state overheating of its winding, subject to a voltage K,Uy, is @. 
We now consider what the operating voltage must be for winding temperature Samp so that the relay will give 
failure-free operation in use. By the proper reasoning it is easily established [3] that this voltage, which we 


denote by equals 
U 
Here, Kg= a is the coefficient of possible decrease of the network voltage, K, is the coefficient 


of construction tolerance, characterizing the accuracy with which the given run of relays was prepared (tolerances 
on the number of turns, the impedance of the conductor, the residual magnetism, the friction of the moving parts, 


4 
etc.); for relays in ordinary use, Ky mw 1.2; Kg= _ is a coefficient which characterizes the possible 


increase in the opposing force due to the action on the relay of the accelerations to which it may be subject in use. 


* Here, the coefficient of core charging, fy, and the heat-transfer coefficient y, are assumed to be constant. 
In practice, the almost linear dependencies of the charging coefficient on the winding's density of impedance 
(i,e., on the impedance per unit volume of winding) and of the heat- transfer coefficient on temperature [2] 
only slightly complicate the expression without changing the qualitative nature of the relationship. Therefore, 
only for definitive precision of the computations should the nonconstancy of fy and yp be taken into account. 
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Obviously, then the condition for reliable relay operation will be Ucy > Ucg. 


By equating Expressions (6) and (7) we can find an interesting relationship between any of the winding 
dimensions and its steady-state overheating @ for given conditions of use. 


For the majority of dc relay types in use, a decrease in dimensions is easily made by shorting the core 
length 1,, since with this the forming of the contact groups is almost unaffected, but the resistance of the 
magnetic circuit is sharply decreased, which permits an increase in the tractive force on small gaps. After some 


transformations we obtain, for this case * 
L.=K V a+ (8) 


FPF 1, Sip 
8 mp ) iA, 


where 


K 


(9) 


The expression obtained shows graphically that an increase in working overheating does not lead to an 
indefinite decrease in winding dimensions, but tends to some finite value, Thus, the core length 1, of an 
electromagnet (the other dimensions being given, and invariant) has the limiting value 


lim = lim KV a. (10) 


If we now solve the same problem as it applies to the winding of an electrical element in the general 
case, where it is required to provide a given value of magnetic flux, ©, under given conditions of usage then, 
as is easily seen, only Maxwell's formula drops out of consideration, and we shall be led to the same results, 
In Expression (9), instead of the product FS;, we should in this case place the quantity 6", 


However, in some cases (for example, in the control winding of a magnetic amplifier) it may be necessary 
to decrease, not the length, but the height, h,, of the core, Then, since 


do= 2he+dj , 
and in conjunction with Relationship (5), we will have 
S, = 2nl(hi+d;). 


By taking this fact into account, we obtain, for the height of the core, the formula 


d 

/ +e (a+ 4), (11) 
where 
K' = (12) 
With this, 
—d 

hem = limke= +V (4) + Ka. (13) 


* The magnetic flux # corresponds to the maximum armature gap when almost all the circuit's reluctance is 
concentrated in the air-gap. It is therefore possible to consider R constant in a first approximation, Strictly 
speaking, decreasing the core length will affect (favorably, it is true, in the sense of computational reliability) 
the magnitude of the leakage coefficient. 
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By using the method presented, one may obtain corresponding expressions for other cases arising in practice, 
and make the necessary refinements. * 


It should be remembered that a shortening of the winding dimensions will be accompanied by two adverse 
circumstances, First, a significant change in the winding resistance in the process of establishing its overheating 
will entail a corresponding change in the relay parameters (operating and releasing voltages, operation time, 
etc,) which might adversely affect the operation of the device as a whole. The second circumstance is the 
additional energy requirement of the winding with decreased dimensions, in the sense that this winding will re- 
quire power equal to that of the original winding only in the state of establishing the heating, but in the state 
when current is not heating it, its impedance will be correspondingly less, However, both these disadvantages are 
less if there is the possibility of preheating the apparatus, 


If the character of the temperature dependency of the winding's conducting material is taken into account 
more accurately [4], Expression (4) must be given in the form 


(4a) 
= 


In practice, even for overheating of the order of 500° C, the correction introduced by taking the coefficient 
6 into account will, for copper, amount to some 10%, However, it is interesting to note that, in this case, with 
increasing overheating the dimension /, will not tend asymptotically to some limit, but will have some mini- 
mum value, 


In fact, instead of (8), we shall have the expression 


L=KY a+5+pe, 


whose minimum value, found by varying 6, will be 


le min= KVa+ 2V 8. 
However, this dimension occurs at the overheating: 


VB -V0.43-10-* 


while the melting point of copper is about 1000° C. 
Thus, for relays, only Expressions (8), (9), and (10) are of practical interest. 


As was shown above, the winding dimensions depend, not only on the size of the unit impedance of its 
conducting material, p, but also on the magnitude of the resistance temperature coefficient of this material, a. 
As follows from Expressions (8) and (9), the ratio of the dimensions of two windings with identical overheating, 
@, but made of different materials, will equal 


l 
on pi (1 + a8) 
pe (1 + aad) (14) 
Of the widespread conducting materials which can compete with copper, in the sense of lowering dimen- 


sions, only silver qualifies, since it has both a lower specific impedance and a lower resistance temperature 
coefficient, But this difference is very small, and allows the dimensions to be decreased only by a few percent. 


6 1500°C, 


* For arbitrary computed values, particularly those with design significance, the coefficients K or K" are easily 
found from known values of a, @, and Ic orh, and dj by means of (8) or (11). Maintaining of the magnitude 
of K and K* guarantees the previous degree of reliability of relay operation in the given conditions of use. 
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At the same time, the specific gravity of silver exceeds that of copper by almost 20%, i.e,, the weight of the 
winding would not be decreased by the substitution. The transition to aluminum (although its resistance tempera- 
ture coefficient is lower than that of copper) might give rise to an increase of dimensions by 20 to 30% or more, 
but permits the weight of the winding to be at least halved. However, the increase in winding dimensions 
ordinarily requires a corresponding increase in the dimensions of other parts of the electric element (for example, 
the magnetic circuit of a relay), as a result of which the over-all weight of the element is not reduced. 


SUMMARY 


1, The effect of the resistance temperature coefficient is such that an increase in the temperature of 
winding overheating causes its dimensions to tend rather quickly to some minimum limiting values. 


2. Due to the closeness of the resistance temperature coefficients of the widely used conducting materials 
(copper, silver and aluminum), from the point of view of decreasing winding dimensions, the upper limits of 
overheating most effective for them are practically identical and, as shown by a number of planned designs, 
equal approximately 200° C. A sharpening of this limit can be carried out on the basis of the relationships given 
in this paper, starting with the concrete conditions of each individual case, 


3. In view of the rapid development of heat-resistant types of insulators, it should be expected that in the 
future the possibilities of heat-resistant insulators may be used for decreasing dimensions and weights, in view 
of the circumstances mentioned in 1 and 2 above, The question therefore, arises as to the development of con- 
ducting materials with sufficiently small specific resistances and, moreover, small resistance temperature co- 
efficients, As follows from Relationship (14), these materials must possess, for high temperatures, lower specific 
resistances than existing conducting materials, The specific weights of these materials must be small. 
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INITIAL PRODUCTION RELIABILITY OF DEVICES 


Ya. A. Rips 


(Moscow) 


A method is suggested for the determination of the initial production 
reliability, the magnitude of which is determined by design peculiarities of 
the device and by production factors, Analytical expressions for computational 
purposes are deduced, and the relationship is shown between the initial production 
reliability and the rated safety factors. The effects of additional inspection and 
finishing are considered, and an engineering method for taking these into account 
is provided, this method approximating the distribution law by a series of normal 
distribution laws and their linear combinations, 


The reliability of any device is a function of its lifetime and the conditions under which it is used. It may 
be defined as the probability of failure- free operation of the device in given conditions of use during a fixed 


time, t , and is measured, as an average, as the probability of the relative quantity of devices operating without 
failure up to timet. 


Insofar as one speakes of relative quantities, one must specify with some exactitude which quantity is taken 
as a reference value for the failure-free operation of the devices considered. Sometimes, this reference value is 
taken to be the number of devices which are produced in working condition. This, however, excludes from con- 
sideration the failures attributable to devices which are inherently unable to operate, the fraction of which is de- 
termined by design peculiarities and peculiarities in the technological preparation of the devices, by the choice 
of materials, by the character and quality of the inspection given, and by other production factors, Failures of 
these devices occur in the first cycle of their operation, and their number is one of the factors (in the case of 
one-shot devices, one of the basic factors) which determine the reliability in use. Therefore, the index of re- 
liability, R(t), should be taken relative to the number functioning at time t out of all the devices produced in 
the productive process (i.e., put into use) independently of the time of their being taken out of service. 


The probability that a device will function correctly we define as the initial production reliability, R(0), 
the probable relative number of new devices working without failures at the time of their production (t=0). Such 
a definition of this concept is advantageous because it generalizes the concept of reliability, and provides the 
possibility of taking productive factors into account. 


We note that various production measures— control, inspection, finishing, etc. — can be considered as 
methods of increasing the initial production reliability. In particular, the continuous testing of devices before 
their delivery is a measure which increases the initial production reliability to unity. 


In addition to the fact that the reliability R (0) is one of the factors which determine the reliability of the 
device in use, its size is interesting in itself. In planning, it is advantageous to compare design variants from the 
point of view of their initial production reliability, which might be used as one of the criteria in evaluating 
devices for a given productive process, Asialysis of the expression for initial production reliability allows one to 
determine which components lower it to the greatest extent, thereby indicating at which steps in the productive 
process it is necessary to introduce additional finishing, inspection, etc, Starting from a desirable magnitude of 
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R(0), one might also calculate and rationally choose the corresponding production (rated) safety factors (tolerances) 
which, up to now, are usually determined only experimentally. 


In the present paper, we consider a general methodology for computing the initial production reliability, 
This methodology allows the quantity R(0) to be determined even in the design development stage, thus allowing ‘ 


a predetermination of the initial production reliability of a planned device even before its actual production be- 
gins, which is very important. 


Method for Computing the Initial Production Reliability 


We shall agree to understand by the failure of a device, not only its going out of order, but also such a 
state in which at least one of its properties or working characteristics is beyond some previously established limit, 


The reasons for the failure of a new device may be divided into two groups, The first of these consists of 
reasons having a random character, which give rise to breakdowns and malfunctions of devices under the influence 
of random actions in the productive process, It is impossible to predict these failures, Statistical data, obtained 
by the author with respect to approximately 70,000 relays and 90,000 contactors, showed that the relative frequency 
of such failures, with respect to each reason |, retained an approximately constant value, 1— y,. For a given 


productive process, these quantities may be considered as probabilities of failure of new devices due to reasons in 
the first group, 


In the second group fall those causes which give rise to deviations, beyond the established permissible limits 
of item dimensions, material properties and characteristics of the device's elements, Depending on the principle 
of operation and the device's design peculiarities, various of these reasons will have various degrees of influence. 


Due to the deviations of a number of parameters on the unfavorable side, the properties and working 
characteristics of a device may go beyond the admissible limits: the probability of failure-free operation of the 
device (considering only causes in the second group). In dependence on the construction and principle of opera- 
tion of the device, this circumstance can be expressed, for each characteristic, by the probability that an in- 
equality of the following type hold 


A>B A—B>O, (1) 


where A and B are parameters characterizing, respectively, the state of the characteristic and the admissible 
limit of its translation, 


Relationships of the type given in (1) may be called reliability criteria for the device's characteristics, 
For various characteristics, such criteria may easily be established from the physical idiosyncracies of the device's 
operation [1]. It should be stated that these criteria must be necessary and sufficient conditions that the operating 
characteristic lie within admissible limits, For example, the criterion for reliable operating and holding of an 
electromagnetic relay is the relationship Fe), mag > Fmech» etc. 


Taking what has already been presented into account, we may find the magnitude of the initial production 
reliability of a device, R(0), in a definite ratio, in the following way: 


R (0) = toP I(Ao— Bo) > 9), (2 


r 
where Yo = ut! —r, and the subscript 0 indicates that the parameters A and B are considered at time, 
1 


t=0, 


These parameters are continuous random variables, realized by some value or another in each given device, 
If the distribution laws are g(Ag) and y(B,)* and their parameters are uncorrelated, then the distribution law 
for the difference of the random variables, g(Ag— By), can be found from the relationship 


* The same symbols are used for the random variables and the current variables corresponding to them. 


789 
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Ba) = | 9(Bo) {(4o— Bo) + Bol smal 


men 
and the magnitude of the initial production reliability will equal 


R (0) = %o (Ap — Bo)d (Ag — Bo) =o ? (Bo) [((Ao—Bo)+ Bo] dBo| d (Ay — Bo). 


0 —c 


It is clear from Expression (4) that the computation of R(0) requires that one find the distribution laws of the ¢ 
the parameters A, and B,, these being functions of a number of independent random arguments, x, Xg, . . . , X,. Ax, 
Physically, these arguments are either parameters characterizing the geometric form of the device, or parameters meat 
of the items’ materials or characteristics of the individual elements of the device. Each of the arguments, in its 
turn, is a function of a large number of other independent random variables. The values of these quantities can 
be different within the limits established for their admissibility, which are determined by a large number of 
productive factors, the appearance and intensity of which depend on random circumstances and causes which may 


lie far outside the sphere of the phenomena considered. and 


In the sequel, it will be convenient to consider, as the random variables, the deviations from nominal values, 
For the nominal values, we take those values with respect to which the limits of admissibility will be symmetrical. 
Then 


Arn = In— tym, AAg= AB= By — By. 


In this case, each oj the random variables, Ax, Axy,..., Ax, being a function of a large number of 
other random variables, is formed by a summing operation. If we assume that the Lyapunoy conditions hold, i.e,, 
that none of the terms essentially dominates the others, and that the number of terms is sufficiently large, then 
we may consider that the distribution laws of the quantities Ax;, Ax, . . . , Ax, are sufficiently close to the 
normal law. 


The data given in[2, 3] attest to the fact that, in the majority of cases, the distribution laws of the random 
variables of the type considered here differ very little from the normal, since the normal distribution is character- 
istic of those production processes in which all measures to inspect for normal values are taken, and this is just 
the case in practice. This also was confirmed by the distribution laws obtained by us for the winding impedance, 
time and current of operating and releasing, armature travel and core diameter with respect to 500 specimens of 
the new relay type RKN of the “Krasnaya zarya® ("Red Dawn" ) factory. Moreover, in solving the problem 
of finding the initial production reliability of actually nonexistent designs, it is not possible to start from concrete 
conditions, nor to get an approximate fix from statistically obtained data, Therefore, it is necessary to start from 
general, highly probable, assumptions which will permit, with an accuracy sufficient for practical purposes, the 
problem posed to be solved. These assumptions are just that the random variables Ax;, Ax, ..., Ax, are dis- 
tributed normally, which are the initial assumptions which we shall make. 


in te! 


The normal distribution law is characterized by two parameters — the mathematical expectation M, and 
the mean-square deviation ox, The random variable Ax, is centered, and its mathematical expectation is zero, 
An approximate value of the mean-square deviation may be obtained from the fact that the area bounded by the 
normal distribution curve and the straight lines Ax) = + 3 ox, equals 0.9973, and an outcome with this prob- 
ability is practically certain, Therefore, it is possible to consider the absolute value of the field of the given 


4 
tolerance, , as equal to 6 ox, and ox, Then 


where 
(11), 


(Az,,)* 
= ———— exp (— ). (5) 
( ) 
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We now attempt to determine the distribution law of the parameter Ay. Since the magnitude of AA, is 
small in comparison with Ap, the latter may be linearized in the region of possible values of the random argu- 
ments, where the mean values of the arguments are best taken as their nominal values, xnp. 


Then, 


1 


It is obvious from Expression (6) that the random variable AAg is a linear function of the random variables 
AXxy, Axg, .. . , AX, and its distribution law can be found by means of well-known expressions for determining 
the distribution law (composition) of a sum of random variables [4, 5). If the distribution laws of the variables 
Axy, Axg, . . . , AX; are normal, then the distribution laws of the variables Ay and AA, will also be normal, with 


mean-square deviations equal to 
cA = —) (1) 
6 


and mathematical expectations equal, respectively, to AN, or zero, Thus, 


The distribution law for the parameter By is determined in an analogous fashion, 
By computing g(Ag— By) from Expression (3) and substituting it in (4), we obtain 


Bo) — (Any — 


ROO) = | ( Bob 


By means of the substitution oO nas et is z the integral in Formula (9) may be expressed 


Voix 
in terms of the tabulated Laplace function @ (x): 
V hte} 


This expression defines the magnitude o/ the initial production reliability of the device. 


A 
If we introduce the safety margin (tolerance), 8% = ag , into Formula (10), we obtain 


— ~ are the relative mean-square deviations. By means of Expression 


(11), one may determine the production tolerance necessary to guarantee a given initial production reliability. 


where = and oR * 
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Taking Additional Inspection and Finishing into Account 


It was mentioned earlier that various production measures may be considered as methods of increasing the 
initial production reliability. The necessity for taking these measures arises in those cases when the established 
tolerances do not guarantee the required characteristics of the device, but a tightening of the tolerances is either 
impractical (for technological reasons) or impossible (if standard materials, made to GOST tolerances, are used), 
In these cases, the components and element characteristics of the device are subject to inspection after their 
preparation, the admissible tolerance 6» here being set less that 60, (Fig. 1). The inspection point rejects that 
portion of the elements where the magnitude of Ax goes beyond 
the more rigid limits (this portion corresponds to the stripped 
area under the curve of Fig. 1), The rejected portion of the 
elements may either be simply discarded or may be subject 
to finishing (reworking), It may be assumed that the amount 
of effort expended in reworking is random and may, with equal 
probability, provide any value of Ax in the interval 55. There- 
26, hy fore, after reworking, the distribution law of the variable Ax 

will, for this fraction, be uniform in the interval 6». In many 
Fig. 1 cases, the magnitude of 5) « 60, , and practically all the 
elements are subject to reworking. 


Thus, as a result of additional inspection and finishing (rework), the distribution laws for some of the variables 
Axy, OX, . . . , AX, may differ from the normal law, being either truncated normal, a sum of truncated normal 
and uniform, or uniform, The distribution laws of AA, and AB, will then be compositions of a number of variables, 
parts of which have the distribution laws just cited, the rest of which will be normally distributed. 


The sought-for distribution laws for AA, and AB, (and, consequently, also for the parameters A, and By), 
necessary for computing the magnitude of R(0), can be produced by successive operations of determining k— 1 
compositions of centered random variables of the form 


These operations can be carried out most simply in the case when each composition is a normal law. 
However, not all these compositions are normal laws, and their successive determination by means of the general 
formulas for finding the distribution law for a sum of random variables entails great mathematical difficulty, 
which is a serious obstacle to practical application. With the 
aim of avoiding these difficulties, the attempt was made to 


2 find approximations for the composition of distribution laws 
aT most frequently arising in practice, This was successfully 

sd carried out with an accuracy sufficient for engineering applica- 
tions, 

ib VA The computations which served as a basis for these approxi- 

mations were very arduous and will be omitted here. Their 

I results are given in the table, which was supplemented by a 

e a number of simple cases well-known from probability theory 

a A {4, 5]. The table includes the majority of compositions which 
mal A a a occur in practice, and can be used to advantage in computing 


the initial production reliability of devices, 


mg.'S The coefficients k, and kg, functions respectively of the 


parameters m, and mg, are defined by the graph of Fig. 2, For the computations, it is convenient to use the di- 
mensionless relative quantities, of the form 
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After the determination of the distribution laws for the parameters Ay and By, the computation of the initial 
production reliability is carried out in accordance with Formulas (4), (10), or (11). 


For example, let us consider the case when the distribution laws for parameters Ay and By are compositions 
of random variables distributed normally and of variables distributed by a law which is the sum of a truncated 
norma! distribution and a uniform distribution (table, line 8): 


9 (Ay) = Hyn(Apg, Ayo, Sex) + (1— Hy)n(Ao, Ano, 201), 
(Bo) = Hyn(Bo, By, + (1 — H2)n(By, Rio, 202), 


H, =1—0(5"), H,=1—0(32 


2c.) ° 


After substituting @(A,) and ¢(B,) in (4) and integrating, we find that 


RO) = + + | 
+H +o + 
De 


[1 + o( 
D1 D2 


In this case, the expression defining the magnitude of the initial production reliability consists of four terms, 


- By means of the method presented here, the effect on the device's initial production reliability of the size 
of the established tolerances, the quality of the materials used, the technological means and other productive 
factors may be easily determined, 


As an example, we adduce here the results of computing 


& ae) Ee the magnitude of the initial production reliability of operating 
ee an electromagnetic relay of type RKN with three contact groups, 
C SN The example is only of an illustrative nature, since approximate 
Py case “b* numerical values of the parameters were used, 
if The computation was made for two casesy a) the relay 
9 Pin, coils, after preparation, were inspected for impedance value, 
tN and those coils which were beyond the established tolerance on 
ee number of ampere-turns were rejected; b) the relay coils were 
™! not subject to this inspection. It was assumed, moreover, that 
nv an in the preparation of the relays, the air-gap was inspected and 
va a7 7) finished, and the contact springs were regulated, 


The computed results are shown in Fig. 3, which gives 
Fig. 3 the dependence of the initial production unreliability, Q,)= 
= 1—R(@), on the selected tolerance, s», for the cases con- 
sidered, An analysis of the curves shows the great influence of the additional inspection on the quality and re- 
liability of the end-product. For example, for a tolerance of sy = 1.5, of each 100,000 relays produced, the 
probable number of unsatisfactory relays was 6 in the first case and increased to 135 in the second. 


The author wishes to express his gratitude to Professor B, C. Sotskov for his many good counsels, which were 
taken to heart, 


Contr 


Engli 


Contr 


N. M, 


Auto! 


where 
* See 
196 
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SEVERAL SIMPLIFICATIONS IN THE APPLICATION OF THE D-PLOT METHOD 
FOR DETERMINATION OF CRITICAL VALUES OF A REAL PARAMETER 
WHICH OCCURS LINEARLY IN THE CHARACTERISTIC EQUATION 


Vv. A. Andreyuk 


(Leningrad) 


Several simplifications are presented in the application of the D-plot 
method for the determination, according to stability criteria, of critical values 
of a parameter which enters linearly into the characteristic equation. 


We presume that it is required to find the critical values of a real parameter k which enters linearly into 
the characteristic equation of a system which we will write as 


P(A) + kQ(A) = 0, (1) 


where P(,) and Q(\) are polynomials with real coefficients of arbitrary degree in ). 


We will call those values of the parameter k critical values which map a region of variation of the param- 
eter k within subregions in which Eq, (1) has a definite number of roots with negative real parts. 


To solve this problem, it is natural to make use of the method of constructing D- regions [1, 2]. 
According to this method, it is necessary to carry out the following operations: 

a) substitute into Eq. (1) = jw such that P (jw) = Py (w*) + (ws and Q(jw) =Qy + (ws 
b) multiply the result of the substitution by the complex conjugate of Q(jw)s 

c) separate out of this product the real and imaginary parts and equate them to zeros 


(— + = 0, (2) 
PQ: + + k (Q? + w7Q?) = 0. (3) 


The critical values of the frequency are determined from Eq. (2) while the corresponding values of the 
parameter k are determined from (3): 


The critical value of k given by 
(5) 


198 


cor 
(w 
firs 
to se 


corresponds to the null root of (2) (w =0). 


It is, however, more convenient to evaluate those values of k which correspond to nonzero roots of Eq, (2) 
(w #0) by two equations which are simpler than Eq, (4). To get these equations we multiply Eqs. (2) and (3) 
first by wQ, and Q,, respectively, then again by Q, and wQs,, respectively, and add each time the results. 

As a result of this, we find 


+ °Q3) (Pi + Qi) = 0, @ (Qi + (Pa + = 0. 


Since Q? + w*Q} #0 and the roots of Eq. (2) are nonzero (w #0), then the factors Py+kQ, and Py+kQ, 
are equal to zero, Hence, we obtain 


(6) 
(7) 


These expressions have a much simpler form in comparison with Eq, (4). In concrete problems, it remains 
to select the simplest of the two equations (6) and (7). 
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ON THE CALCULATION OF THERMISTOR CHARACTERISTICS 


N. P. Udalov 


(Moscow) 


An analytical method, which follows from equations of thermal equilibrium, 
is given for calculating the volt-ampere characteristics of thermistors. 


A series of articles [1-5] is devoted to the development of methods of constructing static thermistor volt- 
ampere characteristics for any given temperature. However, the problem cannot be considered to be completely 
solved, since the proposed methods either are too cumbersome or are based on the use of a large amount of data 
which must be obtained experimentally for each thermistor. 


A method is given below which makes it possible to compute, in a very simple manner, the volt- ampere 
characteristics of thermistors, on the basis of their temperature characteristics. 


For arbitrary points of a static volt-ampere characteristic curve the equation of thermal equilibrium holds: 


PR,= b0, 


where I is the steady-state current through the thermistor, R, is the resistance of the active element in the given 
steady state, b is the diffusion coefficient, and 6 is the overheat temperature (the difference between the 
temperature, @, is of the active element in the given steady state and the temperature, 6, of the surrounding 
medium). 


The diffusion coefficient b which enters into the equation takes into consideration all forms of heat pro- 
pagation from the active element of the thermistor (thermal conductivity, convection, and thermal radiation). 
Its magnitude is a function of the material, dimensions, and state of the surface of the active element of the 
thermistor and of the current-carrying parts, as well as of the properties of the surrounding medium. By comput- 
ing the values of the diffusion coefficient at a series of points of a family of volt-ampere characteristic curves, 

it is possible to satisfy oneself that it is a function of the overheat temperature [4]. Proceeding from the assump- 
tion that the deviations of the dimensions from the rated values which take place in serial and mass production 
of thermistors may not have a significant influence on the value of the diffusion coefficient, we will consider 

the curve b= f (@), which is constructed for some particular sample thermistor, as valid for any thermistor of the 
same type. Consequently, by setting up this functional dependence and the thermal characteristics of a given 
thermistor, it is possible to compute its volt-ampere characteristics for any given temperature which does not 
exceed the maximum possible for the given thermistor type. | 


We will assume that it is necessary to find points of the volt-ampere characteristic curve for some tem- 
perature 6, of the surrounding medium. It is known that to each point of the static volt-ampere characteristic 
curve there corresponds a value of the steady-state temperature of the active element, @,, which may take all 
values in the range 0, to 6,4, where @y is the maximum possible temperature for the given thermistor, Having 
been given the necessary number of values 9, in the given interval, it is possible to find by means of the 
temperature characteristics the corresponding values Rg and by means of the curve b= f (@) it is possible to 
find the diffusion coefficient, calculating in advance the overheat temperature 6 =, — 0. 
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Substituting the values Rg, @, and b which were thus found into Eq. (1), it is easy to calculate the current 
value: I=, | ~ , and then the voltage drop across the thermistor: U = IR,. 

For example, suppose it is desired to construct the volt- ampere characteristics of a thermistor of type MMT-4 
with temperature of the surrounding air 6,=16°C. The temperature characteristics of the given thermistor 


sample are described by the equation R, = 0,458 e , 4 model curve b= f (6) is given in Fig. 1, 
We will find on the volt-ampere characteristic of 


of thermistors. 


Using this data, we find the desired 


= 7.95 v. 
uv 
60 


0 2 4 i ma 


Fig. 3. Volt-ampere characteristics of 
various types of thermistors, calculated 
for an averaged curve b= f (@). 


the thermistor the coordinates of the point K which BL 
corresponds to temperature = 80° C of the active 6°16" 
° 
mw/* C ° 
| 
0 
\\ 
Fig. 1. Curves which show the diffusion co- Fig. 2. A family of volt-ampere characteristics 
efficient as a function of the overheat of the type MMT-4 thermistor which was cal- 
temperature for several widely used types culated by means of the temperature character- 


element. The resistance of the active element is Rg =0.410 kohm, the overheat temperature is 6 =64° C, and 
the diffusion coefficient b=6.6 mw/* C (which we find on the curve of Fig. 1). 


If the volt-ampere characteristic is computed for @» = 60° C, then for the same temperature 6, = 80° C 
of the active element, 6 = 20°C, b=7.7 mw/*C and the coordinates of the point K* will be 1=19.4 ma, U= 


istic and curve b= f (6). 


coordinates of the point K; 1=32 ma, U =13.1 v. 


In Fig. 2 three volt-ampere curves are shown which were 
constructed by the method described above together ~ith plotted 
points which were found experimentally, The agreement of the \ 
calculated results with the experimental appeared to be quite 
satisfactory: for@ > 10° the relative error did not exceed 5 %. 
It remains to observe that in the construction of the volt-ampere 
characteristics of the given thermistor, the curve b= f (@) was 
used which was constructed on the basis of a family of volt-ampere 
characteristics of another sample thermistor of the same type which ; 


has a considerably different temperature characteristic (R,=2.9e T . 
In this fashion, the proposition stated above about the universality 


of the curve b= f (6) for a definite type of thermistor was | 
confirmed. 


/ 
| 
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The possibility also was verified of using the curve b= f (@) for a thermistor with the same geometrical 
dimensions, but made out of different materials, Types MMT-1 and KMT-1 thermistors were investigated, 
Their diffusion coefficients for identical overheat temperatures differed by approximately 20%. A curve for 
the average value of the diffusion coefficients as a function of the overheat temperature given in Fig. 1 was con- 
structed which was used in the formation of the volt-ampere characteristics ior KMT-1 having temperature 


4100 /T 2380/T 
characteristic R,=0,255e and MMT-1 with temperature characteristic R, = 3.56 € . The results of 


the construction and the corresponding experimental points are presented in Fig. 3. 
In all examples given, the temperature characteristics were computed at two points [6]. 


SUMMARY 
1. For the calculation of the static volt-ampere thermistor characteristics it is sufficient to have a tem- 


perature characteristic of the given sample and a curve giving the diffusion coefficient as a function of the over- 
heat temperature. 


2. The relation b= f(@) may be found by calculating values of b = at points of the volt- ampere 
characteristic taken for known constant temperature of the surroundings for any thermistor of the given type. 

3, In thermistor handbooks the curves b= f (@) or their analytical expressions should be given. 

M, M, Sverdlov and E. V. Zaitsev participated in the experimental part of the work presented here. 
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ON THE QUESTION OF ESTIMATING THE SWITCHING CAPACITY 
CONDUCTORS FOR DC CIRCUITS 


B. S. Sotskov 


(Moscow) 


The parts of a switch most subject to wear are its contacts. 


We will consider just what takes place upon breaking the contacts in a circuit with voltage E, in which a 
current I is flowing. 


Let U, be the minimal voltage and I, the minimal current necessary for arcing. IfE > U, and I> ky, 


then upon breaking the circuit there appears both a liquid metallic bridge which effects a transfer of metal from 
anode to cathode and an arc which causes wear of the contacts. 


The energy given out at the cathode during the arcing process is equal to 


te 


A= Upc tdt = idt=Uo.9, 
0 0 


where Upc is the cathode voltage drop in the arc, Ugce = Up, i is the current in the arc, t- the time during 


which the switch is closed, and q the quantity of electricity flowing through the arc. The extent of the wearing 
of the contact is proportional to the cathode energy output: 


where y is the erosion coefficient, 
The current i in the contact circuit is determined by the equation 


di 


where U, is the voltage drop across the arc, R is the resistance, and L is the inductance of the switching circuit. 


The voltage drop across the arc is a function of the time, U, = f (t); on the basis of experimental data, it is 
possible to assert that 


Ug = Ugg + hit + 


where Uy is the initial voltage drop across the arc, Frequently, U, may be considered to be approximately 
equal to Ugg Ug Ugg. Applying the initial conditions t =0, i=1= E/R, we obtain 


ta 


G=kA=hUg I= 14 
wh — 
ere r= R | 
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The value of the closure time t,. may be determined from the condition i =I, for t =t, in the form 


“a0 
In the case of contactors, it commonly happens that 1k « E and it is possible to assert that 
U, 
=t 
The quantity of electricity which is flowing through the arc after time tg is given by of 
3) 
Consequently, the wear of the contacts after N, breaks may be determined by the expression = 
N N 
G =G'N, = N, BI = 
2 (Yao 2) (Us 
the; 
Knowing the material of the contact (and consequently, the value of y, the maximum amount of wear var 
of the contact (G = G,) and the value of Ugg, which depends upon the construction of the contact, and the param- 
eters of the connected circuit (£, I, r), it is possible to determine the service life of the contact (i.e., the 
maximum number of breaks) with wear, occurring only at the 
WISI time of breaking as 
N + 
Ic 5c 
01 
I 2 U E 
N N Go 
‘ 
= 
contact materi Lp 
Oga 
If vibration of the contacts occurs when the circuit is 
closed and it arcs, then the wear of the contacts at the make 
time is defined as 
= al’ tay 
where I’ is the current at the moment of making, ky = = is the ratio of the time of arcing tg to the vibra- and | 
v 


tion time tay, kp = 0,3-0.7. The maximum number of breaks is N, = S- . In the figure the relationship 


N= ora is illustrated for various values of Gp. Having determined the amount of electricity which is flowing 


through the arc upon opening the circuit, qy » = ui? — » we find the value N, on a curve on the figure. 
ao 

Having determined the amount of electricity which is flowing through the arc upon closing the circuit, qz = I"tayky, 

we find the value of N,, The actual service life of the contacts Nx, taking into consideration the wear of making 


and breaking upon them, may be defined by the equation 


804 


N, +0") = (Ft = Co, 


and consequently, 
4 i 


The quantity ™ = Nq = K may be called the switching capacity of the contact. With good construction 


of contactors and magnetic starters, the vibration time t,, is small; hence, G" < G' and 


oly E\ U E 2 
3) ao” 3 


1 
N 


The formulas obtained for the service life and switching capacity are, naturally, approximate. However, 
they allow a first approximation evaluation of these two basic characteristics of contact apparatus, as well as : 
varying the service life as a function of the working state of the contacts (i.e., of the values I and r). 
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ERRATA 
(Automation and Remote Control, No, 2, 1959) 


In *To the theory of y-telay devices with gas discharge counter,” by Vasiley and Klempner, for Eqs. (19") 
and (20°), read 


Tine Ng) +2Y + + > tM, (19°) 


2? — V + At (nm, + nq) 2* 
Ne == Mg +- 


In bibliography on switching theory, under Zemanek;, H.; should read: 
Arch, elektr, Ubertrag., 12,1 (1958). 
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SIGNIFICANCE OF ABBREVIATIONS MOST FREQUENTLY 
ENCOUNTERED IN SOVIET TECHNICAL PERIODICALS 


AN SSSR 
FIAN 

GITI 

GITTL 

GOI 

GONTI 
Gosénergoizdat 
Gosfizkhimizdat 
Goskhimizdat 
GOST 
Gostekhizdat 
GTTI 

IAT 

IF KhI 

IFP 

IL 

IPF 

IPM 

IREA 

ISN (Izd. Sov. Nauk) 
[Yap 


SPB 

Stroiizdat 
URALFTI 
TsNHTMASh 
VNIIM 


Academy of Sciences, USSR 

Physics Institute, Academy of Sciences USSR 
State Scientific and Technical Press 

State Press for Technical and Theoretical Literature 
State Optical Institute 

State United Scientific and Technical Press 
State Power Press 

State Physical Chemistry Press 

State Chemistry Press 

All-Union State Standard 

State Technical Press 

State Technical and Theoretical Press 

Institute of Automation and Remote Control 
Institute of Physical Chemistry Research 
Institute of Physical Problems 

Foreign Literature Press 

Institute of Applied Physics 

Institute of Applied Mathematics 

Institute of Chemical Reagents 

Soviet Science Press 

Institute of Nuclear Studies 

Press (publishing house) 

Leningrad Electrotechnical Institute 

Leningrad Institute of Physics and Technology 
Leningrad Institute of Metals 

Leningrad Institute of Precision Instruments and Optics 
State Scientific-T echnical Press for Machine Construction Literature 
Moscow State University 

Metallurgy Press 

Moscow Regional Pedagogical Institute 
Scientific Research Association for Physics 
Scientific Research Institute of Physics 


_ Scientific Research Institute of Mathematics and Mechanics 


Scientific Institute of Motion Picture Photography 
People’s Commissariat of the Heavy Machinery Industry 
State Press of the Defense Industry 

Joint Institute of Nuclear Studies 

United Scientific and Technical Press 

Diwision of Technical Information 

Division of Technical Science 

Radium Institute, Academy of Sciences of the USSR 
All-Union Special Planning Office 

Construction Press 

Ural Institute of Physics and Technology 

Central Scientific Research Institute of Technology and Machinery 
All-Union Scientific Research Institute of Metrology 


NOTE: Abbreviations not on this list and not explained in the translation have been transliterated, no 
further information about their significance being available to us — Publisher. 


Izd 
LETI 
LFTI | 
LIM 
LITMiO 
Mashgiz 
MGU 
Metallurgizdat 
MOPI 
NIAFIZ 
NIFI 
NIIMM 
NIKFI 
NKTM 
Obrongiz 
olyal 
ONTI 
OTI 
OTN 
RIAN 


Publication of a “Soviet Instrumentation and Control Translation Series” by the Instru- 
ment Society of America has been made possible by a grant in aid from the National 
Science Foundation, with additional assistance from the National Bureau of Standards 
for the journal Measurement Techniques. 


Subscription rates have been set at modest levels to permit widest possible distribution 
of these translated journals. 


The Series now includes four important Soviet instrumentation and control journals, 
The journals included in the Series, and te subscription rates for the translations, are 
as follows: 


MEASUREMENT TECHNIQUES — Izmeritel’naya Tekhnika 
Russian original published by the Committee of Per year (12 issues) starting with 1959, No.1 


Standards, Measures, and Measuring Instruments General: United States and Canada . . $25.00 
of the Council of Ministers, USSR. The articles in Elsewhere .. - 28.00 
this journal are of interest to all who are engaged Libraries of nonprofit academic institutions: 

in the study and application of fundamental meas- United States and Canada. . $12.50 
urements. Both 1958 (bimonthly) and 1959 15.50 


(monthly) available. 


INSTRUMENTS AND EXPERIMENTAL TECHNIQUES 
Pribory i Tekhnika Eksperimenta 


Russian original published by the Academy of Per year (6 issues) starting with 1959, No. 1 
Sciences, USSR. The articles in this journal relate General: United States and Canada. . $25.00 jj 


to the function, construction, application, and op- Elsewhere ..... 

eration of instruments in various fields of experi- Libraries of nonprofit academic institutions: 

mentation. 1958 and 1959 issues available. United States and Canada. . $12.50 
Elsewhere ... 


- AUTOMATION AND REMOTE CONTROL — Avtomatika i Telemekhanika | 


Russian original published by the Institute of Per year (12 issues) starting with Vol. 20, No.1 5 
Automation and Remote Control of the Academy General: United States and Canada. . a 


of Sciences, USSR. The articles are concerned Elsewhere .. . 88005 
with analysis of all phases of automatic control Libraries of nonprofit academic institutions: 
theory and techniques. 1957, 1958,and 1959 issues United States and Canada. . $17.00 
available. Elsewhere ...... -. 20.50 


INDUSTRIAL LABORATORY — Zavodskaya Laboratoriya 


Russian original published by the Ministry of Per year (12 issues) starting with Vol. 25, No.1 
Light Metals, USSR. The articles in this journal General: United States and Canada. . $36.00 5 


relate to instrumentation for analytical chemistry Elsewhere .. . B880F 
and to physical and mechanical methods of mate- Libraries of nonprofit academic institutions: . 
rials research and testing. 1958 and 1959 issues United States and Canada. . $1750 9 
Single issues of all four journals, toeveryone,each ... $60 


' Prices on 1957, 1958 issues available upon request 


SPECIAL SUBSCRIPTION OFFER: 


One year’s subscription to all four journals of the 1959 Series, as above listed: 


General: United States and Canada. . $1100) 
... . 12200 

Libraries of nonprofit academic institutions: ; 
United States and Canada . 
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